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i.  imnucrioM 

A.  BACKGROUND 

The  development  of  electronic  integrated  circuits  and  their  appli- 
cation to  our  defense  systems  have  led  to  increased  reliability  in  life; 
reduced  cost,  size  and  weight;  and  reduced  primary  power  requirements. 

The  example  set  by  this  development  has  provided  the  justification  for 
an  increasing  effort  addressed  to  an  integrated  optical  circuit  technology 
with  the  expectation  of  realizing  some  of  the  same  benefits  for  coherent 
optical  signal  processing,  sensors  and  communication  systems  while  simul- 
taneously providing  an  increased  information-handling  capacity. 

One  important  application  of  integrated  optics  is  in  the  field  of 

coherent  optical  signal  processing  executing  the  Fourier  transform. '' s 

Optical  signal  processing  is  admirably  suited  to  perform  linear  analog 

operations  such  as  matrix  products,  the  Fourier  transform  integration 
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and  related  correlation  and  convolution.  The  process  employs  a coher- 
ent source,  modulation  function  (s)  and  the  detection  function  to  be 
performed  throughout  the  signal  spatial  field  with  one  or  more  Fourier 
transform  lenses  between.  The  optical  signal  function  is  usually  per- 
formed for  a two-dimensional  image,  whereas,  an  integrated  optics  for- 
mat will  limit  the  processing  to  a single  transverse  dimension  as  a 
function  of  time.  For  an  application  such  as  spectral  analysis  of  micro- 
wave  signals  extending  over  a broad  expanse  of  the  spectrum,  this  one- 
dimensional characteristic  is  not  a limitation. 

The  optical  Fourier  transform  applied  to  spectral  analysis  of 
microwave  signals  is  attractive  because  it  provides  a real-time,  wide- 
open,  continuous  response.  The  schematic  diagram  of  Figure  1 illustrates 
the  essential  components  of  an  optical  Fourier  transform  process  applied 
to  spectral  analysis.  Implementation  in  an  integrated  optics  format  of 
the  spectrum  analyzer  will  contain  a coherent  source  (presumably,  a DFB 
(GaAl)As  laser)  coupled  to  a thin-film  optical  waveguide,  an  assortment 
of  thin-film  waveguide  lenses0,,£erforming  the  beam  expansion  and  the 
transform  functions  together  with  a spatial  modulator  (presumably,  a 
surface  elasto-optic  Bragg-type  interaction,2within  a common  thin  film) 
and  a waveguide-coupled  photodetector  array  responsive  to  the  laser  with 
a temporal  bandwidth  encompassing  the  down-converted  microwave  signal. 
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Figure  1.  Coherent  Optical  Fourier  Transform  Signal  Processing  Schematic 

The  microwave  signal  is  temporarily  "stored"  as  a waveform  in  the  spatial 
modulator  which  is  simultaneously  accessed  by  the  optical  field  and  trans- 
formed into  its  spectral  components  spatially  separated  by  the  lens  and 
detector  in  an  array.  Use  of  an  integrated  optics  format  for  the  spec- 
trum analyzer  on  a single  substrate  provides  the  necessary  mechanical 
structure  to  obtain  rigidity  between  the  various  optical  elements.  It 
provides  the  potential  for  a substantial  reduction  in  size  and  weight  of 
optical  signal  processing  systems  and  provides  the  potential  for  a sub- 
stantial reduction  of  cost,  when  the  benefits  of  electronic  integration 
are  exploited  by  optical  integration. 

The  optical  Fourier  transform  process  requires  a series  of  lenses 
with  sufficient  aperture  and  resolution  to  encompass  the  time-spatial- 
bandwidth-product  system  requirements.  Near  diffraction-limited  lens 
performance  and  low  scattering  are  necessary  to  obtain  a useful  resolu- 
tion and  dynamic  range.  Lens  systems  designed  to  fulfill  these  require- 
ments using  conventional  optical  elements  in  three  dimensions  is  a highly 
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developed  art.  The  corresponding  implementation  of  lenses  in  two-dimensional 
thin-film  waveguides  constitutes  a developing  technology  of  which  this  re- 
port constitutes  the  third-year's  effort  sponsored  by  the  Air  Force  Avionics 
Laboratory.  The  first-year's  effortMaddressed  various  approaches  to  lenses, 
lens  systems  and  reflectors  as  may  be  required  for  implementation  of  the 
optical  Fourier  transform  function  in  an  integrated  optics  format.  The 
optical  waveguide  Luneburg  lens  and  geodesic  analogy  received  most  of  the 
attention.  The  second-year's  effort’addressed  the  waveguide  Luneburg  lens 
implemented  on  a thermally-oxidized  silicon  substrate  and  developed  a com- 
putational procedure  to  predict  its  performance  as  an  inhomogeneous  re- 
fractor. A portion  of  this  second-year's  effort  addressed  a companion 

component  of  the  spectrum  analyzer,  specifically,  waveguide-coupled  photo- 
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detector  array  with  charge-coupled  device  readout  which  was  carried  out  at 
the  University  of  Cincinnati.  The  third-year's  effort,  reported  herein, 
has  continued  to  address  waveguide  Luneburg  lenses  on  thermally-oxidized 
silicon  substrate  successfully  achieving  diffraction-limited  performance, 
the  reduction  of  waveguide  and  lens  scattering,  the  implementation  of  a 
spatial  filter,  implementation  of  a channel  waveguide  fan-out  structure 
and  improved  waveguide  reflectors. 

B.  PROGRAM  OBJECTIVE 

Basic  to  the  development  of  an  integrated  optical  technology  are 
the  requirements  for  (1)  a coherent  source,  (2)  the  general  modulation 
function,  both  temporal  and  spatial,  (3)  the  detection  function  and  (4) 
a waveguide  technology  to  appropriately  interconnect  all  of  these  com- 
ponents. All  of  these  components  must  be  implemented  on  a planar  struc- 
ture or  hybrid  coupled  together  using  a compatible  processing  technology. 

Optical  lens  elements  within  thin-film  waveguide  medium  have 
several  additional  design  freedoms  and  some  restrictions  when  compared 
with  conventional  optical  systems.  First,  the  lens  system  is  entirely 
immersed  in  the  thin-film  waveguide  which  has  an  effective  refractive 
index  determined  by  the  choice  of  materials  and  film  thickness.  The 
ambient  waveguide  effective  refractive  index  may  be  greater  than,  or 
less  than,  the  index  employed  for  the  lenses.  The  available  range  of 
refractive  indices  is  restricted  to  values  greater  than  that  of  the  sub- 
strate bulk  refractive  index  to  preserve  wave  binding  within  the  thin- 
film  waveguide.  Second,  the  use  of  nonplanar  substrates  permits  the 
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formation  of  lenses  of  uniform  refractive  index  by  changes  of  the  path 
length  because  the  rays  traverse  geodesics.  The  maximum  refractive  in- 
dex is  usually  limited  by  the  thin-film  waveguide  material  except  for 
geodesics  where  the  path  length  increase  synthesizes  the  larger  value. 
Third,  inhomogeneous  lens  elements  are  feasible  using  homogeneous  di- 
electrics and  may  be  produced  by  graded  waveguide  film  thickness.  Con- 
sideration of  these  various  ramifications  was  included  in  the  initial 
effort.  Optical  Waveguide  Lenses,  I.  One  important  conclusion  derived 
from  this  initial  effort  suggested  that  generalized  thin-film  Luneburg 
lenses  were  a favored  approach  for  optimization  in  integrated  optics 
format. 

Referring  to  Figure  1,  the  essential  lens  elements  include  a 
beam  expansion  system  and  a Fourier  transform  lens.  In  addition  to  the 
requirement  for  near  diffraction-limited  performance  and  low  scattering, 
the  design  of  these  lenses  must  consider  the  length  of  the  optical  axis 
because  of  the  Fourier  transform  "chip"  substrate  area  limitations. 
Reasonable  substrate  dimensions  are  comparable  to  that  employed  in  the 
integrated  electronic  circuit  industry.  Therefore,  consideration  of  the 
telephoto  effect  so  as  to  reduce  the  physical  focal  length  of  the  Fourier 
transform  lens  received  attention  during  the  initial  program.  Optical 
Waveguide  Lenses,  I.'4  The  required  large  compression  of  the  physical 
focal  length,  field  of  view  and  field  curvature  negated  further  consider- 
ation of  the  telephoto  effect  in  preference  to  a folded  optical  system. 

It  was  also  necessary  to  develop  plane  refractors  for  folding  so  as  to 
minimize  the  substrate  area  employed. 

The  second  effort,  Optical  Waveguide  Lenses,  II1,5  addressed  in- 
homogeneous thin-film  waveguide  lenses,  lens  systems  and  refractors  for 
the  waveguide  material -structure  and  exploited  thermally-oxidized  silicon 
as  an  isolation  layer  between  the  waveguide  and  supporting  silicon  sub- 
strate. Typical  Luneburg  lenses  employed  Ta205  as  an  overlayer  or  under- 
layer in  conjunction  with  Corning  7059  glass  employed  as  the  waveguide. 
The  thermally  oxidized  silicon  substrate  with  appropriate,  careful  pro- 
cessing has  proven  to  yield  the  lowest  scattering  substrate  thus  far 
identified. 
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Interest  in  the  thermally-oxidized  silicon  substrate  was  moti- 
vated by  the  potential  advantage  of  integration  of  both  optical  and 
electronic  components.  This  structure  is  of  particular  interest  because 
the  removal  of  the  isolation  oxide  provides  for  waveguide  coupling  to 
a photodetector  array  in  the  substrate.  The  requirement  for  the  spatial 
modulator  for  an  optical  Fourier  transform  "chip"  can  be  satisfied  with 
a piezoelectric  ZnO  overlay  film  to  launch  Rayleigh  waves  which  may 
Bragg  diffract  in  the  optical  waveguide. 

The  basic  refractive  element  considered  is  the  positive  generalized 
Luneburg  lens  (with  arbitrary  optical  focii  equal  to  or  greater  than  the 
lens  diameter).  Two  of  these  lenses  with  comparative  short  focal  length 
were  employed  to  demonstrate  a beam  expansion  system.  A single  lens  with 
a longer  focal  length  was  folded  using  reflectors.  This  second  effort 
included:  (1)  optimization  of  Luneburg  lens  radial  refractive  index 
profile  to  more  closely  approach  diffraction-limited  performance,  (2) 
reduction  of  scattering  losses  due  to  waveguide  films,  substrate  isola- 
tion layer  and  optical  element  boundaries,  (3)  establishment  of  lens 
reproducibility,  alignment  and  positioning  without  regard  to  the  specific 
profiles  and  (4)  suppression  of  mode  conversion  within  optical  elements. 

A part  of  this  second  program  was  directed  to  a waveguide-coupled 
photodetector  array  for  use  in  the  spatial  frequency  plane  of  the  opti- 
cal Fourier  transform  "chip."  This  waveguide-coupled  photodetector  array 
investigation  was  carried  out  by  the  Department  of  Electrical  and  Computer 
Engineering  at  the  University  of  Cincinnati  under  a subcontract  under  the 
direction  of  Professor  J.  T.  Boyd. 

C.  TECHNICAL  APPROACH 

Based  upon  the  results  of  the  preceding  programs.  Optical  Waveguide 
Lenses,  I,  and  Optical  Waveguide  Lenses,  II,  this  effort.  Waveguide  Optics 
for  Coherent  Optical  Processing,  reported  herein,  has  addressed  generalized 
waveguide  Luneburg  lenses  that  yield  diffraction-limited  performance. 
Corning  7059  waveguide  and  Ta20,-  overlay  lenses  were  employed  using  a 
thermally-oxidized  silicon  substrate.  The  thermally-oxidized  silicon 
provides  a smooth,  low-loss  isolation  layer  between  the  optical  waveguide 
and  the  silicon  substrate  to  prevent  absorption  losses  of  the  waveguide 
evanescent  field  in  the  substrate.  Generalized  waveguide  Luneburg  lenses 
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are  circular  symmetric  inhomogeneous  refractors  that  can  be  void  of  all 
aberrations  except  circular  field.  These  waveguide  Luneburg  lenses  are 
formed  by  deposition  through  a circular  mask  where  the  edge  of  the  mask 
is  shaped  to  control  the  waveguide  Luneburg  lens  radial  refractive  in- 
dex profile.  To  obtain  diffraction-limited  performance  requires  that  the 
radial  refractive  index  profile  fit  a specific  function.  A significant 
achievement  of  the  effort  reported  herein  has  been  the  synthesis  of  the 
mask  edge  shape  so  as  to  control  the  lens  radial  refractive  index  profile 
to  yield  essentially  diffraction-limited  results  within  the  capability 
of  its  measurement.  Other  problems  addressed  in  this  program  have  in- 
cluded waveguide  Luneburg  lenses  with  both  short  and  long  focal  lengths, 
namely,  F/4  and  F/25.  Further,  the  effort  has  investigated  various  means 
to  reduce  lens  and  waveguide  scattering  and  has  provided  means  to  measure 
both  the  in-plane  and  orthogonal  out-of-plane  scattering  transverse  to 
and  longitudinal  with  the  optical  beam.  A tolerance  analysis  for  wave- 
guide Luneburg  lenses  and  the  effect  of  lens  focal  length  correction 
has  been  included. 

Because  the  integrated  optical  spectrum  analyzer  must  conserve 
substrate  area  and  because  some  designs  require  a long  focal  length, 
means  to  realize  improved  reflectors  were  considered.  Cleavage  of  the 
single  crystal  silicon  substrate  through  the  isolation  layer  and  wave- 
guide has  yielded  good  reflection  employing  the  total  reflection  phenomenon. 

The  integrated  optical  spectrum  analyzer  requires  a high  degree  of 
transverse  mode  purity.  Two  types  of  spatial  filters  to  obtain  mode 
purity  were  investigated  and  demonstrated  as  an  alternative  to  the  re- 
quirement of  the  high  degree  of  mode  purity  from  the  mode  source.  Ab- 
sorptive and  dispersive  types  were  included. 

Because  the  RF  spectrum  analyzer  may  employ  a short  focal  length 
path  for  the  Fourier  transform  with  a beam  width  substantially  less  than 
a convenient  waveguide  coupled  photodetector  array  period,  a channel  fan- 
out waveguide  coupler  was  also  included  in  this  effort. 

The  following  Technical  Discussion  is  organized  by  subdivisions, 
reporting  on  the  computational  synthesis  results  for  various  Luneburg  len- 
ses, radial  refractive  index  profiles  and  focal  lengths,  together  with 
the  measured  lens  diffraction  patterns  and  data  for  the  sundry  components. 
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As  a result  of  this  effort,  several  scientific  papers  have  been  pre- 
pared for  publication  in  the  open  literature.  These  manuscripts  are 
included  as  appendices  as  follows: 

A.  D.  B.  Anderson  and  R.  R.  August,  "Progress  in  Waveguide 
Lenses  for  Integrated  Optics,"  Preprint  for  the  Trans- 
actions of  the  Institute  of  Electronics  and  Communications 
Engineers  of  Japan,  March  1978:  Special  Issue  Integrated 
Optics  and  Optical  Fiber  Communication. 

B.  S.  K.  Yao,  D.  B.  Anderson,  C.  M.  Oania,  and  V.  G.  Kreismanis, 
"Mask  Synthesis  for  Diffraction-Limited  Waveguide  Luneburg 
Lenses,"  Reprint  presented  at  the  Topical  Meeting  on  Inte- 
grated and  Guided  Wave  Optics,  January  16,  1978,  Salt  Lake 
City,  Utah. 

C.  J.  T.  Boyd  and  D.  B.  Anderson,  "Effect  of  Optical  Waveguide 
Scattering  on  the  Integrated  Optical  Spectrum  Analyzer  Dynamic 
Range,"  To  be  published  IEEE  Journal  of  Quantum  Electronics, 
June,  1978. 

D.  D.  B.  Anderson  and  R.  R.  August,  "Fourier  Transform  Signal 
Processing  in  an  Integrated  Optical  Format,"  Presented  at 
the  IEEE  IOSA  Conference  on  Lasers  and  Electro-Optical 
Systems,  February  8,  1978,  San  Diego,  California. 

An  effort  has  been  made  to  avoid  duplication  between  the  Technical  Dis- 
cussion and  the  various  appendices.  Appendix  A reports  the  status  of 
thin-film  waveguide  Luneburg  lenses  and  geodesic  lenses  as  of  mid-1977. 
Appendix  B delineates  progress  through  only  the  first  step  of  the  syn- 
thesis procedure  for  diffraction-limited  lenses.  The  exposition  con- 
tained in  Appendix  C is  a full  treatment  of  optical  waveguide  scattering 
and  its  effect  upon  the  integrated  optical  spectrum  analyzer  dynamic 
range  while  Appendix  D delineates  Rockwell's  "first-cut"  approach  to  an 
integrated  optical  spectrum  analyzer  utilizing  waveguide  Luneburg  lenses 
as  developed  under  this  effort.  Appendix  E,  "Integrated  Optical  Spectrum 
Analyzer  Signal  Level  Requirements,"  by  Professor  J.  T.  Boyd  addresses  the 
spectrum  analyzer  photodetector  and  CCD  image  array  dynamic  range  to  deter- 
mine the  Fourier  transform's  laser  power  requirements. 


II.  TECHNICAL  DISCUSSION 


A.  MULTILAYER  OPTICAL  WAVEGUIDES 

To  provide  a foundation  for  discussion  of  waveguide  Luneburg  lenses',7 
we  will  briefly  summarize  some  of  the  more  important  properties  of  bound 
modes  in  dielectric  waveguide.'8  In  particular,  we  will  consider  the  four- 
layer  structure  illustrated  in  Figure  2.  It  consists  of  a substrate 
having  a refractive  index  of  n4  and  superstrate  having  a refractive  index 
of  np  each  with  an  infinite  extent  with  one  or  two  thin  films  between 
which  serve  as  the  optical  waveguides  or  as  a waveguide  with  isolation 
layer.  The  planar  structure  considered  also  has  infinite  extent  in  y with 


Figure  2.  Cross  Section  of  Asymmetrical  Multilayer  Planar  Dielectric 
Waveguide  Showing  Zero  - Order  Mode  Bound  in  Both  Layers  with 
Evanescent  Field 

optical  propagation  in  the  z-direction.  For  wave  binding,  it  is  necessary 
that  ^ and/or  n3)  be  greater  than  (n^  and  n^).  We  will  not  consider 
the  modes  bound  to  the  substrate  when  it  is  of  limited  extent.  The  modes 
bound  to  this  type  of  structure  may  be  either  transverse  electric  (TEm)  or 
transverse  magnetic  (TMm)  of  order  m.  A transverse  field  component  bounded 
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to  both  layers  (overlayer  and  underlayer)  is  illustrated  in  Figure  2,  showing 
the  zero  -order  with  the  evanescent  field  extending  into  the  superstrate 
and  substrate.  It  is  possible  that  the  refractive  index  and  the  film  thick- 
ness of  these  layers  be  such  that  wave  binding  occurs  only  in  that  region 
which  has  the  highest  index  with  an  evanescent  component  in  the  other  layers. 
1 . Dielectric  Waveguide  Modes 

The  transverse  electric  modes  (TE)  are  the  eigen  solutions  of  the 
field  equation  where  the  time  variation  eia)t  has  been  suppressed 
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subject  to  the  following  boundary  conditions:  Ey  -*•  0 as  x approaches  + ®, 
and  the  tangential  components  of  E and  H are  continuous  across  the  boundaries. 
The  z dependence  will  be  the  same  in  all  four  regions  for  a given  mode  and 
can  be  written  as 
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where  m is  the  mode  order.  Additionally  e(x)  can  be  expressed  as  follows 
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Applying  the  boundary  conditions  on  E to  the  solutions  of  (1)  the  fields 
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where  Am  is  a normalization  constant,  and  Pm,  hm,  1_  and  qm  are  the  trans- 
it! m m m m 

verse  components  of  the  propagation  vector  in  regions  1,  2,  3,  4 respec- 
tively. Pm,  h , !l  , qm  and  are  related  by  the  following  dispersion 
equations 
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ilication  of  the  boundary  conditions  on  H yields 
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These  equations  when  combined  lead  to  the  following  transcendental  equation 
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The  roots  of  this  equation  are  the  allowed  values  of  the  propagation 
constants. 

The  transverse  magnetic  modes  (TM)  are  the  eigen  solutions  of  the 
field  equation 
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The  development  of  the  TM  modes  is  identical  to  that  for  the  TE  modes  and 
again  the  allowed  values  of  the  propagation  constants  are  roots  of  a trans- 
cendental equation,  which  is 
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2.  Thin-Film  Waveguide  Dispersion  Data 

Tantalum  pentoxide  and  Corning  7059  glass  have  been  employed  to 
form  waveguides  and  lenses  on  thermally-grown  Si02  substrates  because 
they  exhibit  comparatively  low  scattering  loss  and  may  be  deposited  with 
reproducible  properties  and  thickness  profiles.  Therefore,  the  dis- 
persion for  each  has  been  calculated  and  the  dispersion,  when  combined 
to  form  refractive  lens  elements,  has  also  been  calculated.  Dispersion 
data  is  also  presented  for  the  multilayer  structure  when  the  field  is 
trapped  in  only  the  dense  layer  and  when  it  is  trapped  in  both  layers. 


The  single  layer  thin-film  waveguide  dispersion  of  Ta205  struc- 
ture on  Corning  7440  glass  is  shown  in  Figure  3 for  both  TE  and  TM  modes 
as  a function  of  the  normalized  film  thickness.  The  thin-film  waveguide 
dispersion  of  a Corning  7059  glass  structure  on  an  Si02  substrate  is 
shown  in  Figure  4 for  both  TE  and  TM  modes  as  a function  of  the  normal- 
ized thickness.  The  effective  refractive  index  approaches  the  bulk 
value  for  thick  films  and  approaches  the  substrate  for  thin  films.  Each 
of  these  thin-film  waveguide  modes  exhibit  a cutoff  where  ne(m)  = 1.47 
for  Figure  3 and  ng(m)  = 1.47  for  Figure  4 in  a dielectric  waveguide 
structure.  Cutoff  means  the  lack  of  wave  binding  to  the  film.  An 
asymmetric  dielectric  structure  at  cutoff  leads  to  refraction  into  the 
substrate. 

The  dispersion  for  a two-layer  structure  for  only  the  transverse 
electric  modes  is  shown  in  Figures  5(a)  and  5(b)  where  Ta20g  is  em- 
ployed as  the  overlayer  of  variable  thickness  T and  where  the  underlayer 
normalized  thickness  is  constant  (kt  = 2.37  ir)  which  is  equivalent  to  a 
7059  layer  having  a thickness  of  0.75  microns  and  an  optical  wavelength 
of  0.63  microns.  Each  of  the  dispersion  curves  of  Figures  5(a)  and  5(b) 
exhibit  an  inflection  representing  the  transition  region  where  wave  bind- 
ing occurs  to  only  the  dense  overlayer  and  where  wave  binding  occurs  in 
both  layers.  Figure  5(b)  is  an  expansion  of  Figure  5(a)  in  the  region 
where  wave  binding  occurs  in  both  layers  showing  the  effect  of  substrate 
refractive  index. 

The  dispersion  for  a two-layer  structure  for  only  the  transverse 
electric  modes  where  7059  is  employed  as  the  overlayer  of  fixed  thick- 
ness and  where  the  underlayer  of  Ta205  of  variable  thickness  T is  shown 
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Figure  3.  Effective  Refractive  Index  Vs.  Normalized  Film  Thickness 

for  Ta^Og  Film  on  Si02  Substrate  TE  Modes  TM  Modes  
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Figure  4.  Effective  Refractive  Index  Vs.  Normalized  Film  Thickness  for 
7059  Glass  Film  on  SiO,  Substrate  TE  Modes  TM  Modes  - 


Figure  5.  Effective  Refractive  Index  of  Composite  Structure  as  a Function  of  Normalized  Thickness 
of  the  Dense  Overlay  Film  (Lens)  Where  the  Overlay  Film  Normalized  Thickness  kt  = 2.37tr 
for  TE  Modes  Only  (a).  Expanded  Scale  Showing  Effect  of  Three  Substrates  for  TE  Modes  Oi 


Figure  6.  Effective  Refractive  Index  of  Composite  Structure  as  a Function  of  Normalized  Thickness 
of  the  Dense  Underlay  Film  (Lens)  Where  the  Underlay  Film  Normalized  Thickness  kt  = 2.37ir 
for  TE  Modes  Only  (a).  Expanded  Scale  Showing  Effect  of  Three  Substrates  for  TE  Modes  Only  (b). 


in  Figures  6(a)  and  6(b).  The  fixed  thickness  waveguide  (kt  = 2.37  w)  is 
equivalent  to  a 7059  layer  having  a thickness  of  0.75  microns  for  the 
optical  wavelength  of  0.63  microns.  Each  of  the  dispersion  curves  of 
Figures  6(a)  and  6(b)  again  exhibit  an  Inflection  representing  the  trans- 
ition region  where  wave  binding  occurs  to  only  the  dense  underlayer  and 
wave  binding  occurs  In  both  layers.  Figure  6(b)  is  an  expansion  of 
Figure  6(a)  in  the  region  where  wave  binding  occurs  in  both  layers  show- 
ing the  effect  of  substrate  refractive  index  for  the  same  three  values  of 
refractive  index. 

This  dispersion  data  is  sufficient  for  design  of  homogeneous  and  inhomo- 
geneous refractors  and  reflector  lens  systems.  In  general,  only  transverse 
electric  modes  are  employed  in  the  lowest  order  to  minimize  loss.  Res- 
triction to  the  TE  mode  is  achieved  by  input  coupling  polarization.  The 
thickness  of  some  lenses  will  permit  the  presence  of  a multiplicity  of 
mode  orders,  however,  most  transitions  are  sufficiently  gradual  to  minimize 
excitation  of  higher  order  modes.  Location  of  the  highest  index  layer  ad- 
jacent to  the  substrate  maximizes  conversion  to  the  next  higher  order  mode 
which  is  odd.  It  is  also  possible  to  prohibit  mode  conversion  by  symmetry, 
that  is,  by  interleaving  the  denser  film  between  two  separate  halves  of 
the  waveguide. 

The  following  sections  consider  thin-film  waveguide  generalized 
Luneburg  lenses  having  circular  symmetry  requiring  that  the  radial  refrac- 
tive index  profile  follow  a prescribed  function  to  obtain  a focus  of  the 
radiation.  The  following  numerical  procedure  will  utilize  the  above  analysis 
in  the  computations  to  derive  the  overlay  waveguide-lens  thickness  such  that 
it  will  satisfy  the  radial  refractive  index  dependence.  It  will  also  include 
the  effect  of  errors  arising  from  the  waveguide  thickness  and  the  lens 
radial  functional  thickness  dependence. 
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B. 


INHOMOGENEOUS  REFRACTIVE  WAVEGUIDE  LENSES 


Conventional  optical  lens  systems  employ  homogeneous  dielectrics 
with  spherical  surfaces  in  various  combinations  to  correct  aberrations. 
Inhomogeneous  materials  with  weak  perturbations  randomly  distributed  are 
avoided.  Generally,  materials  with  large  inhomogeneous  perturbations 
are  not  available  and,  thus,  lenses  requiring  large  refractive  index  vari- 
ations such  as  Maxwell's  fish-eye’9and  the  Luneburg  lens20have  not  been 
realized  in  the  optical  portion  of  the  spectrum.  Thin-film  optical  wave- 
guide provides  the  opportunity  to  synthesize  thin-film  inhomogeneous  di- 
electric lens  using  homogeneous  materials  by  variation  of  the  film  thick- 
ness. 

1 . "Classical"  Luneburg  Lens20 

Luneburg  lens  is  an  inhomogeneous  positive  refractor  with  radial 
symmetry  which  images  a plane  wavefront  to  a hemispherical  spot  located 
on  the  opposite  boundary  of  the  refractor.  Ray  traces  through  the  "clas- 
sical" Luneburg  lens  are  illustrated  in  Figure  7.  Interest  in  the  Luneburg 
lens  arises  because  of  its  4tt  field  of  view  and  because  it  is  free  of  all 
aberrations  except  field  curvature.  The  cross  section  of  a dielectric 
waveguide  Luneburg  lens  is  illustrated  in  Figure  8 where  a circular  over- 
lay film  is  employed  with  prescribed  radial  thickness  function.  For 
integrated  optics,  its  simplicity  of  fabrication  in  a thin-film  structure 
realized  reproducibly  is  an  important  attribute.17  The  circular  symmetry  is 
of  particular  interest  where  the  substrate  area  is  limited  because  fold- 
ing of  the  optical  axis  may  use  the  same  lens  for  several  purposes.  The 
waveguide  Luneburg  lens  has  a further  advantage  that  all  refractive  index 
variations  are  continuous,  smooth  and  slow  functions  with  respect  to  the 
optical  wavelength.  Mode  conversion  in  the  lens  is  suppressed  by  the 
graded  thickness  profile  and  by  symmetry. 

The  Luneburg  lenses  may  be  generalized  to  form  a perfect  geometric 
image  of  two  given  concentric  circles  on  each  other.2'  The  perfect  focusing 
property  may  be  achieved  in  various  regions.  If  one  of  the  circles  is  of 
infinite  radius,  the  lens  may  be  designed  to  focus  on  another  circle  outside 
of,  at  the  boundary,  or  inside  of,  the  Luneburg  lens.  Luneburg  found  a 
solution  for  the  refractive  index  function  which  had  a conjugate  foci  for 
two  given  points  outside  of  the  lens.  The  solution  for  a parallel  beam 
array  which  focuses  on  the  opposite  lens  interface  requires  that  the 
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Figure  7.  Ray  Trajectory  Through  Luneburg  Lens  Where  f-Number  = 0.5, 
which  is  Equivalent  to  Morgan's  S = 1. 


Figure  8.  Cross  Section  of  a Waveguide  Luneburg  Lens  using  a Dense 

Ta^O^  Overlay  Lens  on  Corning  7059  Waveguide  on  a Thermally- 
Oxidized  Silicon  Substrate. 


refractive  index  n(r,m) 

n(r,m)  = ne(m)  2 - (r/rQ)2  r < rQ  (12) 

where  ng  is  the  refractive  index  of  the  medium  surrounding  the  lens,  rQ 
is  the  lens  radius,  and  r is  the  radial  coordinate  for  mode  (m).  Southwell 
has  developed  a technique  for  ray  tracing  through  Luneburg  lenses  which  has 
been  employed  to  generate  Figure  7 (The  corresponding  diffraction  pattern 
is  depicted  in  Figure  14.). 


Generalized  Luneburg  Lenses 

The  "classic"  Luneburg  lens  has  been  generalized  by  Morgan2,devel op- 
ing a radial  refractive  index  profile  which  will  focus  one  circle  on 
another  in  the  region  beyond  the  edge  of  the  lens  radius  r . The  increased 
focal  length  reduces  the  range  of  refractive  index  variation  required  to 
form  the  lens.  Extending  the  focus  outside  of  the  lens  region  increases 
the  f-number  and  thus  the  diffraction  pattern  width  in  the  focal  region. 


Figure  9.  Ray  Trace  Geometry  Through  Luneburg  Lens  and  Ray  Trace 
Intercept  Plot. 

The  geometry  of  a generalized  Luneburg  lens  is  depicted  in  Figure  9 
showing  a perfect  ray  trace.  An  imperfect  ray  trace  is  also  included  in 
Figure  9 from  which  a ray  intercept  plot  is  developed  which  is  also  in- 
cluded in  Figure  9.  In  the  following  section,  numerous  ray  traces  are 
calculated  yielding  the  ray  intercept  plot  from  which  the  phase  wavefront 
has  been  derived  followed  by  the  calculated  intensity  diffraction  pattern. 
This  procedure  has  provided  the  means  to  calculate  generalized  Luneburg  lens 
performance  starting  from  waveguide  thickness  profiles  with  various  func- 
tional errors  to  determine  the  allowed  tolerances  as  determined  by  degrada- 
tion of  the  diffraction  pattern. 
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The  functional  radial  refractive  index  profile  for  some  generalized 
Luneburg  lenses  can  be  obtained  from  a paraxial  approximation  in  the  region 
where  the  f-number  is  >1.5.  The  normalized  refractive  index  profile  is 
given  by 


^ne(m)2r02  _ f2n2(r1m) 


’^ne(m)r1 


(13) 


In  the  above  equation,  r is  the  radial  coordinate  from  the  lens  center 
along  the  waveguide  surface,  rQ  is  the  radius  of  the  lens  region,  the  focal 
length  is  r^=F  and  ng(m)  is  the  effective  refractive  index  for  the  wave- 
guide mode  (m)  of  interest.  Note  that  the  above  equation  is  not  solved  for 
n(r,m),  but  n(r,m)  is  readily  determined  from  Eq.  (13)  numerically  as  a 
function  of  r-j.  Using  Eq.  (13),  Figure  10  depicts  the  normalized  effective 
refractive  index  profile  as  a function  of  the  normalized  lens  radius  for 
several  values  of  f-numbers.  It  indicates  that  the  transverse  profile  be- 
comes more  uniform  having  a lower  refractive  index  over  a larger  central 
region  of  the  lens  as  the  focal  length  is  increased. 

Note  that  Equation  (13)  depends  upon  the  waveguide  mode  order  and, 
therefore,  a waveguide  which  will  support  more  than  a single  mode  will 
have  a separate  foci  for  each  of  the  modes  of  which  only  one  can  yield  a 
diffraction-limited  pattern. 

For  the  classic  Luneburg  lens  n(0,0)/ne(0)=  %/2  , a high  index  layer 
like  Ta2°5  or  Nb2C)5  1S  reclu''recl  on  Corning  7059  for  F/D  = 0.5.  Generalized 
Luneburg  lenses  with  longer  focal  lengths  may  be  formed  with  these  materials 
using  thinner  films  or  using  films  with  lower  index  and  thicker  layers. 
Generalized  Luneburg  lenses  have  been  realized  using  an  "all"  7059  thin 
film  when  F/D>6.  Lenses  which  are  formed  from  a single  composition  by 
changes  of  the  layer  thickness  exhibit  substantially  less  scattering  and 
further  suppress  the  excitation  of  higher  order  modes. 
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Generalized  Luneburg  Lens  Analysis 

22  23 

Southwell's’  analytical  approximation  is  derived  for  Luneburg's 
integral  equation  from  which  the  index  profile  for  the  generalized 
Luneburg  lens  may  be  extracted  to  sufficient  accuracy  for  ray-tracing 
to  diffraction-limited  performance.  This  expression  was  derived  as  an 
expansion  based  on  a consideration  of  known  behavior  at  certain  limits. 

The  expression  consists  of  five  parameters  which  depend  on  the  generalized 
Luneburg  lens  focal  length.  Two  of  these  parameters  are  determined 
from  the  "boundary"  limits  and  the  others  are  determined  using  a least- 
squares  procedure.  This  analytical  expression  provides  a convenient  means 
for  obtaining  the  radial  index  n(r)  for  ray-tracing  purposes.  It  also 
provides  a means  of  conveniently  determining  integrated  optical  Luneburg 
lens  thickness  profiles. 

A spherically  symmetric  index  distribution  n(r)  which  gives  a 
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perfect  image  of  an  infinite  object  satisfies  the  integral  equation 


n = expMp,s)L 
where  p=rn  and 


u(p 


•*>-i  / 


arcsin(x/s)  . 
(x2.p2)l/2  • 


(14) 


(15) 


in  which  the  optical  guided  wave  index  of  refraction  variable  n is  normalized 
to  the  optical  guided  wave  index  of  refraction  ng  of  the  homogeneous  wave- 
guide medium  surrounding  the  lens  such  that  n>l . Both  the  radius  variable 
r and  the  image  distance  parameter  s as  measured  from  the  center  of  the 
lens  are  normalized  with  respect  to  the  generalized  Luneburg  lens  radius 
rQ  such  that  0<r<l  and  s>l.  The  full  aperture  f-number  for  these  lenses 
turns  out  to  be  approximately  s/2  for  s>l . 


For  the  case  s=l , the  above  integral  may  be  evaluated  in  closed  form 
yielding  the  "classical"  Luneburg  lens.  To  obtain  a closed  form  analyti- 
cal approximation  for  the  generalized  Luneburg  lens  given  by  Equation  (15), 
a change  of  variable  is  required.  Letting  x'  = x/s  and  p=0,  the  value  of 
the  generalized  Luneburg  lens  refractive  index  maxima  which  occurs  at  the 
center  is  given  by  Equation  (14)  and 
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TABLE  I.  THE  FUNCTION  w(o,s)  USING  EQ.  (16)  n = expMo.s)] 

MlaX 


s 

w(o,s) 

"max(rs0> 

1 

0.34657359 

1.41421356 

2 

0.161532974 

1.17531121 

3 

0.1067/8890 

1.11268820 

4 

0.0798585731 

1.08313387 

5 

0.0638050027 

1.06588453 

6 

0.0531341389 

1.05457109 

7 

0.0455246840 

1.04657684 

8 

0.0398234213 

1.04062701 

9 

0.0353921041 

1.03602586 

10 

0.0318487205 

1.03236132 

17 

0.0187277137 

1.01890418 

a)(o,s) 


1/s 

arcsinx 


I 


dx 


(16) 


The  refractive  index  maxima  for  Luneburg  lenses  are  given  in  Table  I as 
well  as  by  the  ordinate  of  Figure  10. 

To  investigate  the  Luneburg  refractive  index  profile  for  p^o,  a 
new  variable  y=x-p  is  employed  for  Equation  (15)  which  becomes 


co(p,s) 


(17) 


We  may  now  remove  the  y"^2  singularity  by  performing  an  integration  by 
parts , 
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-1/2 


w(p,s) 


1 2 arcs  in  (j) 

* W (Hp)177 

arcsin  *±£\  l/2d 


(18) 


From  Equation  (18),  it  is  observed  that  as  p+l , which  is  near  the  edge 
of  the  lens  and  in  the  limit  as  y-+0,  the  integrand  of  Equation  (18)  be- 
haves as  y1 72 . Thus,  it  is  found  that  the  integral  portion  behaves  as 
(1-p) 3/1,2 , which  is  a higher  order  term  compared  to  the  first  term  on  the 
right  side  of  Equation  (18).  Thus,  in  the  limit  as  pr*T  the  integral  may 
be  ignored  and  oj(p,s)  behaves  like 


w(p,s) 

P+1 


(19) 


Furthermore,  the  integral  of  Equation  (18)  indicated  that  the  next  order 
term  will  have  a (1-p)372  behavior.  Therefore,  it  is  reasonable  to  sug- 
gest a power  series  for  the  function  w(p,s)  as  follows 


a)(p,s)  = p-,  (1  -p)1/2  + p2(l-p)3/2  + p3(l-p)5/2  + p4(l-p)7/d  + p5(l-p)9/2  (20) 

Obviously, 

pi  = arcsin  m (21) 

and 


P5  = ai(o,s)  - p1  - p2  - p3  - p4  . (22) 

Thus,  only  three  of  the  p coefficients  are  as  yet  unspecified.  To  deter- 
mine these  remaining  coefficients,  we  fit  the  numerical  values  of  data 
given  in  Equation  (18)  to  the  function  oa(p,s)  given  by  Equation  (20). 

More  points  are  included  near  the  p=l  region  in  order  to  more  heavily  weight 
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TABLE  II 


RESULTS  OF  FITTING  EQ.  (20)  TO  THE  DATA  IN  TABLE  I. 

The  values  listed  for  SUM  are  the  sum  of  the  squares  of  the 
errors  for  the  29  points.  p5  was  determined  by  Eq.  (22). 


s=2 

s=3 

s=5 

$=9 

SUM 

3.9101x1 0_ 1 ^ 

5.81 34x1 0_ 1 ^ 

1.7893xl0"12 

5.3512xl0"13 

pi 

0.235687835 

0.152976535 

.0906399959 

.0501194645 

P2 

-7.47500358x10"^ 

-4. 22494061 xlO"2 

-2. 34300068x1 0"2 

-1.26356394xl0"2 

P3 

6. 72894476xl0~3 

-1 .93175172x1 0-3 

-2.511 0001 7x1 0-3 

-1 .6393441 8x1 0"3 

P4 

-5.1 4447054x1 0"3 

-8.25244897xl0"4 

-1.49192458xl0"4 

-4.07947091xl0"5 

-9.89299661x10"^ 

-1 .191 24255x1 0~  3 

-7 . 4479371 7x1 0"4 

-4.11 5821 22x1 0-4 

the  region  when  the  greatest  change  is  occurring.  Table  II  is  a tabulation 
of  the  results  of  the  curve-fit  for  the  29  "data"  points  determined  from 
the  numerical  integration  of  Equation  (18).  The  deviation  of  this  approxi- 
mation from  the  exact  formula  has  been  examined  and  is  found  to  be  less 
than  10"®. 

Having  an  analytical  expression  for  co(p,s),  we  are  now  in  a posi- 
tion to  extract  n(r)  numerically  for  any  specified  value  of  r.  To  do 
this,  we  define  a function  F, 

F(n)  = exp[w(p,s)]  - n.  (23) 

For  a given  r and  s,  we  extract  the  value  of  n that  makes  F(n)  vanish  using 
Newton's  iterative  method.  A tabulation  of  results  is  given  in  Table  III. 
Plots  of  these  profiles  are  shown  in  Figure  10. 

In  addition  to  a numerical  value  of  n,  the  ray-trace  algorithms” 

2 2 

also  require  values  for  dn/dr  and  d n/dr  . The  analytical  expressions 

2 2 

for  dn/dr  and  d n/dr  in  terms  of  r and  the  numerical  value  for  n can  be 

derived  using  implicit  differentiation  of  Equation  (23)  using  Equation  (20) 

for  w(p,s).  The  results  are 
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TABLE  III 

NORMALIZED  REFRACTIVE  INDEX  PROFILES  FOR  GENERALIZED  LUNEBURG  LENSES 


r 

s=2 

s=3 

s=5 

s=9 

0 

1.175311212 

1.112688200 

/ 1.065884534 

1.036025859 

0.05 

1.174999407 

1.112507930 

1.065788087 

1.035976411 

0.10 

1.174071705 

1.111969546 

1.065500332 

1.035828873 

0.15 

1.172521576 

1.111069129 

1.065018948 

1.035581965 

0.20 

1.170340993 

1.109801290 

1.064340603 

1.035233817 

0.25 

1.167519780 

1.108158811 

1.063460759 

1.034781859 

0.30 

1.164044879 

1.106132212 

1.062373426 

1.034222692 

0.35 

1.159899514 

1.103709214 

1.061070836 

1.033551908 

0.40 

1.155062156 

1.100874045 

1.059543028 

1.032763858 

0.45 

1.149505212 

1.097606509 

1.057777265 

1.031851324 

0.50 

1.143193243 

1.093880684 

1.055757219 

1.030805062 

0.55 

1 . 1 36080437 

1.089663045 

1.053461790 

1.029613132 

0.60 

1.128106836 

1.084909633 

1.050863316 

1.028259883 

0.65 

1.119192434 

1.079561599 

1.047924746 

1.026724355 

0.70 

1.109227400 

1.073537799 

1.044594938 

1.024977603 

0.75 

1.098054844 

1.066721643 

1.040800259 

1.022977928 

0.80 

1.085437716 

1.058935566 

1.036428193 

1.020661529 

0.85 

1.070987513 

1.049884988 

1.031290974 

1.017921692 

0.90 

1.053981685 

1.039009952 

1.025027636 

1.014552365 

0.91 

1.050159236 

1.036524794 

1.023580201 

1.013768729 

0.92 

1.046155078 

1.033902395 

1.022045349 

1.012935510 

0.93 

1.041942017 

1.031119755 

1.020407545 

1.012043672 

0.94 

1.037484456 

1.028146025 

1.018645667 

1.011080865 

0.95 

1.032734249 

1.024938133 

1.016729693 

1.010029408 

0.96 

1.027623448 

1.021432662 

1.014614324 

1.008862349 

0.97 

1.022050589 

1.017529080 

1.012225136 

1.007534739 

0.98 

1.015851602 

1.013049431 

1.009422596 

1.005960167 

0.99 

1.008726930 

1.007614146 

1.005874179 

1.003921693 

1.00 

1 

1 

1 

1 
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(24) 


dn  = -n  G 
W 2(l-rn)1/2+nrG 


d2n  _ 2o3[g2  (l-rn)1/2+H  . 2 

— 5"  “ r T77 + — 

dr' 


f dn  i 


where 


2(l"-rn)1/2+nrG  3 n W 


G = p1  + 3p2(l-rn)  + 5p3(l-rn)2  + 7p4(l-rn)3  + 9p5(l-rn)4 


(25) 


(26) 


and 


H =-p1  + 3p2(l-rn)  + 15p3(l-rn)2  + 35p4(l-rn)3  + 63p5(l-rn)4  . (27) 


The  design  of  the  waveguide  Luneburg  lens  requires  a specification 
of  the  kind  of  waveguide  mode  and  its  order.  For  a dense  overlay  lens  on 
a waveguide  on  an  isolation  substrate,  the  equations  governing  the  rela- 
tionship between  the  waveguide  effective  index  of  refraction  n and  the 
overlay  lens  thickness  t are  given  by  Equations  (9)  and  (11)  for  the  TE 
modes  and  TM  modes,  respectively.  When  rewritten  as  explicit  form  for 
the  lens  overlay  thickness  required  to  produce  certain  waveguide  effective 
index.  Equation  (9)  becomes 


tan 


-1 


h [l  - I tan(£d)]  + |^-  | + tan(£d) 


I [q  + tanUd)]  - P [l  - ^ tan(£d)  J 


n<n. 


ht 


(28) 


tan 


-1 


h 

2£d 

(q  + £)e  + q - £ 

+ EtL 

£ 

2£d  -j 

(q  + £)e  - q + £ J 

h2 
. £ 

2£d 

(q  + £)e  - q + £. 

- P 

2£d  I 

(q  + £)e  + q - £ J 

n>n 


3' 


for  the  TE  modes.  If  all  the  parameters,  p,  h,  q,  are  replaced  by 
p h i q Equation  (28)  can  be  employed  for  the  TM  modes. 

7T’  7T*  “T*  77* 


"l  - n2  n3  n4 
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Figure  11.  Waveguide  Overlay  Thickness  Profile  for  Generalized  Luneburg 
Lenses  Using  Ta?0,-  on  Corning  7059  (d=1.0665y)  on  Si0o  Which 
is  Tabulated  in^Table  IV.  * 


Consequently,  a Luneburg  lerrs  waveguide  thickness  profile  can  be 
obtained  from  the  above  numerical  determination  of  n(r).  Some  numerical 
examples  have  been  calculated  for  a Ta205  lens  overlay  having  a refractive 
index  of  2.1  using  a Corning  7059  waveguide  having  a refractive  index  of 
1.565  and  a thickness  of  1.0665y  on  a Si02  substrate  having  a refractive 
index  of  1.47  and  using  a wavefront  of  0.9y.  Figure  11  illustrates  the 
waveguide  lens  thickness  profile  with  the  corresponding  numerical  data 
tabulated  in  Table  IV  for  s=2,  3,  5,  9. 

To  test  the  accuracy  of  the  approximation  of  Eq.  20,  we  have  ray 
traced  the  s=2  case  using  Southwell's  procedure  which  provides  the  wave- 
front  phase  errors  and  the  intensity  diffraction  patterns.  Figure  12, 
Section  (a),  is  a plot  showing  six  rays  traced  through  an  s=2  Luneburg 
lens.  Nineteen  rays  were  traced  and  the  resulting  wavefront  phase  error 
at  X = 0.9y  was  determined.  This  wavefront  phase  error  result  is  shown  in 
Figure  12,  Section(b),  with  a corresponding  intensity  diffraction  pattern 
shown  on  a logarithmic  scale  in  Section (c)of  this  figure.  Note  that  the 
refractive  index  profile  as  derived  from  Eq.  20  produces  diffraction- 
limited  results. 
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OVERLAY  THICKNESS  IN  MICRONS  FOR  OPTICAL  WAVEGUIDE  LUNEBURG  LENSES 

The  waveguide  parameters  are  n,=l,  n?*2.1,  n~=1.565,  n.s1.47,  d=1.0665y, 
and  X=0.9y  1 


r 


(cm) 

s=2 

s=3 

s=5 

s=9 

0 

0.2191 

0.1591 

0.1223 

0.0979 

0.05 

0.2187 

0.1590 

0.1222 

0.0979 

0.10 

0.2177 

0.1585 

0.1220 

0.0978 

0.15 

0.2160 

0.1578 

0.1216 

0.0975 

0.20 

0.2135 

0.1567 

0.1211 

0.0972 

0.25 

0.2105 

0.1554 

0.1204 

0.0968 

0.30 

0.2068 

0.1537 

0.1196 

0.0963 

0.35 

0.2024 

0.1518 

0.1186 

0.0957 

0.40 

0.1975 

0.1495 

0.1174 

0.0949 

0.45 

0.1920 

0.1469 

0.1160 

0.0940 

0.50 

0.1860 

0.1440 

0.1141 

0.0930 

0.55 

0.1794 

0.1406 

0.1126 

0.0918 

0.60 

0.1722 

0.1370 

0.1105 

0.0905 

0.65 

0.1646 

0.1328 

0.1081 

0.0888 

0.70 

0.1563 

0.1282 

0.1054 

0.0869 

0.75 

0.1472 

0.1229 

0.1022 

0.0846 

0.80 

0.1374 

0.1169 

0.0983 

0.0817 

0.85 

0.1262 

0.1097 

0.0935 

0.0780 

0.90 

0.1130 

0.1006 

0.0870 

0.0727 

0.91 

0.1099 

0.0984 

0.0853 

0.0713 

0.92 

0.1067 

0.0960 

0.0835 

0.0697 

0.93 

0.1031 

0.0933 

0.0814 

0.0679 

0.94 

0.0993 

0.0903 

0.0790 

0.0657 

0.95 

0.0949 

0.0869 

0.0762 

0.0632 

0.96 

0.0898 

0.0827 

0.0728 

0.0601 

0.97 

0.0835 

0.0774 

0.0682 

0.0559 

0.98 

0.0749 

0.0699 

0.0616 

0.0499 

0.99 

0.0597 

0.0562 

0.0496 

0.0395 

1.00 

0 

0 

0 

0 
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Figure  12.  Ray  Traces  of  an  s=2  Generalized  Luneburg  Lens  Using  the  Index 
Profile  Derived  from  Eq.  14  , the  Wavefront  Phase  Error  in 

the  Pupil  Plane  (b),  and  the  Corresponding  Intensity  Diffrac- 
tion Pattern  Using  a Logarithmic  Scale  in  the  Image  Plane  (c). 

4.  Inhomogeneous  Lens  Analysis 

Southwell's  procedure  for  gradient-index  ray-tracing  techniques  to 
evaluate  the  performance  on  inhomogeneous  optical  lenses  is  derived  herein 
for  reference  because  of  its  application  to  generalized  waveguide  Luneburg 
lenses.  Using  the  thin-film  waveguide  parameter  including  the  lens  thick- 
ness profile,  the  wavefront  phase  error  and  intensity  diffraction  patterns 
are  derived.  The  procedure  is  applied  to  a classical  Luneburg  lens  index 
profile  for  which  exact  results  are  known.  The  procedure  yields  better 
than  diffraction-limited  accuracy  within  reasonable  computer  running  time. 
Restricting  our  interest  to  lens  index  profiles  having  circular  symmetry, 
the  ray  path  position  vector  r shown  in  Fig.  13  can  be  described  in  a 
coordinate  system  with  its  origin  at  the  center  of  symmetry  of  the  index 
profile.  The  position  of  the  ray  is  described  by  the  components  of  r. 
Defining  the  ray  parameter  s as  a scalar  distance  along  the  ray  path  mea- 
sured from  an  arbitrary  reference  point  on  the  ray,  the  unit  ray  direction 
vector  S,  which  is  tangent  to  the  ray  curve,  is  therefore 
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Figure  13.  Ray  Path  Geometry  Showing  Ray  Position  Vector  r and  the  Ray 
Vector  S.  The  Scalar  s is  the  distance  Along  the  Path  and 
the  Vector  S is  a Unit  Vector  Tangent  to  the  Path. 
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= x i + y 
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j 

(29) 

dr  j 

= 3?  " ai 

A 

+ Pi 

(30) 

where  a and  3 are  the  directional  cosines  of  the  ray.  The  characteristics, 
namely  the  position  and  direction  vectors,  of  the  ray  at  any  point  (s+As) 
may  be  obtained  from  the  ray  characteristics  at  the  previous  point  (s)  by 
Taylor's  series  expansion.  That  is 


r(sQ  + As)  = 


S(sQ  + As)  = 
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(31) 


(32) 


dr 


Noting  that  is  the  ray  direction  vector,  the  second  derivative  of  r is 

-►  o-t-  2 -»■ 

the  curvature  vector  K = d r/ds  = dS/ds.  With  these  two  Taylor  expansions, 
the  ray  can  be  traced  point  by  point  to  any  distance,  provided  that  the  ray 
characteristics  of  an  initial  point  and  the  necessary  derivatives  are  given. 

In  order  to  obtain  the  necessary  derivatives,  use  can  be  made  of  the 
equation  of  motion  for  the  ray  position  in  a gradient  index  medium,  which  is 
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■ ... 


d_  I"  dr  I 
ds  [n  ds  J 


vn 


where  vn  is  the  index  gradient 


n„  _ 3n  ; . 9n  q 

vn  - 37  1 + 37  J 


(33) 


(34) 


Expanding  Eq.  (33),  we  obtain  the  curvature  vector 
it  = 1 [vn  - (vn  • 1)  1 ] , 

Before  obtaining  the  last  required  derivative  dK/ds,  we  rewrite  Eq.  (35), 
making  use  of  the  symmetry  in  n 

It  = f (r)  [r  - (r  • s)  S ] 

where 


f (r) 


_ 1 dn 
nr  dr  ’ 


(35) 


(36) 

(37) 


Differentiating  Eq.  (36),  we  have 

jf  - 1 % (I  • ?)[?  - (r  • DS  - f[(r  • t) J t (f  • s)lt]  (38) 

where  we  have  used  the  general  relationships, 

jjf  - S • vg,  and  Vg(r)  = £ ^ r (39) 

The  ray  tracing  algorithm,  therefore,  is  as  follows: 

(1)  From  the  known  position  and  directions  (x0,  y0,  a0,  B0)  and  index  gradient, 
evaluate  the  curvature  K and  its  derivative,  Eqs.  (35)  and  (38)  at  that 
point. 

(2)  Select  an  incremental  step  size  as  and  evaluate  the  new  position  and 
angle  from  Eqs.  (31)  and  (32). 

(3)  Considering  our  new  position  and  angle  as  known,  go  back  to  Step  (1)  and 
repeat  the  sequence  until  the  exit  surface  is  reached  or  until  the  region 
is  encountered  when  the  index  gradient  vanishes. 


♦* 
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We  note  that  this  algorithm  requires  the  index,  the  gradient  dn/dr, 

2 2 

and  the  second  derivative  d n/dr  at  each  point.  In  case  of  ray  tracing 

for  an  ideal  generalized  Luneburg  lens,  the  power  series  for  the  index  n, 

dn/dr,  and  d n/dr  are  given  by  Eqs.  (20),  (32)  and  (33),  using  the  five 

P parameters.  In  case  of  ray  tracing  from  either  a measured  or  a computed 

thin  film  Luneburg  lens  thickness  profile,  the  index  profile  is  obtained 

from  thickness  profile  as  numerical  data  points.  The  derivatives  dn/dr 
2 2 

and  d n/dr  can  be  readily  differentiated  numerically  from  the  index  profile 
data . 

It  is  necessary  to  consider  the  problem  of  precisely  locating  the 
ray  position  as  it  exits  the  lens  region  because  any  error  here  will 
result  at  a focal  plane  error;  i.e.,  aberration;  of  the  same  magnitude. 

For  computational  efficiency,  the  step  size  AS  is  large  compared  to  the 
required  tolerances.  Thus,  an  interpolation  is  required  once  the  ray  has 
stepped  outside  the  lens. 

The  procedure  we  take  for  locating  the  edge  is  as  follows.  We  first 
define  the  function  F 

F = r02  - r2  (40) 

2 2 2 

where  rQ  is  the  lens  radius  and  r = x + y as  determined  by  the  ray  position. 
By  considering  the  last  point  inside  the  lens  as  the  base  point,  we  see 
that  r becomes  a function  of  the  step  size  As.  The  function  F then  may  be 
considered  a function  of  as,  F = F(as).  The  problem  of  locating  the  edge 
is  then  solved  by  finding  the  root  as  that  makes  F(as)  = 0.  This  may  be 
done  using  Newton's  method 


(As)e  = (as)q 


Aso 


(41) 


Since  the  functional 


dependence  of  F on  as  is  given  through  Eq.  (31),  we  have 


aF 

3(as) 


-2r 


dr 

dlAiT 


-2r  • s . 


(42) 


L 
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Initially  (as)q  = 0 at  our  base  point  inside  the  lens  and  the  derivative 
in  Eq.  (41)  has  already  been  evaluated.  It  is  thus  a simple  matter  to 
apply  Eq.  (41)  until  |F|  < 10  or  some  prescribed  small  quantity. 

There  are  two  factors  to  be  considered  when  determining  a value  for 
the  step  size  as.  They  are  computing  time  and  accuracy.  When  tracing 
diffraction-limited  systems,  the  ray  trace  must  provide  an  exit  height 
accurate  to  Some  fraction  of  a wavelength.  This  may  require  extremely 
small  step  size  and,  consequently,  long  computing  time. 

With  a classic  Luneberg  lens  profile,  the  ray  tracing  error  due 
to  finite  incremental  step  size  as  can  be  evaluated  because  the  theoretical 
result  will  have  no  aberrations.  Southwell  has  shown  that  for  diffraction 
limited  ray  tracing,  the  ray  trace  error  must  be  smaller  than  a few 
wavelength  and  that  the  ray  trace  step  size  may  be  approximately  a few 
tens  of  wavelength.  To  reduce  the  computing  time,  a compensation  scheme 
based  on  extrapolation  of  ray  trace  error  due  to  finite  step  size  has  been 
developed  by  Southwell. 

The  spot  diagrams  determined  by  ray  tracing  may  be  employed  to 
determine  the  optical  phase  front  error  through  the  system.  These  phase 
errors  are  then  expressed  in  a pupil  function  (which  is  a function  of  the 
pupil  coordinates)  and  is  then  Fourier  transformed  and  squared  to  form  the 
image  diffraction  pattern. 

To  derive  the  wavefront  errors  from  ray  trace  intercept  data,  we 
draw  upon  the  theory  of  Hamilton's  angle  characteristic  function  T(e), 
which  represents  the  optical  path  length  along  the  path  through  the  system 
described  by  the  entrance  height  y0.  This  function  is  also  expressed  as  a 
function  of  the  exit  angle  e.  The  theory  states  that  the  ray  intercept 
error  y in  the  image  plane  is  given  by 

y • - f . (43) 

To  utilize  this  theory,  we  assume  that  the  optical  path  length  may  be  expressed 
as  a power  series  in  the  exit  angle  coordinate  &, 

T = c^2  + c2e4  + c3b6  + . . . (44) 
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Since  we  are  dealing  with  symmetric  systems  we  consider  only  even  powers  of 
3.  The  ray  trace  intercepts  y^  and  the  associated  angles  3-j  can  be 
obtained  by  ray  tracing.  After  a few  traces,  the  yi , Bi  are  fitted  to  a 
least  square  deviation  fashion  to  the  function 

f|  = (2c]B  + 4c263  + 6c3B5  +...).  (45) 

This  procedure  determines  the  c coefficients  for  the  optical  path  length  in 
Eq.  (44). 

Prior  to  fitting  the  intercept  data  to  the  expression  (45),  it  is 
convenient  to  first  normalize  the  b’s.  By  doing  this,  the  wavefront  error 
function  (44)  is  then  automatically  expressed  in  terms  of  normalized  aperture 
coordinates.  The  preceding  statement  is  valid  strictly  for  systems  with  a 
linear  relationship  between  aperture  and  exit  angle.  Although  this  is  only 
approximately  true,  it  is  sufficiently  close  for  our  case. 

The  phase  error  is  2tt/x  times  the  optical  path  error.  Figure  14(a) 
is  a plot  of  the  phase  error  of  the  Luneberg  lens  of  Figure  7 using  extrapola- 
tion on  the  rays  with  As  = 0.1  cm  for  the  initial  step  size.  Since  we 
expect  zero  phase  error,  the  amount  shown  represents  the  residual  ray 
tracinq  error.  We  note  that  it  is  considerably  less  than  tt/4  and  consequently 
represents  a diffraction-limited  wavefront. 

The  pupil  function  is  given  bv  the  expression 

P = exp[i27rT(y/y0)/x],  (46) 

where  the  function  T is  expressed  in  terms  of  normalized  aperture  position 
y/Yo  • The  square  of  the  Fourier  transform  of  P yields  the  diffraction  pattern 
in  the  image  plane.  This  is  performed  using  the  discrete  fast  Fourier 
transform  algorithm.  Figure  14(a)  is  the  diffraction  pattern  for  the  "classic" 
Luneburg  lens  as  traced  in  Figures  7 and  14(a). 

To  further  test  the  system,  we  moved  the  detector  plane  slightly 
away  from  the  edge  of  the  lens  and  repeated  the  evaluation.  The  corresponding 
wavefront  phase  error  and  its  diffraction  pattern  are  depicted  in  Figure  14(b). 


34 


Phase  Error  in  Pupil  Plene 
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Figure  14.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  "Classic”  Geodesic  Lens  Ray  Traced  in  Figure  7 
for  Perfect  Radial  Profile  (a)  and  Where  the  Image  Plane  Has 
Been  Displaced  by  X Yielding  the  Defocus  Aberration  (b). 


5.  Computed  Performance  of  Generalized  Luneburg  Lenses 

Using  the  above  procedure,  the  results  for  a series  of  generalized 
Luneburg  lenses  having  pen^ct  refractive  index  profiles  and  various  types 
of  defective  index  profiles  are  presented  as  a series  of  figures  from 
Figure  15  through  71.  The  defects  have  been  introduced  as  changes  of 
the  best  of  n parameter.  The  data  has  been  calculated  for  the  s parameter 
value  of  2,  5 and  9 where  the  lens  f-number  is  one-half  of  s. 

These  figures  depict  the  overlay  waveguide  lens  profile  for  a 
standard  set  of  waveguide  parameters.  An  overlay  lens  of  Ta20g  having 
a refractive  index  of  n^  = 2.1  has  been  assumed  for  all  cases.  The  wave- 
guide effective  refractive  index  has  been  assumed  to  be  1.535  for  all  cases 
which  can  be  derived  from  a Corning  7059  thin  film  nw  = 1.565  having  a 
thickness  t = 1.0665y  using  0.9u  radiation.  The  substrate  in  all  cases 
is  Si02  having  a refractive  index  of  n$  = 1.47. 

Table  V provides  an  index  for  the  figures.  They  are  arranged  in 
groups  of  S=2,  S=5  and  S=9.  Each  group  starts  with  a perfect  profile 
with  Figures  15/16,  34/35  and  53/54,  respectively.  The  table  has  been 
arranged  for  symmetry  about  the  perfect  profile  starting  with  the  largest 
p-j  perturbation  normalized  to  the  perfect  profile  descending  through  the 
perfect  profile  to  an  equivalent  normalized  distortion.  Some  additional 
distortions  have  been  introduced  into  p4  and  p2.  The  serial  sequence 
starts  with  a perfect  profile,  with  normalized  steps  up  and  down  and  then 
the  results  for  an  order  of  magnitude  lesser  step.  A standardized  format 
has  been  employed  for  these  figures.  The  overlay  waveguide  thickness  pro- 
file is  included.  An  insert  also  indicates  the  ray  traces.  Ray  trace 
calculations  and  subsequent  data  use  only  0.8  of  the  lens  aperture.  The 
ray  intercept  in  the  image  plane  where  the  ray  intercept  is  defined  in 
Figure  9 is  included.  The  ray  intercept  utilizes  "x"  marks  correspond- 
ing to  individual  ray  trace  calculations  to  which  a smooth  curve  has  been 
fitted  for  subsequent  calculation.  Deviations  from  the  abscissae  repre- 
sent the  blur  width  in  the  image  plane  corresponding  to  the  ray  traces. 

Together  with  each  tabulation  of  Luneburg  lens  parameters  included 
in  these  figures  is  a companion  set  of  figures  containing  wavefront  phase 
error  and  the  corresponding  intensity  diffraction  pattern.  The  wavefront 
phase  error  and  intensity  diffraction  pattern  are  presented  for  the 
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TABLE  V.  INDEX  OF  FIGURES  FOR  PERFECT  AND  DISTORTED  LUNEBURG  LENS  PROFILES 
(Arranged  for  Synuietry  About  the  Perfect  Profile  with  Decreasing  P^) 


Lens 

Distortion 

S=2 

5s  5 

S=9 

1.05P1, 

1.35P2 

Fig.  29/30 

Fig.  48/49 

Fig.  67/68 

1 . 04P-J , 

io.op4 

Fig.  23/24 

Fig.  42/43 

Fig.  61/62 

1.0004P1 

. 1-1P4 

Fig.  27 

Fig.  46 

Fig.  65 

1.02Pn 

Fig.  17/18 

Fig.  36/37 

Fig.  55/56 

1.0002Pn 

Fig.  21 

Fig.  40 

Fig.  59 

Perfect  Profile 

Fig.  15/16 

Fig.  34/35 

Fig.  53/54 

0.9998Pn 

Fig.  22 

Fig.  41 

Fig.  60 

0.98Pn 

Fig.  19/20 

Fig.  38/39 

Fig.  57/58 

0.9996P1 

. -1.1P4 

Fig.  28 

Fig.  47 

Fig.  66 

0. 96P-| , 

-10.0P4 

Fig.  25/26 

Fig.  44/45 

Fig.  63/64 

0.9995P-, 

, 0.997P2 

Fig.  33 

Fig.  52 

Fig.  71 

0.95Pr 

0.7P2 

Fig.  31/32 

Fig.  50/51 

Fig.  69/70 

paraxial  focus.  Data  for  lenses  containing  a distorted  profile  include 
the  perfect  overlay  waveguide  lens  thickness  profile  together  with  the 
defective  profile.  For  some  cases  where  S=5  and  S=9,  the  distortion  is 
not  apparent,  however,  is  listed,  and  its  physical  effect  will  be  ap- 
parent by  referring  to  the  corresponding  S=2  case.  Inspection  of  Figures 
16,  35  and  54  for  S=2,  S=5  and  S=9,  respectively,  shows  the  wavefront 
phase  error  and  the  corresponding  intensity  diffraction  pattern.  The 
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Overlay  Lens  Thickness  (Microns) 


Figure  15.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 

Corresponding  Overlay  Waveguide  Lens  Thickness  Radial  Profile, 
Ray  Traces  through  this  Lens  and  the  Intercept  Error  in  the 
Image  Plane,  S = 2 (Perfect  Luneburg  Lens). 
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Overlay  Lens  Thickness  (Microns) 


Overlay  Waveguide  Thickness  Profile  GENERALIZED  LUNEBURG  LENS 


Perfect  Profile 


S = 2 


/ 

Defective  Profile 


.240401592 


P2  ■ -.0762450365 
P3  = 6.863523657xl0'3 
P4  = -5.247359951x10-3 
P5  = -4.239745156x10" 


Normalized  Radius  r/r 


Lens  Ray  Traces 


\ + 300xf/D 


Blur  Width 


Full  Aperture  of  Ray  Traces 
- 300  Xf/D 


Lens  Overlay  Ta205  n^  =2.1 
Waveguide  ne^  = 1.535 
Corning  7059  nw  = 1.565 
Thickness  t = 1.0665p,  X - 0.9y 
Substrate  Si0o  n$  = 1.47 


Ray  Trace  Intercepts 
in  Image  Plane 


Figure  17. 


Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 2 (1.02  P, , 1.02  P,,  1.02  P-, 
1.02  P4,  4.2856  P5).  1 1 i 


Figure  18.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  17. 
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Overlay  Lens  Thickness  (Microns) 


\ 

i 


Figure  19.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 2 (.98  P, , .98  P,,  .98  P-, 
.98  P4,  -2.2856  P5).  1 1 3 
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0 5 10  15  Radians 

Figure  21.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  17, 
S = 2 (1.0002  Plf  1.0002  P2,  1.0002  P3>  1.0002  P4>  1 .0328  P5). 
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Figure  22.  Wavefjront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  19, 
S = 2 (.9998  Pr  .9998  P2,  .9998  P3,  .9998  P4,  .9671  Pg). 
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Overlay  Waveguide  Thickness  Profile  GENERALIZED  LUNEBURG  LENS 


P]  - .2451153437 
P2  = -.07475003583 
P3  =■  6.728944762x10-3 
P4  = -5.14447054x10-2 
P5  = .0358834217 


Figure  23.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 


Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 2 (1.04  P, , 10  P.,  -36.2715  PR). 
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Figure  24.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  23. 
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Overlay  Waveguide  Thickness  Profile  GENERALIZED  LUNEBURG  LENS 

S = 2 = .2262603219 

P2  = -.07475003583 

/ p3  = 6.728944762x10 

Defective  Profile  / P4  = -6.14447054x10 

Perfect  Profile  p = -.0481509622 


Perfect  Profile 


Normalized  Radius  r/r 


Blur  Width 


+ 300xf/D 


Full  Aperture  of  Ray  Traces 


Lens  Overlay  Ta205  n-j  = 2.1 
Waveguide  ng^^  = 1.535 
Corning  7059  nw  = 1.565 
Thickness  t = 1.0665y,  X = 0.9y 
Substrate  SiO^  n$  = 1.47 


Ray  Trace  Intercepts 
in  Image  Plane 


Figure  25.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 2 (.96  P, , P0,  P,,  -10  P., 
48.6717  P5).  123  4 
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Phase  Error  in  Pupil  Plane 


GENERALIZED  LUNEBURG  LENS 

S = 2 Pj  = .2357821104 
?2  = -.07475003583 
P3  = 6 . 728944762x10"^ 
P4  = -5. 65891 7594xl0“3 
Pc  = -5.69127768x10"^ 


Aperture 


Diffraction  Pattern  T o dB 


0 5 10  15  Radians 

Figure  27.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  23, 
S * 2 (1.0004  Pr  1.1  P4,  .5752  Pg). 
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Figure  28.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  25, 
S = 2 (.9996  Pr  P2,  P 3,  -1.1  P4,  11.8249  Pg). 
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Overlay  Waveguide  Thickness  Profile  GENERALIZED  LUNEBURG  LENS 


Perfect  Profile 


Defective  Profile 


P]  = .2474722271 
P2  = -.1009125483 
P3  = 6. 728944 762x1 0"3 
P4  = -5. 14447054xl0"3 
Pc  = .0133888210 
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Lens  Overlay  Ta205  n]  =2.1 
Waveguide  ng^f  = 1.535 
Corning  7059  nw  = 1.565 
Thickness  t = 1.0665p,  x = 0.9p 
Substrate  Si02  ns  = 1.47 


Full  Aperture  of  Ray  Traces-'- 


Ray  Trace  Intercepts 
in  Image  Plane 


Figure  29.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 2 (1.05  P] , 1.35  P2,  -13.5336  Pg) . 
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Figure  30.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  29. 
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Figure  31.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 2 (.95  P, , .7  P9,  11.7550  Pc) 

Figure  32.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  31. 
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Figure  33.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  31 
S = 2 (.9995  P,,  .997  P0,  1.1075  P,). 


wavefront  phase  errors  for  the  perfect  profile  are  due  to  the  computational 

residues  arising  from  limitation  of  the  computational  time.  In  each  case, 

2 

the  intensity  diffraction  pattern  is  the  expected  sine  x function  which 
has  the  expected  single  lobe  and  side  lobes  due  to  use  of  a uniform  elimi- 
nation function. 

Some  of  the  computed  results  do  not  show  a single  principle  lobe 
type  of  diffraction  pattern.  This  is  due  to  the  introduction  of  exces- 
sive errors  in  the  waveguide  thickness  profile.  There  are  two  optical 
means  to  suppress  this  effect,  i.e.,  stopping  down  the  aperture  or  the 
introduction  of  a tapered  illumination  function.  The  computer  program  can 
be  altered  to  accommodate  various  aperture  distributions  and  their  limits. 

Some  of  the  figures  for  the  S=2  case  are  identical  to  those  contained  in 
the  preceding  Optical  Waveguide  Lenses  II  Report  AFAL-TR-77-60;  however, 
they  have  been  included  herein  again  for  comparison  with  the  S=5  and  S=9 
cases. 

Now  referring  to  Figure  17,  it  does  not  produce  a well-defined  single 
principle  lobe  diffraction  pattern  because  of  the  defective  profile  caused 
by  the  introduction  of  a 2%  increase  of  all  the  Pn  values  except  Pg.  Inspec- 
tion of  the  ray  intercept  clearly  indicates  that  the  computational  plane  is 
displaced  from  the  true  focusing  because  of  its  tilt  which  is  also  apparent 
from  the  curvature  of  the  wavefront  phase  error.  Using  a new  computational 
image  plane,  a much  better  well-defined  principle  lobe  diffraction  pattern 
can  be  obtained.  Figure  19  shows  the  results  for  a corresponding  2%  decrease 
of  all  Pn  values  except  Pg  leading  to  a distorted  diffraction  pattern  which 
can  be  refocused  by  shifting  of  the  image  plane  in  the  opposite  direction. 
Figures  21  and  22  also  show  these  same  effects  to  a much  lesser  degree  because 
the  distortion  introduced  is  0.02 % which  has  been  refocused  yielding  results 
nearly  identical  to  Figure  16.  Inspection  of  Figures  40  and  41  for  S*5  and 
Figures  59  and  60  for  S=9  show  similar  results;  however,  these  have  not 
required  refocusing.  This  is  because  the  central  lens  thickness  has  not 
been  changed  appreciably. 

It  is  also  apparent  that  stopping  down  the  lens  aperture  will  remove 
some  marginal  rays  and  refocusing  will  partially  negate  the  effect  of  defect 
tive  profiles.  This  can  be  observed  by  inspection  of  the  lens  ray  traces 
and  ray  intercepts  such  as  in  Figures  23/24,  25/26,  29/30,  and  31/32.  This 
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Figure  34.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 

Corresponding  Overlay  Waveguide  Lens  Thickness  Radial  Profile, 
Ray  Traces  through  this  Lens  and  the  Intercept  Error  in  the 
Image  Plane,  S = 5 (Perfect  Luneburg  Lens). 
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Figure  35.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  34. 
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Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 

‘ ^ , 1.02  P2,  1.02  r3, 


Error  in  the  Image  Plane,  S = 5 (1.02  P, , 1.02  P,,  1.02  P, 


1.02  P„ 
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Figure  37.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  36. 
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Figure  38.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 5 (.98  P, , .98  P9,  .98  Pv 
.98  P4,  -.7333  P5).  J 
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Figure  40.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  36, 
S = 5 (1.0002  Pr  1.0002  P2,  1.0002  1.0002  P4>  1.0173  Pg). 
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Figure  41.  Wavefront  Phase  Error  and  Corresponding  Intens-ity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  38 
S = 5 (.9998  Plf  .9998  POJ  .9998  P,,  .9998  P.,  .9826  Pc). 
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Figure  42.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 5 (1.04  P. , 10  P.,  4.0653  Pc). 
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Figure  43.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  42. 


67 


Overlay  Lens  Thickness  (Microns) 
Using  Ta?Oc  on  7059  Waveguide 


r 


Overlay  Waveguide  Thickness  Profile 


GENERALIZED  LUNEBURG  LENS 
S = 5 P1  = 0.08701439606 

P2  = -2.34300068  x 10'2 
P3  « -2.51100017  x 10'3 
P4  - 1.49192458  x 10'3 
P5  = 1.23968933  x 10'3 


Figure  44.  Tabulation  of  Generalized  LuneDurg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thicxness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 5 (.96  Pr  P2,  p -io  P4 
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Figure  45.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  In  Figure  44. 
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Figure  46.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  42, 
S = 5 (1.0004  Pr  1.1  P4,  1.0287  Pg). 
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Figure  48.  Tabulation  of  Generalized  Luneburg  Lens  Parameters 
with  Corresponding  Overlay  Waveguide  Lens  Thickness 
Radial  Profile,  Ray  Traces  through  this  Lens  and  the 
Intercept  Error  in  the  Image  Plane,  S = 5 (1.05  P,, 
1.35  P2,  P3,  P4,  3.9257  Pg). 
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Figure  49.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  48. 
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Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 5 (.95  P^  .7  P2,  4.3528  P5). 
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effect  is  more  pronounced  for  the  S = 5 and  S =9  lenses,  however,  requires 
a more  severe  limit  imposed  upon  the  aperture. 


As  the  S value  is  increased,  these  lenses  trend  toward  a more 
homogeneous  lens  where  its  refractive  power  is  concentrated 
near  the  edge.  This  dictates  that  the  mask  contour  employed  to  form 
these  lenses  should  accurately  conform  to  a circle  measured  to  a 
fraction  of  the  optical  wavelength. 

Inspection  of  the  S = 2 perfect  profile  and  various  types  of  defec- 
tive profiles  illustrate  that  the  defects  considered  represent  zonal 
errors.  For  Figures  17/18  and  19/20,  the  error  is  concentrated  in  a 
single  circle  zone.  For  Figures  23/24  and  25/26,  the  errors  are  con- 
centrated as  a circular  ring  located  near  r/rQ  = 0.5.  For  Figures  29/30 
and  31/32,  zonal  errors  are  introduced  of  opposite  effect  near  the  center 
and  near  r/rQ  =0.7.  These  types  of  errors  are  not  pictorially  apparent 
for  the  S = 5 and  9 sequence;  however,  the  pattern  has  been  preserved 
and  should  be  inferred  by  reference  to  Table  V. 

Subsequent  experimental  data  will  depict  ray  traces  through  various 
waveguide  Luneburg  lenses  and  diffraction  patterns.  The  resulting  diffrac- 
tion patterns  contain  a variety  of  defects  which  may  be  attributed  to 
spherical  aberration  in  the  coupling  of  convergent/divergent  beams  for 
measurements,  Mie  and  Rayleigh  scattering  in  and  on  the  waveguide  and 
lens,  and  imperfections  of  the  input  beam  due  to  scattering,  particularly 
at  the  input  coupling  prism  edge.  For  some  of  the  lens  data,  the  lens 
loss  is  not  insignificant.  Data  for  scattering  out  of  the  waveguide  above 
the  substrate  is  presented,  while  scattering  into  the  isolation  layer  and 
substrate  is  largely  absorbed.  All  of  these  real  effects  tend  to  mask 
the  inherent  lens  diffraction  pattern.  Thus,  achievement  of  near  diffraction- 
limited  performance  is  emphasized  for  attainment  of  the  appropriate  beam- 
width,  suppression  of  near-in  side  lobes  below  13-17  dB  and  suppression 
of  wide  angle  side  lobes  as  far  as  possible. 
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Overlay  Lens  Thickness  (Microns) 
Using  18305  on  7059  Waveguide 
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Overlay  Waveguide  Thickness  Profile  PERFECT  GENERALIZED  LUNEBURG  LENS 


Figure  53.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 

Corresponding  Overlay  Waveguide  Lens  Thickness  Radial  Profile, 
Ray  Traces  through  this  Lens  and  the  Intercept  Error  in  the 
Image  Plane,  S = 9 (Perfect  Luneburg  Lens). 
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Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 9 (1.02  P, , 1.02  P,,  1.02  Pv 
1.02  P4,  3.1015  P5).  1 J 
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Figure  57.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  In  the  Image  Plane,  S = 9 (.98  P,,  .98  P,,  .98  P-, 
.98  P4,  -.7398  P5).  1 1 3 
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Figure  58.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  57. 
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Figure  59.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  55, 
S = 9 (1.0002  Pr  1.0002  P2,  1.0002  P3>  1.0002  P4,  1.0174  Pg). 
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Figure  60.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  57, 
S = 9 (.9998  Pr  .9998  P2,  .9998  P3,  .9998  P 4,  .9826  Pg). 
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Figure  61.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 9 (1.04  P1 , 10  P4,  4.9788  P&). 
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Figure  62,  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  61. 
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Figure  63.  Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S = 9 (.96  P, , P9,  P,,  -10  P«, 
-2.7806  P5). 
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Figure  64.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  63. 
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Figure  65.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  61, 
S = 9 (1.0004  Pr  1.1  P4,  1.0387  Pg). 
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Figure  66.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  63 
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Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 
Error  in  the  Image  Plane,  S * 9 (1.05  P, , 1.35  P,,  Pv  Pd, 
-3.6564  Pc).  1 c J 
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Figure  69. 


Tabulation  of  Generalized  Luneburg  Lens  Parameters  with 
Corresponding  Overlay  Waveguide  Lens  Thickness  Radial 
Profile,  Ray  Traces  through  this  Lens  and  the  Intercept 


Error  In  the  Image  Plane,  S = 9 (.95  P1 , .7  P2,  4.1214  Pg). 


94 


o 3 10  15  Radians 


Figure  70.  Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  69. 


Figure  71 


Wavefront  Phase  Error  and  Corresponding  Intensity  Diffraction 
Pattern  for  Generalized  Luneburg  Lens  Delineated  in  Figure  69, 


S = 9 (.9995  P,,  .997  P9,  1.0312  P,). 
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6.  Computer  Simulation  of  the  Luneburg  Lens  Deposition  Mask 

The  preceding  analysis  has  accurately  defined  the  overlay  Luneburg 
lens  thickness  profile  necessary  to  produce  diffraction-limited  results 
and  the  effects  of  defective  profiles  as  exemplified  by  their  diffraction 
patterns. 

One  method  to  form  thln-fllm  waveguide  lenses  is  by  deposition 
through  an  appropriate  mask.  The  RF  sputter  process  has  been  employed 
throughout  for  this  purpose.  The  cross  section  of  the  RF  sputter  depo- 
sition sample  is  depicted  in  Figure  72.  Argon  Ions  are  accelerated  to- 
ward the  Ta  or  Nb  cathode  target  where  bombardment  releases  molecules 
that  are  oxidized  during  their  flight  until  collision.  In  the  sputter- 
ing environment  employed,  the  mean-free-pathlength  is  greater  than  the 
thickness  of  the  mask  and,  thus,  effects  of  collision  and  diffusion  are 
avoided.  A Lambertian  distribution  has  been  assumed  and  verified.  Thus, 
the  molecular  beam  arrives  at  the  mask  entrance  with  a predictable  angular 
momentum  distribution  A(0,v).  Where  the  mask  produces  a shadow,  inte- 
gration over  the  particle  velocity  range  provides  an  angular  distribution 
function  A^ ( 0 ) . The  distribution  function,  Ad(0),  was  experimentally 
measured  with  a knife-edge  mask  confirming  the  Lambertian  source.  Based 
upon  this  knowledge,  the  effect  of  a shadow  mask  upon  sputtering  has  been 
predicted  with  the  aid  of  a computer. 

The  shadow  mask  computation  is  relatively  simple  when  the  mask  has 
circular  symmetry.  The  coordinates  employed  to  compute  the  distribution 
are  illustrated  in  Figure  73.  For  an  arbitrary  point  PQ,  the  thickness 
distribution  becomes 


where 

0 = cos'1  D , R=  N/^-Xo^+^-yo)^2 

^(xrx0)2+(y,-y0)2 

and  (dxdy)  is  an  Infinitesimally  small  area  of  the  mask  entrance  window, 
and  S is  the  portion  of  the  entrance  window  which  is  visible  from  point  P. 
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The  distance,  D,  in  Figure 73  extends  from  the  plane  containing  the  top 
of  the  mask  to  the  substrate  surface. 

The  determination  of  the  effective  window  function,  S,  for  a 
mask  contour  of  a general  cylinder  is  also  illustrated  by  Figure  73 
The  edge  contour  is  conveniently  described  by  r = f (z)Q.  The  line  of 
sight  between  points  PQ  and  P1  is  given  by  the  solution  of 

*,y  ♦O'0-y,)x  - ■ o (48) 

z,x  ♦ (x0-x,)z  - z,x0  = 0 

which  defines  two  planes  each  parallel  to  the  x or  the  y axis  while  con- 
taining the  line  of  sight  PQpr  any  Point  on  line»  the  distance 
to  the  axis  of  rotation  of  the  mask  is,  therefore, 

r'  - ((VX')Z  - x»)2  +((y°'y1)Z  - )2  • (49) 

If  the  line  of  sight  is  not  interrupted  by  the  edge  of  the  shadow  mask, 
the  following  function  must  not  have  real  solution. 


(r'l2  - r2» 


o . |(vxi)z  - x0j  ^(vyi); 


f2(z) 


For  f(z)  described  by  quadratic  functions,  analytic  result  may  be  easily 
obtained.  For  f(z)  of  higher  order  functions,  numerical  evaluation  must 
be  employed.  Anytime  when  the  above  function  becomes  greater  than  or  equal 
to  zero,  it  means  that  interference  with  line  of  sight  is  occurring.  Thus, 
the  deposition  distribution  function  becomes 


T(x0.y0>  = II  a(9)cos9G  dxdy  ( 

FT 

where  W is  the  entrance  window  and  G is  a logic  function  which  is  one 
(or  zero)  when  the  line  of  sight  is  clear  (or  interrupted). 


z 


Sputtered  Molecule 
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Figure  74.  Coordinate  System  for  Knife-Eage  Shadow  Mask. 

The  anglular  distribution  of  sputtered  material  depends  upon  the 
ion  angle  of  arrival,  the  target  and  sputtering  system  geometry,  and  film 
growth  properties.  For  the  vacuum  environment  employed,  the  ions  experi- 
ence a mean-free  path  length  greater  than  the  shadowing  mask  thickness, 
while  the  target  source  extends  over  the  full  angular  range  of  the  shadow 
mask  acceptance  angle.  The  angular  distribution  of  sputtered  particles 
has  been  measured  with  a small  pinhole  and  with  a knife  edge  for  calibra- 
tion, and  has  been  found  to  confirm  the  assumed  cosine  or  Lambertian 
distribution. 

Assume  that  the  knife  edge  is  located  at  a distance  d above  the 
substrate,  as  shown  in  Figure  74,  and  that  the  deposited  film  thickness 
profile  is  B (x)  which  can  be  measured.  From  Figure  74  it  is  noted  that 
the  relationship  between  the  coordinate  x and  the  deposition  acceptance 
angle  4>(x)  is  x = d cot<f>.  Because  the  physical  deposition  rate  at  any 
point  is  actually  dependent  on  the  deposition  acceptance  angle  available 
to  that  point,  i.e., 


where  x is  the  deposition  time,  the  angular  deposition  rate  may  be 


tm 
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obtained  by  taking  the  derivative  of  B (x)  followed  by  a change  of 
variable,  i.e.. 


dB(x) 

dx 


-rd  A^(x) 


Ad(x) 


-Ad(e) 

x=d  cot© 


(53) 

(54) 


Thus,  from  the  above,  the  angular  deposition  rate  is  found  to  closely 
approximate  the  cosine  angular  distribution. 

With  the  computer  model  we  have  iterated  a sample  conical  mask  and 
computed  the  resulting  deposition  thickness  profile.  A cross  section  of 
the  conical  mask  and  its  descriptive  parameters  is  depicted  in  Figure  75 
(a).  The  effect  of  changing  the  conical  angle  is  depicted  in  Section  (b) 
of  the  Figure.  The  effect  of  the  mask  spacing  above  the  substrate  is 
depicted  in  Section  (c),  and  the  effect  of  mask  aperture  diameter  is  depicted 
in  Section  (d).  The  accuracy  of  these  computations  has  been  confirmed  by 
a stylus  probe,  including  measurements  of  the  slope  and  steps  for  a number 
of  shadow  masks. 

From  these  variational  effects  depicted  in  Figure  75,  a first 
iteration  approximation  to  a Luneburg  lens  thickness  profile  was  developed 
using  two  conical  sections  as  illustrated  by  the  cross  section  in  Figure  76. 
The  resulting  calculated  thickness  profile  is  indicated  together  with  an 
experimentally  measured  segment  near  the  edge  of  the  deposition.  Using  the 
preceding  computer  simulation  programs,  the  thickness  profile  is  converted 
to  a refractive  index  profile  with  ray  traces  through  the  lens  calculated 
to  find  the  ray  intercept  errors  at  the  focus.  Figure  77  depicts  the 
calculated  ray  intercepts  for  the  thickness  profile  of  Figure  76  derived 
from  the  conical  mask.  The  insert  illustrates  the  ray  traces.  The  experi- 
mental points  concerning  the  calculations  have  been  taken  from  the  measured 
ray  traces  photographed  for  Figure  78  using  a waveguide  Luneburg  lens  formed 
with  mask  SK-1.  The  mask  SK-1  is  only  a first  approximation  and  referring 
to  Figure  77,  can  only  be  used  out  to  0.5  of  its  radius.  Diffraction  data 
for  a number  of  waveguide  Luneburg  lenses  fabricated  with  mask  SK-1  are 
included  in  the  next  section. 
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(d) 


Figure  75.  Cross  Section  of  Mask  (a)  and  Calculated  Thin-Film  Deposition 
Profiles  for  Variable  Cone  Angle  (b)  Variable  Mask  Spacing 
Above  Substrate  (c)  and  for  Variable  Mask  Aperture  Diameter  (d). 
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Figure  78.  Measured  Ray  Traces  for  Waveguide  Luneburg  Lens  Formed  With 
Mask  SK-1 . 

A second  mask  (SK-2)  to  more  closely  approximate  the  ideal  lens 
thickness  profile  has  been  developed  which  employs  12  conical  segments  to 
shape  the  mask  edge.  The  resulting  computed  thickness  profile  is  depicted 
in  Figure  79  and  compared  with  the  profile  generated  from  mask  SK-1  and 
the  Luneburg  lens  ideal  thickness  profile.  A much  closer  approximation 
has  been  achieved.  The  computed  enlarged  cross  section  of  the  improved 
mask  SK-2  is  depicted  in  Figure  80.  Its  thickness  is  6 mm  and  its  dia- 
meter is  8 mm.  This  computer  generated  profile  was  transferred  directly 
for  fabrication  of  the  mask.  Diffraction  patterns  using  mask  SK-2  are 
presented  in  the  following  section  including  a lens  where  an  error  was 
made  by  deposition  through  the  wrong  side  of  the  mask.  These  masks  have 
been  scaled  to  a smaller  size;  however,  scaling  to  a larger  diameter  will 
require  a corresponding  increase  in  the  target  size  in  order  to  extend 
coverage  entirely  over  the  full  mask  aperture. 


Ideal  Thickness  Rwfile 

SKI-Mask  Profile  

SK2-Mask  Profile  


Figure  80.  Computed  Cross  Section  of  SK-2  Mask  Employed  for 
Optical  Waveguide  Luneburg  Lenses  Which  yield 
Diffraction-Limited  Results  (Mask  Thickness  6mm 
and  Diameter  8mm). 


Figure  79.  Sputtered  Deposition  Profiles  for  Waveguide  Luneburg  Lens 
Using  Computer-Generated  Masks. 


Figure  81  Calculated  Ray  Traces,  Wavefront  Phase  Errors  and  Intensity 


Diffraction  Patterns  for  Five  Overlay  Luneburg  Lens  Fractional 
Thicknesses  for  a Fixed  Focal  Length  Showing  Changes  Due  to 
Errors  of  Thickness  for  S=2. 

7.  Luneburg  Lens  Tolerance  Analysis 

The  preceding  computational  results  have  indicated  how  the  shape 
of  the  lens  thickness  profile  alters  the  diffraction  pattern.  The  lens- 
thickness  shape  is  determined  by  the  deposition  process  and,  specifically, 
the  mask-edge  shape.  A mask  which  can  produce  the  perfect  profile,  however, 
may  not  yield  ideal  results  because  of  either  waveguide  thickness  or  lens 
thickness  or  their  incorrect  refractive  indices.  To  determine  the  effect 
of  either  waveguide-thickness  or  lens-thickness  errors,  a set  of  ray  traces, 
wavefront  phase  errors,  and  diffraction  patterns  have  been  calculated  for 
S = 2,  5,  and  9 perfect  lens  thickness  profiles. 

The  lens  thickness  has  been  increased  and  has  been  decreased  to 
represent  an  incorrect  deposition  time.  The  resulting  ray  traces,  wave- 
front  phase  errors,  and  diffraction  patterns  are  depicted  in  the  collage 
of  Figures  81,  82,  and  83.  The  nominal  central  lens  thickness  of  2240A, 

1281 A,  and  1013&  yield  focal  lengths  of  1.040  cm,  2.499  cm,  and  4.01  cm  for 


106 


Lens  Ray  Traces 


Phase  Error  in  Pupil  Plane 


Figure  82.  Calculated  Ray  Traces,  Wavefront  Phase  Errors  and  Intensity 

Diffraction  Patterns  for  Five  Overlay  Luneburg  Lens  Fractional 
Thicknesses  for  a Fixed  Focal  Length  Showing  Changes  Due  to 
Errors  of  Thickness  for  S=5. 

the  S = 2,  5,  and  9 cases  respectively.  These  thickness  profiles  have 
been  increased  and  decreased  ± 100  and  ± 20oX  at  the  center  of  the  lens 
and  in  proportion  to  the  lens  thickness  out  to  the  edge.  Inspection  of 
the  resulting  diffraction  patterns  shows  that  an  error  of  thickness  intro- 
duces a slight  spherical  aberration  which  produces  a slight  increase  of 
the  adjacent  side  lobes  and  fill-in  of  the  first  null.  It  should  be  noted 
that  the  phase  error  introduced  by  a ± 200S  thickness  error  produces  an 
approximate  ff/8  wavefront  phase  error  for  all  of  the  S values  computed 
and,  thus,  similar  diffraction  patterns.  However,  for  a given  error  of 
thickness,  the  effect  upon  the  focal  length  of  S = 9 lenses  is  much  more 
than  that  of  S = 2 lens  design.  This  effect  upon  the  lens  focal  length 
is  apparent  in  Figure  84. 

Although  the  calculations  have  been  performed  by  varying  the  lens 
thickness,  comparable  results  are  obtained  by  varying  the  waveguide  thick- 
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Lens  Ray  Traces  Phase  Error  In  Pupil  Plane 


Figure  83.  Calculated  Ray  Traces,  Wavefront  Phase  Errors  and  Intensity 

Diffraction  Patterns  for  Five  Overlay  Luneburg  Lens  Fractional 
Thicknesses  for  a Fixed  Focal  Length  Showing  Changes  Due  to 
Errors  of  Thickness  for  S=9. 

ness.  A nominal  value  for  the  Ta,0c  deposition  rate  employed  is  approxi- 
mately  17A/minute  and,  thus,  for  a 100a  error,  the  deposition  tolerance 
becomes  approximately  6 minutes  so  as  not  to  affect  the  diffraction  pattern, 
however,  the  deposition  time  must  be  controlled  more  accurately  to  achieve 
the  correct  focal  length. 
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Luneburg  Lens  Focal  Length  (cm) 


Figure  84.  Calculated  Waveguide  Luneburg  Lens  Focal  Length  for  S=2, 
S=5  and  S=9  Perfect  Profiles  Where  a Fractional  Thickness 
Error  Has  Been  Introduced  to  Simulate  the  Effect  of  an 
Error  of  the  Deposition  Time. 


c.  EXPERIMENTAL  WAVEGUIDE  LUNEBURG  LENSES 

All  the  waveguide  lenses  reported  herein  have  employed  Ta205  as 
the  highest  refractive  index  material  used  for  forming  the  lenses.  All 
the  waveguides  have  employed  Corning  7059  glass  on  thermally-grown  Si02 
on  silicon  as  a substrate.  All  of  the  lenses  have  utilized  the  dense 
overlay  and  have  been  deposited  by  an  RF-sputtering  technique  starting 
with  metallic  Ta  and  using  a reactive  process  with  Ar-02  to  form  the 
deposition  through  a circular  symmetric  mask  with  an  appropriately  shaped 
edge. 


Optical  Waveguide  and  Waveguide  Lens  Characterization 


The  performance  of  conventional  optical  elements  is  determined 
by  various  procedures  such  as  measurement  of  the  diffraction  pattern 
for  an  accurately  known  source  or  by  measurement  of  modulation  trans- 
fer function  for  periodic  sources  having  variable  period.  In  each 
case,  access  to  the  image  plane  is  required  for  probing  of  the  field 
intensity  in  the  image  plane.  For  optical  waveguide  lenses,  this 
technique  is  not  practical  because  the  lens  is  immersed  in  an  ambient 
dielectric  waveguide. 


When  the  ambient  dielectric  waveguide  exhibits  uniform  scatter- 
ing, it  is  possible  to  observe  the  field  distribution  spatially  through 
the  scattering  mechanism.  However,  in  general  the  lens  and  waveguide 
which  exhibit  sufficient  scattering  for  this  type  of  observation  also 
represents  an  unsatisfactory  lens  for  use  in  an  optical  Fourier  trans- 
form. 


The  technique  employed  to  measure  the  diffraction  pattern  of 
waveguide  lenses  has  utilized  a reimaging  procedure  with  several  types 
of  coupling.  A schematic  of  the  reimaging  apparatus  is  depicted  in 
Figure  85.  Prism  couplers  are  employed  to  launch  beams  within  the  thin- 
film  waveguide  which  traverses  the  lens  under  investigation.  Four  types 
of  coupling  have  been  employed  to  extract  the  field  which  is  relayed 
and  reimaged  with  magnification  where  it  is  possible  to  probe  the  field 
distribution  in  this  latter  plane.  A narrow-slit  aperture  with  an  opti- 
cal fiber  has  been  employed  to  probe  the  field  distribution.  The  slit- 
fiber  probe  is  scanned  in  any  of  the  three  major  axes  providing  a signal 
to  drive  a recorder  while  the  optical  signal  is  detected,  amplified  and 
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Substrate 


Figure  85.  Schematic  of  Reimaging  Technique  to  Characterize  Thin- 
Film  Waveguide  Lens  Performance. 

compressed  to  give  a logarithmic  scale  for  the  recorder.  All  experi- 
mental diffraction  patterns  reproduced  herein  have  been  measured  in  this 
manner.  The  physical  apparatus  showing  the  input  beam  alignment  and 
coupling  to  the  waveguide-lens  substrate  mounting  fixture,  followed  by 
the  relay  lens  and  the  slit-optical  fiber  probe,  mechanical  scanning 
fixture,  are  depicted  in  Figure  86  with  the  electronics  and  recorder. 

The  relay  lens  is  a Tropel  diffraction-limited  Fourier  transform  system 
having  F/3  aperture  with  a focal  length  of  221mm. 

Various  methods  of  coupling  from  the  waveguide  lens  substrate  to 
the  relay  lens  have  been  employed.  Most  of  the  experiments  have  utilized 
conventional  prism  couplers  with  special  attention  given  to  obtain  a 
sharp  defect-free  edge.  Another  type  of  coupling  which  eliminates  the 
prism  utilizes  a cleaved  exit  edge  for  the  waveguide.  A good  cleaved 
edge  close  to  the  lens  focus  permits  the  use  of  a high-pewer  microscope 
objective  as  the  relay  lens.  A close-up  view  of  such  a configuration 
using  an  input  coupling  prism  on  a waveguide  on  a thermally-oxidized 
silicon  substrate  with  a Ta20g  overlay  Luneburg  lens  is  depicted  in 
Figure87.  A scanning  slit  fiber  optical  probe  is  again  used  to  measure 
the  image  distribution. 
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Figure  86.  Photograph  of  Apparatus  for  Measurement  of  Thin-Film 
Waveguide  Performance  Using  Tropel  F/3  Fourier 
Transform  Lens  Assembly  (FL  = 221mm)  and  Scanning 
Slit-Fiber  Optical  Probe. 


Figure  87.  Close-up  Photograph  of  Apparatus  to  Measure  the  Intensity 
Distribution  From  Waveguide  and  Lenses  by  Cleavage  Through 
the  Waveguide  and  Reimaging  Using  a Zeiss  40X  Microscope 
Objective. 
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The  use  of  a plane  parallel  optical  plate  introduced  into  a 
convergent  or  divergent  beam  introduces  spherical  aberration;  and  thus, 
the  use  of  a prism  coupler  to  extract  radiation  from  waveguide  contain- 
ing lenses  introduces  an  error  into  the  measurement  due  to  this  added 
spherical  aberration.  Several  methods  were  investigated  to  minimize  or 
eliminate  this  effect.  The  most  practical  technique  proved  to  be  posi- 
tioning of  the  coupler  so  that  the  coupling  section  straddles  the  focus 
equally  and  thus,  in  this  position,  cancels  the  spherical  aberration. 

Two  other  techniques  investigated  to  reduce  the  spherical  aberration  due 
to  the  coupler  are  illustrated  in  Figures  88  and  89.  Thin  spinel  wafers 
edge  mounted  are  shown  in  Figure  88.  Effective  coupling  is  obtained  for 
0.5mm  thick  wafers.  An  aplanatic  prism  with  a spherical  hypotenuse  with 
its  radial  center  located  on  the  coupling  edge  can  also  eliminate  spheri- 
cal aberration.  Such  an  aplanatic  prism  is  depicted  in  Figure89.  This 
concept  is  identical  to  an  immersion  microscope  objective.  However,  it 
must  be  accurately  positioned  so  that  the  spherical  surface  matches  the 
emerging  optical  wavefront.  Accurately  positioning  of  the  aplanatic 
prism  is  difficult  to  achieve  routinely.  Elimination  of  the  prism  using 
a cleaved  edge  of  the  waveguide  also  requires  that  the  cleavage  be 
accurately  positioned  at  the  focus  of  the  lens  to  eliminate  a spherical 
aberration  introduced  at  the  waveguide  exit  refraction. 

Lens  diffraction  patterns  are  the  convolution  of  the  "true"  lens 
diffraction  pattern  and  the  observational  instrumentation.  For  the  case 
at  hand,  the  relay  lens  has  been  employed  to  provide  magnification  for  a 
mechanical  scanning  of  a slit  across  the  reimaging  diffraction  pattern. 

As  the  lens  performance  is  improved,  the  reimaging  lens  must  exhibit 
comparable  performance  or  superior  performance.  When  couplers  are  em- 
ployed, their  distortions  are  also  convolved. 

Most  of  the  measurements  reported  herein  employed  a Trope!  F/3  lens 
designed  for  the  Fourier  transform  as  the  relay  lens.  Further,  some  of 
the  measurements  have  been  performed  using  a Zeiss  40X  microscope  objective 
with  a numerical  aperture  of  0.67. 

The  effect  of  this  lens  upon  the  measurements  will  be  considered 
herein.  For  most  of  the  measurements  of  F/4  Luneburg  lens,  the  Tropel  re- 
imaging lens  becomes  a limiting  factor,  while  the  Zeiss  lens  is  not.  How- 
ever, it  requires  a cleaved  exit  face  accurately  positioned  at  the  focus. 
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Figure  88.  Thin  Spinel  Wafers  Edge  Coupled  to  Corning  7059  Waveguide  on 

Thermally-Oxidized  Silicon  Substrate  Used  to  Minimize  Effective 
Prism  Coupler  Thickness. 


Figure  89.  An  Aplanatic  Prism  Coupler  Fabricated  from  SF-6  Glass. 
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The  Tropel  lens  is  designed  for  a magnification  of  2.  Thus,  for  this 
Tropel  lens,  the  effective  f-number  becomes 


Feff  - fO+fi>  ■ «-5  fe5) 

The  effective  f-number  is  the  actual  f-number  of  the  Tropel  lens 
in  the  reimaging  set  up.  Because  of  the  index  of  refraction  of  the  optical 
waveguide,  which  is  1.5,  the  F/4  Luneburg  lens  appears  as  an  F/2.7  lens  in 
air  and  is  severely  truncated  by  the  F/4. 5 Tropel  lens  reimaging  system. 


The  effect  of  truncation  on  the  beamwidth  measurement  for  a Gaussian 
optical  system  will  be  discussed  as  follows.  For  a Gaussian  optical-input 

O 

beam  with  e wid£h  of  "a",  the  optical  input  amplitude  function  can  be 
described  by  e~^x  ^a2.  The  optical  amplitude  distribution  at  the  back 
focal  plane  of  a planar  lens  is  obtainable  by  taking  a Fourier  transformation 
of  the  input  amplitude  function.  The  optical  intensity  distribution  at  back 
focal  plane  is  therefore. 


where  k is  the  optical  propagation  constant  x1  is  the  transverse 
coordinate  in  the  back  focal  plane,  and  F is  the  focal  length  of  a planar 
lens.  The  integral  within  the  bracket  is  a constant  related  to  the 
total  optical  power  of  the  input  optical  beam.  It  is  ^nvenient  to  de- 
fine the  f-number  of  this  focusing  system  as 

f = A F/a  (57) 
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The  beamwidth  (i.e.,  spot  width)  at  the  back  focal  plane  can  be  obtained 
from  Equation  (56).  Converting  Equation  (56)  to  the  dB  notations,  we 


-j-  = -8.69  x-j2  dB 

° (*A2  (5f 

' •2,'44(xr)  dB 

_2 

The  e beamwidth  is  twice  the  value  of  x-j  which  makes  I/IQ  = 8.79  dB, 
while  the  half  intensity  spot  width  is  twice  the  value  of  x-|  which  makes 
I/IQ  = -3  dB.  Thus 


_2 

e spot  width 


1.27  XF/a 


[3  dB  spot  width  : 0.75  XF/a 

In  practical  optical  imaging  system,  all  the  lenses  have  limited 
aperture  which  truncates  the  Gaussian  optical  profile  and  causes  the  beam 
width  to  increase.  Mathematically,  the  focal  spot  for  truncated  Gaussian 
beam  has  an  intensity  distribution  of 


A * sine 


($)’ 


where  the  sine  x function  denotes  sinx/x,  F'  is  the  focal  length  of  the 
reimaging  lens,  D is  the  aperture  of  the  reimaging  lens,  xj1  is  the  coor- 
dinate in  the  reimaged  plane,  and  f'  A F'/D  is  the  effective  f-number  of 
the  reimaging  lens.  To  distinguish  f and  f',  we  shall  call  f the  "optical 
beam  f-number"  and  f'  the  "truncation  f-number".  The  symbol  "*"  denotes 
the  convolution  operation.  Equation  (60)  can  be  readily  evaluated  with 
a computer  program  for  the  simulation  of  actual  coherent  optical  imaging 
system. 

For  systems  involving  multiple  truncations,  a chain  of  convolution 
operations  may  be  performed.  Because  both  the  Gaussian  function  A and 
the  sine  function  resemble  a delta-function  when  plotted  at  enlarged  scales, 
it  is  obvious  that  (1)  with  little  truncation^ i « f,  It»  | A | 2 and  (2) 
with  extreme  truncation,  f « f',  1^.  « |sinc  ^p-  |2.  That  is,  the  effect 
of  truncation  can  be  ignored  when  f'  » f.  For  f'  « f,  the  truncated  beam 
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resembles  an  uniformly  illuminated  aperture  of  D focusing  at  the  trunca- 
tion f-number  f'.  The  3 dB  beamwidth  for  J:j^e  extremely  truncated  Gaussian 
beam  is  therefore  the  3 dB  width  of  sine  and  is  0.885  Af'.  Figure 90 

shows  the  theoretical  beamwidth  size  evaluation  capability  of  the  two 
testing  lens  systems  employed.  The  straight  line  at  a slope  is  the  theo- 
retical  -3  dB  focal  beamwidth  for  untruncated  Gaussian  Beam.  The  -3  dB 
beamwidth  is  used  because  it  is  less  confusing  for  non-Gaussian  beams  as 
well  as  truncated  Gaussian  beams.  The  horizontal  lines  are  the  limiting 
beamwidth  for  the  two  reimaging  lenses  used  when  they  severely  truncate  a 
Gaussian  beam.  The  curved  lines  are  the  theoretical  expectation  for  the 
3 dB  focal  beamwidth  measured  with  the  two  systems.  For  simplicity,  the 
optical  waveguide  lens  is  assumed  not  truncating  the  Gaussian  beam  at  all, 
because  the  truncation  ratio  for  the  waveguide  Luneburg  lens  involved  in 
this  report  is  f/f^  = 4 which  falls  close  to  the  "little  truncation" 
category.  It  is  noted  that  the  convolution  effect  on  the  measured  results, 
i.e.,  the  curved  lines,  can  be  approximated  by  the  commonly  used  square 
law  given  by  Bracewell 

2 2 . 2 ,C1 , 

“an  “a  * a l (6,) 

where  a's  are  the  -3  dB  beamwidth  for  the  convolving  systems.  Eq.  (61 ) 
is  much  easier  to  use  than  Eq.  (56)  and  may  also  be  extended  to  include 
other  system  aberrations,  and  truncation  induced  beam  broadening  effects. 

For  the  beamwidth  measurement  system  using  the  Tropel  reimaging  lens, 
the  -3  dB  beamwidth  for  a uniformly  illuminated  line  aperture  is 

a£=  0.885  (0.6328)  4.5  = 2.52  microns  (62) 

assuming  the  Tropel  lens  is  diffraction  limited. 

The  -3  dB  beamwidth  for  an  F/4.1  waveguide  Luneburg  lens  in  a 
waveguide  having  an  effective  refractive  index  of  1.5,  is 

Og  = 0.75  (0.6328/1.5)  4.1  = 1.3  microns  (63) 

Applying  Snell's  law  for  refraction  into  air,  the  equivalent  beamwidth 
becomes 

a = 0.75  (0.6328)  4. 1/1. 5 = 1.3  microns  (64) 

a 
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Using  Eq.  (61 ) the  measured  -3  db  beamwidth  for  a diffraction  limited 
waveguide  Luneburg  lens  is  expected  to  be 

Vi  = (1-3)2  + (2-52)2  = (2*84)2 

where  our  measured  results  are  between  3.0  microns  and  3.3  microns.  We  have 
tried  to  verify  the  value  of  a ^ by  placing  a 1 micron  diameter  pin-hole  in 
front  of  the  Tropel  lens  and  have  experimentally  measured  = 2.9  micron, 

which  is  15%  greater  than  the  diffraction  limited  value.  Using  the  measured 
instead,  the  expected  beamwidth  becomes 

aa**=  0.3)2  + (2.9)2  = (3.18)2  (65) 

which  is  in  good  agreement  with  the  actual  experimental  result.  It  is 
concluded  that  the  F/4.1  Luneburg  lens  performance  is  near  the  diffraction 
limit  and,  thus,  the  true  beamwidth  is  approximately  1.3  microns  in  wave- 
guide. From  this  discussion,  it  is  clear  that  the  reimaging  lens  is  limiting 
the  true  measurement  of  the  waveguide  lens  performance. 

To  improve  our  confidence  of  the  waveguide  Luneburg  lens  performance, 
the  second  measurement  system  using  a Zeiss  40  X F/0.67  microscope  objective 
lens  as  the  reimaging  lens  and  a (100)  cleavage  near  the  focal  point  instead 
of  prism  coupler,  was  constructed.  Such  system  will  provide  a system 
limited  beamwidth,  o£,  of 

= 0.885  (0.6328)  0.67  = 0.38  micron  (66) 

An  F/5.5  lens  made  of  the  same  sputtering  mask  with  different 
sputtering  time  was  evaluated  with  this  system.  The  theoretical  beamwidth 
is  expected  to  be 

o = cr  = 0.75  (0.6328)  5. 5/1. 5 = 1.74  micron  (67) 

9 a 

The  expected  measurement  result  should  be 

aa*£  = (1-74>2  + (0'38)2  = 0-78)2  (68) 

whereas  our  measured  result  is  1.65  micron. 

This  measured  result  is  7%  less  than  the  calculated  result  probably 
due  to  the  fact  that  the  illumination  beam  could  be  slightly  larger  than 
the  2mm  assumed  which  would  lead  to  a reduced  f-number.  Thus,  it  is 
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Figure  90.  Measured  Diffraction  Pattern  Beamwidth  (3dB)  as  a Function 
of  the  f-Number  in  Waveguide,  Limited  by  Convolution  With 
the  Reimaging  Lens  Showing  Experimental  Results. 

verified  that  our  Luneburg  lenses  are  actually  diffraction  limited.  Other 
measurements  with  waveguide  Luneburg  lenses  of  larger  f-number  have  no 
problem  with  either  of  the  measurement  systems  as  shown  in  Figure  90. 
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Optical  Waveguide  Luneburg  Lens  Performance 


The  diffraction  performance  of  waveguide  Luneburg  lenses  fabricated 
from  masks  SK-1  and  SK-2  are  delineated  herein.  Nearly  all  of  the  measure- 
ments have  employed  a Tropel  F/3  reimaging  lens  at  a magnification  of  2x 
which  results  in  the  Tropel  lens  operating  as  F/4.5.  Several  lenses  have 
been  evaluated  with  a Zeiss  40x  microscope  objective  which  is  used  with 
F/0.67.  The  latter  lens  f-number  does  not  impact  the  observed  diffraction 

v 

pattern,  whereas  the  Tropel  reimaging  lens  significantly  effects  the 
observations. 

The  waveguide  Luneburg  lenses  have  employed  a 2 mm  Gaussian  beam, 
whereas  the  lens  diameter  for  most  cases  is  8 mm.  The  lens  masks  were 
designed  for  F/4.  Lenses  fabricated  using  the  SK-1  mask  with  thicknesses 
for  longer  focal  lengths  also  yielded  near  diffraction-limited  results. 

Lenses  fabricated  using  the  SK-2  mask  yielded  diffraction-limited 
results  with  much  lower  side  lobe  levels,  partially  due  to  a closer  approxi- 
mation to  the  ideal  in  this  lens  profile,  and  partially  due  to  improved 
depositions  that  reduce  scattering.  The  diffraction  patterns  and  photo- 
graphs in  Figures  91  ,92  , 93  , 94  , 95  , 96  , and  97  are  the  result  of  using 
mask  SK-1.  Each  of  the  diffraction  patterns  are  presented  in  a logarithmic 
scale  and  a linear  transverse  dimension  in  microns  measured  at  the  focal 
point.  The  beam  widths  are  indicated.  Representative  patterns  for  various 
focal  lengths  are  included  with  some  corresponding  photographs  of  radiation 
traversing  the  lens  and  waveguide  for  both  low-loss  and  moderate  loss 
waveguides  (0.2  = 1.2  dB/cm).  Figure  91  represents  a diffraction  pattern 
which  is  essentially  diffraction-limited  to  the  degree  that  it  can  be 
measured  with  the  Tropel  lens.  The  diffraction  pattern  of  Figure  92  has 
an  inferior  beamwidth.  Both  of  the  patterns  have  comparatively  high 
shoulder  lobes  which  indicate  that  the  lens  thickness  profile  is  not  fully 
corrected.  The  operation  of  the  lens  measured  in  Figure  92  is  shown  in 
Figure 93  used  as  a single  lens  and  used  as  a beam  expander  in  reverse. 

The  scattered  radiation  and  the  throughput  at  the  edge  of  the  wafer  are 
apparent. 

The  diffraction  patterns  for  Figures  94  and  95  apply  to  F/22  and 
F/23  lenses  with  focal  lengths  of  44  and  46  mm.  The  losses  of  these 
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Intensity  (dB ) 


Wafer:  SK- 1 - 5th 


H I-  3.1  W 


Luneburg  Lens  Date:  10/15/77 

Waveguide:  0.45u  7059  Wafer  I:  SK-1 -5th 
Hask:  Ski  Lens:  F/4,  SmmFL 

Isolation:  ~4.0u  510,  Substrate:  SI 
Resistivity:  10*  t 103u-c» 

Orientation:  [1,0,0] 

Lens  Diameter:  Snn 


Figure  91.  Measured  Diffraction  Pattern  of  Lens  Using  Mask  SK-1  and 
Tropel  F/3  Reimaging  Lens  Where  Lens:  F/4,  FL  = 8mm, 
Wafer:  SK-1 -5th 


Figure  93  Typical  Scattering  Traces  of  6328&  Radiation  Traversing  Ta^O 
Luneburg  Lenses  on  Corning  7059  Waveguide  on  Thermally- 
Oxidized  Silicon  Substrate.  Single  Luneburg  Lens  (a).  Beam 
Expander  (b),  and  Lens  Layout  as  Newton  Rings  (c). 


Luneburg  Lens  Date:  12/8/77 

Waveguide:  0.45m  7059  Wafer  I:  9 
Mask:  SKI  Lens:  F/22*  44mn  FL 
Isolation:  ~4.0m  S10?  Substrate: 
Resistivity:  102  - 103tt-cm 

Orientation:  [1*0,0] 

Lens  Diameter:  8nin 
Waveguide  Loss:  1.123  dB/cm 


Figure  94 


Measured  Diffraction  Pattern  of  Luneburg  Lens  Fabricated 
Using  Mask  SK-1  Using  Tropel  F/3  Reimaging  Lens  Where  Lens: 
F/22  and  FL  = 44mm  With  1.2dB/cm  Waveguide  Loss. 


Luneburg  Lens  Date:  12/7/77 

Waveguide:  0.45u  7059  Wafer  I:  8 

Mask:  SKI  Lens.  F/23,  46m  FL 
Isolation:  ~4.0u  SIO?  Substrate:  SI 

Resistivity:  10z  - 103fl-cm 

Orientation:  [1.0,0] 

Lens  Diameter:  8nm 
Waveguide  Loss:  1.096  dB/cm 
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Figure  95 


Measured  Diffraction  Pattern  of  Luneburg  Lens  Fabricated 
Using  Mask  SK-1  Using  Tropel  F/3  Reimaging  Lens  Where  Lens: 
F/23  and  FL  = 46mm  With  1.1 dB/cm  Waveguide  Loss. 


0 


Figure  96  Measured  Diffraction  Pattern  of  Lens  Using  Mask  SK-1  with  a 
Cleaved  Reflector  Using  Tropel  F/3  Reimaging  Lens  Where  Lens: 
F/15 , FL  = 30  mm. 


lenses  and  waveguides  are  comparatively  high  and  thus  the  wide  angle  side 
lobes  are  correspondingly  comparatively  high.  The  operation  of  an  F/15 
lens  with  a focal  length  of  30  mm  reflected  from  a cleaved  edge  is  depicted 
in  Figure  96  . Its  comparatively  high  side  lobes,  due  to  the  use  of  the 
SK-1  mask,  are  comparable  to  Figure  91  . The  transmission  streak  through 
the  lens  of  Figure  95  is  illustrated  in  Figure  97  , section  (a)  where  its 
throughput  distribution  at  the  edge  of  the  wafer  is  apparent.  The  conver- 
gence of  the  beam  to  the  focus  is  not  apparent  because  of  the  diffused 
scattering  which  corresponds  to  the  high  side  lobe  level  in  Figure  94  . 

The  transmission  streak  through  the  lens  of  Figure  96  with  the  reflector 
converging  to  a focus  diverging  and  reflecting  the  second  time  appears 
in  Figure  97  , section  (b).  The  reflectors  utilized  the  total  internal 
reflection  phenomena  at  the  cleaved  edge. 

Diffraction  patterns  in  Figures  98,99  ,100,101  ,102,103  , and  104  are 
the  result  of  using  mask  SK-2.  The  measured  wide-angle  diffraction  pat- 
tern of  a diffraction-limited  Luneburg  lens  using  the  Tropel  F/3  (oper- 
ating at  F/4.5)  reimaging  lens  is  depicted  in  Figure  98  where  the  lens: 

F/4  and  FL  = 8mm.  Its  near-inside  lobes  are  below  17dB  and  the  decay 
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Figure  97  Propagation  Streak  Through  a Luneburg  F/25  Lens  Fabricated 
With  an  SK-1  Mask  With  Throughput  Distribution  at  the  Edge 
of  the  Wafer  in  the  Divergent  Region  (a)  and  Propagation 
Streak  Through  an  F/15  Lens  Fabricated  With  an  SK-1  Mask 
With  Two  Total  Internal  Reflections  From  Cleaved  Edges  (b). 

to  30dB  within  approximately  200  beamwidths.  Beyond  that  range,  the 
recorded  signal  is  due  to  the  instrumentation  electronic  noise.  Another 
lens  fabricated  using  mask  SK-2  with  F/5.5  which  yields  diffraction- 
limited  results  as  measured  with  the  Zeiss  40X  objective  is  depicted  in 
Figure  99  . Note  that  its  beamwidth  is  1.65  microns  and  that  its  near- 
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Figure  93  . Wide  Angle  Measured  Diffraction  Pattern  of  Waveguide  Luneburg  Lens  Using  Mask 
SK-2  Using  Trope!  F/3  Reimaging  Lens  Where  Lens:  F/4  and  FL  = 8 mm. 


Figure  99  Measured  Diffraction  Pattern  of  Waveguide  Luneburg  Lens  Fabricated  Using  Mask  SK-2  Using 
a Cleaved  Exit  Face  at  the  Focus  and  Zeiss  X40,  F/. 67  Microscope  Objective  Showing  Dif- 
fraction-Limited Performance  Where  Lens:  F/5.5  and  FL  = 11mm  (Expanded  Scale  Insert). 


Figure  100  Measured  Diffraction  Pattern  of  Waveguide  Luneburg  Lens  Fabricated 
Using  Mask  SK-2  Using  Tropel  F/3  Reimaging  Lens  Where  Lens:  F/4 
and  FL  = 8 mm  Positioned  at  the  Best  Focus. 

inside  lobes  are  less  than  15dB,  however,  they  decay  to  approximately 
23dB  in  approximately  25  beamwidths  where  the  instrumentation  electronic 
noise  dominates  the  recording.  The  results  of  Figures  98  and  99  repre- 
sent diffraction-limited  lenses  to  the  degree  that  the  instrumentation  is 
capable  of  measuring. 

Figures  100  .101  » and  102  depict  the  diffraction  pattern  of 
another  waveguide  Luneburg  lens  fabricated  using  mask  SK-2  and  using  the 
Tropel  F/3  reimaging  lens  where  F/4  and  FL  = 8mm.  Figure  100  represents 
the  best  focus  while  Figure  101  is  measured  62  microns  outside  of  this 
focus  and  Figure  102  is  measured  62  microns  inside  of  this  focus.  These 
three  patterns  show  how  the  diffraction  pattern  varies  to  indicate  its 
depth  of  field. 


Intensity  (dB)  E Intensity  (dB) 


e 101  Same  as  Figure  100  Except  Measured  62  Microns  Outside  of  the  Focus 


Figure  102  Same  as  Figure  100 Except  Measured  62  Microns  Inside  of  the  Focus. 


The  diffraction  pattern  of  Figure  103  represents  a scaled 
increase  of  mask  SK-2  by  a factor  of  2.  It  shows  the  diffraction 
pattern  for  both  the  TE  and  TM  modes  where  the  correct  lens  design 
applies  to  the  TM  mode.  The  comparatively  high  side  lobes,  including 
a shoulder  at  -5dB,  is  due  to  the  lack  of  circular  symmetry  of  the 
deposition.  This  is  confirmed  by  the  Newton  rings  indicating  the 
thickness.  It  is  the  result  of  a slight  offset  from  the  center  of 
symmetry  in  the  sputtering  chamber. 

Another  mistake  type  of  diffraction  pattern  is  shown  in  Figure 
104  . In  this  case,  the  mask,  SK-2,  was  employed  upside  down.  In- 
spection of  the  mask  profile  in  Figurel04  shows  that  its  contour  edge 
is  not  symmetrical.  This  type  of  distortion  can  also  be  trace  to  the 
profile  errors  considered  in  Section  II-B-5. 

D.  WAVEGUIDE-LENS  SCATTERING 

Mie  and  Rayleigh  scattering  in  optical  waveguide  and  waveguide  lenses 
is  an  important  consideration  that  may  limit  the  application  of  integrated 
optical  circuits.  Various  techniques  have  been  investigated  to  reduce 
scattering  with  varying  degrees  of  success.  Processing  methods  to  obtain 
consistent  reproducibility  of  low-scattering  waveguide  has  been  included. 
Consideration  of  protective  superstrate  has  also  been  included.  This  in- 
vestigation has  extended  beyond  the  typical  materials  employed  for  inte- 
grated optics  heretofore.  This  effort  has  largely  been  experimental  in 
nature  and  oriented  for  process  optimization. 

Consideration  of  the  effect  of  scattering  upon  the  dynamic  range 
of  an  integrated  optical  Fourier  transform  circuit  is  treated  fully  in 
Appendix  C.  The  best  measured  scattering  is  consistent  with  the  allow- 
able lens  and  waveguide  scattering  for  application  to  an  optical  Fourier 
transform.  Requirements  for  waveguide  scattering  are  an  order  of  magni- 
tude more  stringent  because  of  its  close  proximity  to  the  photodetector 
array.  Waveguides  considered  have  included  Corning  7059  ( purified 
Si02:Ba0),  Si02  (P205),  Si304:Si02,  Ta205  and  Nb205.  Superstrates  in- 
cluded various  photoresists  and  PMMA.  Thermally-oxidized  silicon  sub- 
strates have  been  employed  throughout. 
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Figure  103  Measured  Diffraction  Pattern  of  Luneburg  Lens  Fabricated 

Using  Mask  SK-3  Using  Tropel  F/3  Reimaging  Lens  Where  Lens: 
F/I2.2and  FL  = 25mm  Showing  the  Results  for  Both  TE  and  TM 
Modes. 


Figure  105  Traveling  Fiber  Optical  Probe  for  Measuring  Longitudinal 
and  Transverse  Scattered  Radiation  Above  the  Substrate. 


As  a diagnostic  tool  for  measurement  of  waveguide  and  lens  scatter- 
ing, a traveling  fiber  optical  probe  was  assembled  along  the  optical 
streak  and  transverse  to  the  streak.  It  measures  the  wide-angle  scattered 
radiation  just  above  the  substrate.  This  automated  apparatus,  together 
with  the  scan  drive  and  recording  equipment,  is  depicted  in  Figure  105.  A 
typical  result  of  scanning  across  a Luneburg  lens  and  a section  of  wave- 
guide is  depicted  in  Figure  106.  The  Corning  7059  waveguide  loss  is 
0.4dB/cm  and  the  Ta20^  waveguide  lens  loss  is  1.74dB/cm. 

Another  example  for  Corning  7059  waveguide  on  thermally-oxidized 
silicon  is  depicted  in  Figure  107  for  wafer  #7  in  both  the  best  direction 
and  orthogonal.  The  angulationsof  this  data  are  due  to  fine  scattering 
fluctuations  while  the  line  is  a computer-generated  RMS  slope.  The  streak 
of  the  waveguide  in  Figure  107  is  barely  perceptible  to  the  eye. 


Figure  106  Waveguide  Scattering  Trace  Across  Ta?0r  Thln-Fllm  Luneburg  Lens  and  Subsequent 
Coming  7059  Waveguide. 


Another  kind  of  waveguide  scattering  observation  is  depicted  in 
Figure  108  showing  scattering  normal  to  the  waveguide  and  30°  above  the 
plane  of  the  waveguide.  Diffraction  coupling  has  been  used  at  the  input 
cleaved  edge.  Close  inspection  of  Section  (a)  shows  a number  of  longi- 
tudinal streaks  as  fine  structure.  They  have  been  traced  to  stepping 
of  the  cleaved  edge.  Two  kinds  of  scatters  are  apparent  in  Section  (a). 
The  "snowball"  effect  for  particles  much  larger  than  the  optical  wave- 
length is  pronounced  in  both  (a)  and  (b),  while  the  relatively  smooth 
fine-grained  distribution  is  substantially  reduced  as  the  angles  change. 
Another  example  of  a Corning  7059  waveguide  on  a thermally-oxidized 
silicon  substrate  is  illustrated  in  Figure  109.  Unfortunately,  the  cor- 
responding photograph  is  not  available.  The  overall  waveguide  assembly 
effective  measured  loss  is  2.8dB/cm.  Three  dashed  curves  are  inserted 
for  various  segments  which  have  a slope  representative  of  0.9dB/cm.  The 
effect  illustrated  in  Figurel09  has  been  traced  to  the  cleaning  process 
of  the  oxidized  surface  prior  to  the  deposition  of  the  waveguide.  A 
smear  appears  corresponding  to  each  of  the  dramatic  steps  in  the  scatter- 
ing profile.  Scattering  at  these  points  perturbes  the  confinement  in 
the  7059  so  that  radiation  is  absorbed  in  the  substrate.  Thus,  the 
actual  streak  traverses  the  guide  at  approximately  0.9dB/cm  until  it 
intersects  the  residue  left  from  cleaning  where  a significant  portion 
is  scattered  in  the  substrate  and  then  continued  at  the  0.9dB/cm  rate 
until  it  intercepts  the  next  scattering  section. 

Another  type  of  waveguide  structure  investigated  employed  a 
Corning  7059  waveguide  on  a phosphorous-doped  silicate  glass  and 
thermally-grown  Sit^  on  silicon  substrate.  The  use  of  both  isolation 
layers  permits  etching  steps  into  the  Si 02  as  may  be  required  for  a 
waveguide  transition  into  a photodetector  array.  The  addition  of  phos- 
phorous reduces  the  melting  temperature.  The  surface  is  heated  until 
surface  tension  from  the  partial  liquid  state  forms  relatively  smooth 
transitions  over  the  steps.  The  results  of  this  structure  are  ill  us 
trated  in  the  photographs  of  FigurellO.  Section  (a)  is  taken  at  normal 
Incidence  and  shows  the  typical  streak  and  one  scattering  center  at  the 
first  transition  edge.  The  optical  beam  traverses  three  additional  steps 
which  are  not  apparent  but  each  of  the  steps  are  reduced  in  height  in 


Diffraction  Coupling  to  Cleaved  Edge 


Tilted  Substrate 
to  Show  Output 
Distribution 


(b) 


Figure  108  Photograph  of  Optical  Streak  through  Corning  7059  Waveguide  on 
Thermally-Oxidized  Substrate  at  Normal  Incidents  (a)  and  at  30° 
above  Plane  Looking  at  the  Output  (b).  Input  Coupling  Provided 
Focusing  Radiation  on  Cleaved  Edge. 


142 


Intensity  (dB) 


Figure  109  Measured  Optical  Waveguide  Scattering  for  Corning  7059 
Waveguide  on  a Thermally-Oxidized  Silicon  Substrate. 

Two  Sharp  Steps  are  due  to  Scattering  Loss  and  Absorption 
in  Silicon  Substrate.  Waveguide  Assembly  Effective  Loss 
2.8dB/cm.  Reduced  Slope  Lines  Represent  0.9dB/cm  Loss. 

descending  order  moving  to  the  substrate  edge.  Section  (b)  of  the  photo- 
graph was  taken  at  45°.  It  shows  a first-step  transition  to  produce 
significant  scattering  in  the  forward  direction  and,  thus,  indicates 
that  the  step  height  is  too  much.  The  remainder  of  the  trace  to  the 
end  of  the  substrate  is  now  apparent  at  the  forward  scatter  angle.  The 
use  of  reflowed  phosphorous-doped  silicate  glass  provides  an  important 
design  freedom  for  integrated  optical  circuits  and  leads  to  a reduction 
of  scattering.  It  is  presumably  due  to  the  fact  that  its  process  reaches 
the  liquid  state  so  that  the  surface  tension  produces  a smoothing  effect. 

The  properties  of  silicon-oxynitride  were  investigated  as  an  al- 
ternative to  Corning  7059  and  Ta20g.  The  refractive  index  of  this  material 
can  be  changed  from  approximately  1 .47  - 2.1  depending  upon  the  SiOg  and 
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Figure  no  Photograph  of  Optical  Waveguide  Streak  Normal  to  Substrate. 
7059  Waveguide  on  Reflowed  Phosporous  Silica  Glass  on 
Thermally-Grown  Si O2  with  Stepped  Edges  on  Silicon  Substrate 
(a), and  Photograph  of  Forward  Scattered  Radiation  at  45°  (b). 
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Si 3N4  proportions.  The  first  result  is  shown  in  Figurelll,  together  with 
the  measured  scattering  intensity  profile.  Careful  inspection  shows 


circular  streaks  which  have  been  traced  to  a typical  semiconductor  clean- 
ing procedure,  spinning  the  wafer  with  solvents  and  a swab.  The  measured 

i 

loss  is  5.2dB.  The  film  thickness  is  approximately  3000A.  A i.«uch  better 
result  is  shown  in  Figurell2.  The  corresponding  measured  scattering  pro- 
file has  a slight  slope  indicating  a loss  of  0.5dB/cm.  An  expanded  in- 
tensity scale  has  been  included  to  show  the  minute  fluctuations.  The 
photograph  of  Figure  112  shows  a throughput  refraction  from  the  cleaved 
edge  and  the  transverse  spatial  distribution  surrounding  the  streak  con- 
tained within  the  waveguide.  This  gives  an  encouraging  comparison  with 
Corning  7059.  If  the  process  can  be  adapted  for  the  formation  of  lens 
profiles  and  for  the  formation  of  waveguide,  it  will  be  possible  to  form 
the  optical  circuit  using  a single  composition.  In  the  region  of  the 
lens,  a much  larger  percentage  of  Si 3N4  would  be  used.  The  SigN^  thin- 
film  process  employed  for  the  waveguide  in  Figurell2  has  utilized  a 
chemical  vapor  deposition  technique  at  high  temperature  and  at  high 
pressure.  This  process  was  used  to  obtain  a higher  concentration  of  SigN^ 
and  to  obtain  a larger  thickness  in  a reasonable  time.  A corresponding 
low-pressure  chemical  vapor  deposition  process  reported  in  the  literature 
also  yields  good  results  but  requires  a considerably  longer  time  to  obtain 
the  equivalent  thickness. 

Low-loss,  low-scattering  Corning  7059  waveguide  is  obtainable. 

Figure  113  illustrates  a prism  coupler  with  an  optical  beam  propagating 
through  the  waveguide  which  exits  from  the  cleaved  face.  The  waveguide 
streak  is  not  apparent  to  the  eye  or  discernable  by  the  available  in- 
strumentation. The  wide  collimated  beam  is  refracted  into  the  back- 
ground in  Figurell3  and  is  perfected  with  reflections  photographed 
in  Section  (b)  to  show  a well-defined  diffraction  pattern.  A more 
typical  low-loss  propagation  streak  is  shown  in  Figure  114  with  its 
throughput.  This  streak  has  fine  structure  due  to  a less-than-perfect 
prism  coupler  edge.  A scattering  line  defect  is  apparent  at  the  output 
due  to  prior  coupling  of  the  prism.  A measured  transverse  beam  dis- 
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Figure  112 


Photograph  of  Optical  Waveguide  (Good)  Streak  for  an 
Si 3N4 : Si O2  Waveguide  on  a Thermally-Oxidized  Silicon  Substrate 
(Increased  Exposure)  (a)  Together  with  Two  Corresponding 
Measured  Scattered  Intensity  Profiles  by  Travelling  Fiber 
Probe  Technique  (b)  of  Specimen  in  (a)  and  Expanded  Scale 
of  Profile  (c). 
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Figure  113  A Good  7059  Waveguide  Without  Observable  Streak  With  a 
Portion  of  the  Beam  Refracted  to  the  Background  (a)  and 
the  Beam  Projected  on  a Reflected  Surface  Showing  its 
Spatial  Structure  (b). 


Figure  114  Another  Propagation  Streak  for  7059  Waveguide  Showing 
Structure  Partially  Due  to  a Less-Than-Perfect  Prism 
Edge  With  its  Throughput. 


Figure  115  Transverse  Beam  Distribution  of  6328^  Radiation  Propagating 
Through  10  cm  of  7059  Waveguide  Corresponding  to  the  Photo- 
graph of  Figure  113. 


tributlon  of  a low-scattering,  low-loss  Corning  7059  waveguide  is  de- 
picted in  Figure  115  . This  propagation  loss  is  less  than  O.ldB/cm. 
Loss  measurements  at  this  low  scattering  level  requires  further  tech- 
nique refinement. 


Figure  n6  Transverse  Cross  Section  of  V-Type  Absorptive  Spatial 


Filter. 

E.  SPATIAL  FILTER 

Because  the  optical  Fourier  transform  requires  a high  degree  of 
mode  purity  to  obtain  diffraction-limited  resolution  and  because  the 
laser  source  may  not  emit  a single  pure  transverse  mode  of  the  lowest 
order  and  because  some  scattering  may  be  present,  an  investigation  of 
a spatial  filter  has  been  included.  Typical  spatial  filters  consist 
of  small  aperture  stops  located  at  the  focus  of  inverting  beam  expanders. 
They  are  largely  reflective  in  nature  with  abrupt  boundaries.  The 
spatial  filter  considered  herein  for  inverting  beam-expansion  lens  sys- 
tems employs  an  absorptive  iris  at  the  focus  using  a planar  waveguide 
structure  supported  by  a thermally-grown  oxide  isolation  layer  on  a 
silicon  substrate.  The  absorptive  iris  spatial  filter  is  formed  by 
removal  of  the  thermally-grown  oxide  so  that  the  planar  waveguide  is 
stepped  down  and  is  adjacent  to  the  absorbed  silicon  substrate  in  the 
iris  section.  A transverse  cross  section  through  the  V iris  apex  is 
shown  in  Figure  116  with  a plan  view  of  the  stepped-down  absorptive  sec- 
tion shown  in  Section  (a)  of  Figurell7.  The  V gap  is  5 microns.  Fiducial 
marks  are  used  for  accurate  positioning  in  the  focus  of  a Luneburg  lens 
beam  expansion  system.  The  operation  of  this  spatial  filter  is  depicted 
in  Section  (b)  of  Figure  117  using  a 6328&  radiation.  The  entrance  beam 
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Figure  118  Cross  Sections  of  Near  Field  Intensity  in  Output  of  V-Type 
Absorptive  Spatial  Filter. 

is  apparent  together  with  the  exit  beam  and  a few  scatterers  (out  of 

focus).  The  input  beam  is  intentionally  offset  and  exhibits  a typical 

streak.  Transverse  profiles  of  another  spatial  filter  through  three 

sections  of  the  beam  at  different  displacements  from  the  spatial  filter 

are  depicted  in  Figurell8  where  the  spatial  filter  aperture  is  7 microns. 

The  continuous-wave, power-handling  capacity  of  Corning  7059  waveguide 

6 2 

on  thermally-grown  Si^  at  the  focus  is  greater  than  0.4  x 10  watts/cm  . 

The  V-groove  is  appropriate  for  tightly  focused  beams,  however,  to 
introduce  a more  sharply  defined  beam  truncation  for  a longer  focal  length 
than  the  system,  the  V shape  has  been  altered  to  more  closely  approximate 
the  hyperbolic  beam  contour  near  the  focus.  Figure  119  illustrates  such  a 
spatial  filter  plan  view  with  a magnified  view  of  the  constricted  wave- 
guide section. 
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Figure  119  Macrophotograph  of  Absorptive  Spatial  Filter  With  Contoured 
Gap  Designed  for  Use  With  F/4  Waveguide  Lenses  (a)  and 
Detail  Showing  the  Contoured  Gap  of  5.6  Microns  (b). 


F.  CHANNEL  WAVEGUIDE  FAN  OUT 

For  an  integrated  optical  Fourier  transform  circuit,  an  optical  wave- 
guide fan  out  array  was  investigated  to  accommodate  high  density  spatial 
informatipn  in  the  lens  focal  plane  to  a less  dense  photodetector  linear 
array. 


Figure  120  Transverse  Cross  Section  of  Channel  Waveguide  Fan-Out. 

A transverse  cross-section  of  a channel  waveguide  fan  out  array  is 
depicted  in  Figure  120.  A plan  view  of  the  channel  waveguide  fan  out 
appears  in  the  photographs  of  Figure  121  (a)  and  its  enlargement  in  section 
(b).  Corning  7059  waveguide  has  been  employed  on  a comparatively  thick 
thermally-grown  isolation  layer  on  a silicon  substrate.  Channel 
separation  has  been  provided  by  photolithography  using  a chemical  etch 
to  define  channels.  The  thermally- grown  SiO^  isolation  layer  thickness 
is  typically  4 microns,  while  the  etch  depth  defining  the  channels  is 
approximately  1.5  microns  deep.  The  Corning  7059  waveguide  approximately 

O O 

5000  A thick  will  fill  the  separation  channels  to  approximately  2000  A 
due  to  the  shadowing  of  the  elevated  waveguide.  The  operation  of  two  channel 
waveguide  fan  out  specimens  is  illustrated  in  Figure  122  for  radiation  where 

o 

6328  A radiation  has  been  coupled  into  and  traverses  a single  channel 
waveguide.  Some  residual  illumination  discloses  the  presence  of  the  other 
channels.  Strong  Mie  scattering  from  the  waveguide  edges  is  apparent. 

The  channel  waveguide  losses  are  unfortunately  comparatively  high  (less 
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Figure  121  Plan  View  of  Channel  Waveguide  Fan-Out  (a)  and  Enlargement 
(b). 

than  10  dB/cm),  however,  this  loss  should  not  impair  the  photodetector 
array  because  the  waveguide  fan  out  length  probably  will  not  exceed  2 milli- 
meters. Coupling  of  radiation  into  a single  waveguide  channel  at  the  end 
with  a larger  period  shows  that  there  is  no  significant  coupling  to  ad- 
jacent channels  as  these  waveguides  converge. 
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III.  REVIEW 


A.  CONCLUSIONS 

Integrated  optical  waveguide  positive  inhomogeneous  refractors 
identified  as  the  generalized  Luneburg  lenses  have  been  investigated, 
calculating  the  performance  and  demonstrating  their  operation,  yielding 
diffraction-limited  results.  A method  to  calculate  and  predict  the  per- 
formance and  a method  to  synthesize  the  mask  required  to  realize  thin- 
film  Luneburg  lenses  has  been  developed  showing  that  the  generalized 
Luneburg  lens  for  integrated  optics  is  on  a sound,  well -developed  foun- 
dation. Further,  the  lenses  exhibit  the  requisite  properties  for  an 
optical  Fourier  transform  "chip"  that  may  be  employed  as  a real-time 
microwave  spectrum  analyzer. 

The  approach  has  emphasized  the  application  of  Ta205  as  the  high- 
index  lens  film  used,  together  with  Corning  7059  as  a waveguide,  all  of 
which  are  deposited  upon  thermally-grown  SiCL,  as  an  isolation  layer  on 
a silicon  substrate.  This  structure  permits  the  integration  of  the 
optical  components  with  their  electronic  counterpart  on  a single  sub- 
strate. This  material -structure  and  its  associated  processing  is 
compatible  with  the  processing  of  electronic  circuit  elements  using 
processes  familiar  to  the  semiconductor  industry. 

The  thin-film  Luneburg  lenses  are  formed  by  RF  sputter  deposition 
of  the  lens  material  through  a circular  mask  with  a prescribed  edge  shape 
on  the  waveguide.  The  procedure  to  synthesize  the  mask  shape  is  included. 
Accurate  synthesis  of  the  "true"  Luneburg  lens  profile  yields  diffraction- 
limited  results  verified  experimentally. 

A numerical  procedure  has  been  developed  and  applied  to  compute 
the  performance  of  waveguide  generalized  Luneburg  lenses.  Starting  with 
the  waveguide  lens  thickness  profile  determining  the  ray  traces  in  the 
waveguide  leads  to  the  ray  intercepts  in  the  image  plane.  With  this 
data,  the  wavefront  is  determined  and  the  resulting  intensity  diffraction 
pattern  is  calculated.  The  procedure  allows  the  introduction  of  various 
types  of  waveguide  thickness  profile  errors  to  be  introduced  and  the 
calculation  of  the  effect  upon  the  diffraction  pattern.  Numerous  ex- 
amples have  been  included  showing  the  effects  for  F/l,  F/2.5  and  F/4.5 
lenses. 
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Waveguide  Luneburg  lenses  with  F/4  to  F/25  have  been  realized 
with  diffraction- limited  beamwidths.  The  smallest  beamwidth  observed 
is  1.65  microns  confirming  the  design  value  and  the  procedure.  Decay- 
ing side  lobes  which  can  be  held  below  17dB  for  low-scattering  lenses 
and  waveguide  can  be  achieved.  The  side  lobes  decay  into  the  optical 
scattering  level  and/or  instrumentation  noise  typically  on  the  order  of 
30dB.  These  results  have  confirmed  mask  synthesis  procedure.  It  has 
been  found,  experimentally,  that  a mask  design  for  F/4  is  applicable 
to  longer  focal  length  lens  yielding  diffraction- limited  results. 

The  effects  of  optical  waveguide  scattering  have  been  reduced 
by  improved  processing  procedures  so  that  consistent  reproducibility 
has  been  achieved  <0.3dB/cm  while  some  results  have  broken  the  barrier 
below  O.ldB/cm.  Corning  7059,  4 Si02:Ba0,  Si02(P205),  Si304:Si02  based 
waveguides  were  considered  with  Ta205  and  Nb205  lenses.  Consideration 
of  the  effect  of  waveguide  and  lens  scattering  upon  the  dynamic  range 
and  resolution  of  an  integrated  optical  Fourier  transform  is  treated 
showing  that  the  best  measured  scattering  is  consistent  with  the  allow- 
able lens  and  waveguide  scattering  for  this  application. 

Some  further  effort  reducing  the  scattering  by  understanding 
the  mechanism  and  introduction  of  process  refinements  is  desirable. 

This  is  particularly  true  for  the  Ta205  and  Nb205  lenses  when  the 
thickness  is  on  the  order  of  hundreds  of  Angstroms. 

Other  integrated  optical  circuit  elements  for  the  Fourier  trans- 
form signal  processing  function  include  a spatial  filter,  reflectors  and 
channel  waveguide  fan-out  structures.  Spatial  filters  using  an  absorptive 
V-type  iris  were  demonstrated  for  use  with  inverting  beam  expansion  lens 
systems.  Several  beam  expansion  lens  systems  have  been  demonstrated 
showing  that  the  mask  positioning  can  yield  a collimated  beam  reproducibly. 
Cleavage  through  the  substrate  and  overlaying  films  can  produce  an  inter- 
face of  sufficient  perfection  to  implement  total  internal  reflectors. 
Positioning  of  the  cleavage  to  photolithographic  patterns  was  demonstrated. 
For  some  integrated  optical  circuits,  planar  waveguide  coupling  to  a 
discrete  linear  array  of  photodetectors  is  desirable.  For  this  purpose, 
a high-linear  density  channel  waveguide  fan-out  structure  has  been  demon- 
strated. 
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The  results  of  this  effort  and  its  predecessors  have  established 
a firm  technology  base  for  waveguide  Luneburg  lenses  on  a thermally- 
oxidized  silicon  substrate  providing  an  attractive  approach  to  implement- 
ing an  integrated  optical  Fourier  transform  "chip"  that  may  be  used  for 
spectral  analysis  and  other  signal  processing  functions  where  real-time 
bandwidth  operation  is  essential.  The  approach  uses  well-known  established 
manufacturing  processes  with  some  refinement  and,  thus,  the  "chip"  can  be 
implemented  economically  as  quasi -monolithic  optical  circuit  integrated 
with  its  associated  electronics. 
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SUMMARY 

Thin-film  waveguide  Luneburg  and  geodesic  lenses  are  baing  developed 
for  near-infrared  integrated  optics.  Considerations  which  yield  a prefer- 
ence for  these  structures  are  enumerated,  together  with  specifics,  for 
generalized  Luneburg  lenses  employing  Ta^Og  (lens)  with  Corning  7059  (wave- 
guide) on  thermally-grown  SiO^  on  Si  substrates.  Spheric  geodesic  lenses 
formed  with  the  same  waveguide,  isolation  layer,  and  substrate  materials 
are  included. 

INTRODUCTION 

Optical  waveguide  lenses,  reflectors,  and  lens  systems  [1],  together 
with  surface  elasto-optic  spatial  modulators  [2]  and  waveguide  coupled  de- 
tectors [3],  are  being  developed  in  an  integrated  optics  format  to  perform 
the  Fourier  transform  signal  processing  function. [4,5]  Coherent  optical 
signal  processing  requires  near  diffraction-limited  performance,  and  the  in- 
tegration of  these  elements  on  a single  substrate  requires  compatible  mater- 
ials-processing  technology.  For  compactness  and  wavelength  stability,  a 
distributed-feedback  electron  injection  laser  is  preferred.  Economics, 
alignment  stability,  and  small  size  compel  the  use  of  a single  substrate 
which  requires  lenses  of  comparatively  small  f-number  and  reflectors  for 
folding  of  the  optical  system  so  as  to  be  contained  within  the  available  sub- 
strate area. 

Two  waveguide-material  structures  are  receiving  attention. [6]  They 
are  (1)  in-  or  out-diffused  LiNbOg  waveguide  [7]  because  of  its  low-loss  and 
because  of  large  piezoelectric  coupling  which  leads  to  efficient  wideband 
elasto-optic  interactions  and  (2)  Corning  7059  glass  waveguide  on 
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thermally-grown  Si02  on  silicon  because  of  minimum  cost  considerations 
where  sputtered  ZnO  is  employed  as  the  piezoelectric  to  launch  Rayleigh 
waves. [8-10] 

WAVEGUIDE  LENS  CONSIDERATIONS 

The  relatively  large  refractive  index  of  LiNbO^  and  the  lack  of 
thin-film  materials  with  substantially  larger  index  forces  the  appli- 
cation of  geodesic  lenses  with  LiNbOg.  Geodesics  employ  curved  surface 
topology  which  alters  the  propagation  length  because  rays  propagate 
according  to  Fermat's  principle. [11-13]  These  surfaces  do  not  conform 
to  the  normal  concept  of  an  integrated  format;  however,  simple  spherical 
depressions  iri  the  substrate  are  easily  fabricated  as  exemplified  by  bi- 
focal glasses.  Fortunately,  metal  in-  or  out-diffusion  on  LiNb03  leads 
to  low-loss  waveguide  with  minimum  scattering,  the  latter  of  which  may 
compromise  the  transform's  resolution  and/or  dynamic  range.  Experimental 
spheric  depression  geodesic  lenses  are  illustrated  in  Figure  1,  together 
with  use  of  two  lenses  as  a beam  expander.  Corning  7059  waveguide  is  em- 
ployed on  the  thermally-grown  oxide  as  the  isolation  layer  on  the  silicon 
substrate.  Three  rays  traversing  the  beam  expander  are  shown  with  scatter- 
ing at  the  sharp  geodesic  edges.  Radiation  (6328fl)  traverses  the  sharp 
geodesic  edge  because  of  the  inherent  radius  due  to  the  thermal  oxide 
growth.  Mode  conversion  and  scattering  at  the  sharp  planar-curved  sur- 
face intersection  can  be  minimized  by  a toroidal  lip. [14] 

Spherical  geodesic  lenses  exhibit  spherical  and  curvature  of  field 
aberrations,  the  former  of  which  can  be  corrected  and  the  latter  is  of  no 
consequence  to  the  Fourier  transform.  Spiller  and  Harper  [15]  have  shown 
that  use  of  a thin-film  additional  layer  coincident  with  the  geodesic  can 
suppress  the  spherical  aberration.  Southwell  [16]  has  shown  that  a toroidal 
lip  may  be  used  to  reduce  the  spherical  aberration.  Computed  ray  traces 
through  a spherical  geodesic,  with  and  without  a toroidal  lip,  are  shown  in 
Figure  2 using  Southwell's  procedure.  Vahey  and  Wood  [17]  have  shown  that 
the  oblate  spheroid  can  also  be  used  to  correct  the  spherical  aberration. 

Perfect  geometrical  focusing  can  be  achieved  with  the  Rinehart[18]-Luneburg 
aspheric  geodesic  surface  which  focuses  on  the  edge.  For  an  arbitrary  fo- 
cal length,  Morgan  [19]  has  generalized  the  Luneburg  class  [20]  of  inhomo- 
geneous lenses, and  Southwell  [16]  has  computed  its  geodesic  analog.  In 
general,  the  aspheric  configurations  are  difficult  to  generate,  polish 
sufficiently,  and  test  when  compared  with  spherics. 
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Spheric  geodesic  lenses  have  been  fabricated  using  Ti-indiffused 
waveguide  on  LiNbOg,  Ag+ion-exchange  in  glass,  thermally-oxidized  Si,  and 
anodized  GaAs  substrates.,  Measured  diffraction  patterns  clearly  show  evi- 
dence of  spherical  aberration  and  scattering  due  to  less-than-perfect  polish 
in  and  at  the  edge  of  the  geodesic. 

Thermally-grown  SiO^  on  Si  has  yielded  the  smoothest  low-scattering, 
low-loss  substrate  surface  for  the  deposition  of  low-loss  Corning  7059 
waveguide  with  refractive  elements.  The  loss  of  this  waveguide  can  be  in 
the  0.1-0. 5 dB/cm  range.  The  low  cost  of  Si  wafers  with  sufficient  dimen- 
sions and  the  current  mature  processing  technology  of  Si  together  with  the 
fact  that  the  optical  circuit  can  be  integrated  with  its  electronic  counter- 
part provides  the  stimulation  for  its  exploitation.  Tantalum  pentoxide  and 
f^O,.  have  been  employed  to  form  homogeneous  and  inhomogeneous  thin-film 
lenses  and  lens  systems.  Their  losses  are  in  the  0.5-3  dB/cm  range  depend- 
ing upon  the  deposition  process.  The  larger  value  is  usually  obtained  for 
very  thin  layers  over  Corning  7059,  and  the  lower  value  may  be  obtained  when 
immersed. 

The  cleavage  of  single-crystal  Si  can  extend  through  the  isolation 
and  waveguide  layers  to  form  perfect  reflectors,  so  as  to  permit  folding  of 
the  optical  system.  A protective  layer  superstrate  remains  to  be  identi- 
fied which  will  be  probably  an  organic  deposited  from  solution. 

Because  the  waveguide  effective  refractive  index  is  related  to  the 
film  thickness,  it  is  possible  to  create  both  homogeneous  and  inhomogeneous 
refractors  using  homogeneous  materials.  The  maximum/minimum  ratios  of  ef- 
fective refractive  indices  are  typically  in  the  range  of  1.4-0. 7.  The  latter 
arises  when  the  ambient  waveguide  is  dense  and  the  lens  element  thickness  is 
reduced.  Thus,  it  is  possible  to  create  positive  or  negative  waveguide  lens 
elements  using  either  convex  or  concave  interfaces.  We  avoid  the  homogeneous 
thin-film  refractors  because  it  is  difficult  to  produce  optically  smooth 
interfaces  lithographically. 

The  circular-symmetric  inhomogeneous  refractor  due  to  Luneburg  [20] 
has  been  demonstrated  in  optical  thin  films. [1,21]  Interest  in  the  Luneburg 
lens  arises  because  it  is  free  of  all  aberrations  except  curvature  of  field 
and  because  of  its  simplicity  of  fabrication  as  a thin-film  waveguide  struc- 
ture. A Luneburg  lens  will  form  perfect  geometric  images  of  two  given 
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concentric  circles  on  each  other, [19]  and,  if  one  is  at  infinity,  the  lens 
will  focus  a collimated  beam  from  any  direction  on  the  opposite  edge.  This 
circular  syumetry  is  of  particular  interest  to  integrated  optics  where  the 
substrate  area  is  limited  because  voiding  of  the  optical  axis  may  use  the 
same  lens  for  several  purposes. 

Southwell  [22,23]  has  developed  a numerical  technique  for  ray  trac- 
ing through  inhomogeneous,  circular-symmetric,  thin-film  waveguide,  general- 
ized Luneburg  lenses.  Two  examples  are  illustrated  in  Figure  3 for  the 
“classic"  [20]  and  the  generalized  Luneburg  lens  for  Morgan's  parameter 
S=2  [19].  Southwell's  numerical  procedure  determines  the  ray  intercept 
and  uses  Hamilton's  theory  to  extract  the  pupil  function  phase  error  fol- 
lowed by  the  Fourier  transform  to  determine  the  intensity  diffraction  pat- 
tern. The  technique  yields  diffraction-limited  results  in  the  Hewlett 
Packard  9830A  calculator.  Computations  and  experiments  show  that  errors 
in  deposition  time  (waveguide  thickness)  lead  to  the  defocused  aberration 
which  may  be  corrected.  Errors  of  the  waveguide  lens  radial  refractive 
index  profile  may,  thus,  be  traced  to  the  resulting  diffraction  pattern 
and  to  the  mask  which  shapes  the  lens  deposition  profile. 

Generalized  Luneburg  lenses  [19]  have  been  fabricated  with  f-numbers 
of  0.5  (classical)  to  greater  than  10  using  a Ta205  overlayer,  underlayer, 
or  interleaved  within  the  Corning  7059  waveguide.  The  interleaved  sym- 
metrical structure  is  preferred  to  suppress  waveguide  mode  conversion  as 
opposed  to  the  overlay  or  underlay  structure.  Luneburg  lenses  with  f- 
numbers  greater  than  6 have  been  realized  using  "all"  7059  thin  films. 

These  lenses  have  been  formed  by  RF  sputtering  [24]  from  a wide-angle 
source  through  a conical  mask  [21].  In  general,  when  the  cone  angle  is 
reduced  and  the  maximum  thickness  is  reduced,  an  increase  of  focal  length 
is  obtained.  The  mask,  to  realize  the  required  inhomogeneous  lens  profile, 
can  be  determined  effectively  from  molecular  geometric  optics  for  the  depo- 
sition process  using  the  mask  to  control  shadowing. 

The  operation  of  two  thin-film  overlay  Ta205  Luneburg  lenses  as  a 
beam  expander  (X3)  is  shown  in  Figure  4 where  the  full  aperture  of  the  ex- 
pander lens  is  uniformly  illuminated  in  one  photograph, and  a single  ray 
is  employed  in  thp  other  to  discern  mode  conversion  in  the  expander  over- 
lay lens.  Unfortunately,  the  effect  of  substrate  damage  extending  into  the 
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thermally-grown  Si02  is  apparent  near  the  focus  of  the  expansion  system. 

Continuous  power-handling  capacity  of  Corning  7059  waveguide  on  thermally- 

6 ? 

grown  Si 02  at  the  focus  is  greater  than  0.4  x 10  watts/cm  . Damage  due 
to  attachment  and  detachment  of  coupling  prisms  prior  to  this  photography 
is  evident  in  both  the  input  and  output.  The  throughput  exiting  the  Si 
wafer  edge  scattered  from  black  paper  is  also  evident.  Control  of  lens 
position  and  thickness  yields  a collimated  beam.  The  single  6328A  ray 
trace  in  Figure  4 enters  the  expansion  lens  as  the  TEq  mode  where  it  is 
partially  converted  to  the  TE-j  mode.  After  this  mode  traverses  the  ex- 
pansion lens,  it  is  reconverted  to  the  TEq  mode,  so  that  two  separate 
rays  continue  through  the  remainder  of  the  expansion  system  as  TEq  modes. 

A typical  measured  diffraction  pattern  of  a generalized  Luneburg 
lens,  like  the  objective  in  Figure  4,  is  illustrated  in  Figure  5.  Its 
beamwidth,  3 dB,  is  5y.  Uniform  illumination  was  employed  and,  thus, 
the  increase  of  side  lobe  level  and  null  fill-in  is  indicative  of  an 
error  in  the  lens  thickness  profile,  producing  a tt/2  quadratic  phase 
error.  The  asymmetric  wide  angle  side  lobe  can  be  traced  to  a Mie 
scatter  in  the  lens. 

CONCLUSIONS 

A thin-film  waveguide  Luneburg  lens  is  considered  the  most  viable 
planar  approach  based  upon  the  use  of  oxides;  and  when  high  refractive- 
index  substrate  materials  must  be  employed,  the  geodesic  becomes  a feas- 
ible alternative.  Accurate  realization  of  the  required  Luneburg  lens 
profile  and  correction  of  geodesic  aberrations  is  continuing  to  receive 
attention.  Further  effort  to  reduce  Rayleigh  and  Mie  scattering  in  the 
waveguide  is  required  because  it  limits  the  resolution  and  dynamic  range 
of  the  Fourier  transform. 
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Figure  2.  Calculated  Ray  Traces  Through  Spheric  Geodesic  Lens  Having 
(a)  Sharp  Edge  and  (b)  Toroidal  Lip  for  Parameters  Listed. 
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Figure  3.  Calculated  Ray  Traces  Through  Luneburg  Lenses  for  (a) 

"Classic"  S=l,  and  (b)  Generalized  S=2  Derived  From  Wave- 
guide Thickness  Profiles  Yielding  Diffraction-Limited  Results. 


Luneburg  Expansion  Lens 


Figure  4.  Luneburg  Lens  Beam  Expander  (X3)  Employing  Ta-O,-  Lenses  With 
Corning  7059  Waveguide  on  Thermally-Oxidized  si3(a)  Operation 
of  Beam  Expander  Showing  Lens  Scatter,  Defects  and  Throughput, 
and  (b)  a Single  Ray  Trace  Through  Same  Beam  Expander  Showing 
Mode  Conversion  in  Expansion  Lens  but  With  Both  Rays  Traversing 
the  Remainder  of  the  Waveguide  as  TEQ  Modes. 
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MASK  SYNTHESIS  FOR  DIFFRACTION-LIMITED  WAVEGUIDE  LUNEBURG  LENSES* 
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An  integrated  optical  circuit  for  RF  signal  spectral  analysis 
requires  a near  diffraction-limited  lens  to  perform  the  Fourier  trans- 
form.' A thin-film  waveguide  Luneburg  lens  using  Ta20g  in  the  lens 
region  with  Corning  7059  waveguide  on  thermally-oxidized  silicon  sub- 

4 

strate  is  a low-cost  approach  providing  for  integration  with  a surface 
elasto-optic  spatial  modulator  and  a waveguide  coupled  photodetector 
array  in  the  substrate.  Waveguide  Luneburg  lenses  are  formed  by  sputter- 
ing high-index  films  through  a circular  mask  with  a shaped  edge. 

The  generalized  Luneburg  lens  is  free  of  all  aberrations  except 
for  curvature  of  field.  It  is  a circular-symmetric,  inhomogeneous  lens 
requiring  a specific  radial  refractive  index  profile  and,  thus,  a wave- 
guide thickness  profile.  Southwell  has  developed  a convenient  technique 
for  gradient-index  ray-tracing  and  has  applied  it  to  generalized  Luneburg 
lenses  determining  ray  intercepts,  the  corresponding  phase  wavefront 
and  the  intensity  diffraction  pattern.  This  ray-tracing  technique  is 
sufficiently  accurate  to  yield  diffraction-limited  results  for  the  "true" 
radial  profile  when  executed  with  an  HP  9830  calculator.  The  practical 
problem  addressed  herein  is  the  synthesis  of  the  mask  edge  shape  such  as 
to  yield  film  thickness  profiles  that  operate  as  Luneburg  lenses  with 
near  diffraction-limited  performance. 
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A computational  procedure  has  been  developed  to  derive  the  film 
thickness  profile  sputtered  through  masks  with  shaped  edge.  It  is 
assumed  that  the  molecules  behave  like  optical  rays,  that  their  source 
distribution  from  the  target  is  Lambertian,  and  that  the  mean-free  path 
length  is  greater  than  the  mask  dimensions;  and,  thus,  the  problem  re- 
duces to  shadowing  from  an  incoherent  extended  source.  The  circular 
mask  edge  is  simulated  by  an  interconnected  series  of  conical  sections. 
The  resulting  radial  thickness  profile  is  converted  into  the  effective 
waveguide  refractive  index  profile  for  the  designated  mode  kind  and 
order.  Southwell's  ray-tracing  program  is  then  employed  to  evaluate 
the  diffraction  pattern  in  the  focal  line.  Thus,  knowing  the  effect 
of  deviations  from  the  "true"  profile  and  recognizing  the  mask  shadow- 
ing effects,  it  is  possible  to  approximate  the  required  "true"  profile 
as  demonstrated  by  the  examples  to  be  cited. 

These  procedures  have  been  employed  to  derive  the  mask  shape 
edge  for  thin-film  waveguide  Luneburg  lenses  having  f-number  of  4, 
which  experimentally  yields  near  diffraction-limited  performance.  The 
mask  shape  edge  for  our  sputtering  system  geometry  and  the  resulting 
calculated  thickness  profile  is  depicted  in  Figure  1 including  an  in- 
sert showing  the  measured  edge  profile.  Using  this  radial  thickness 
profile  as  an  input  to  Southwell's  ray-tracing  program,  the  resulting 
calculated  ray  traces,  wavefront  phase,  and  intensity  diffraction 
pattern  are  shown  in  Figure  2.  Note  that  the  aperture  has  been  limited, 
and  that  the  simulation  has  allowed  a wavefront  phase  error  which  does 
not  seriously  compromise  the  diffraction  pattern.  The  lens  employs  an 
8.2mm  diameter  overlay  with  thickness  of  3000ft  at  the  center  on  a Corn- 
ing 7059  waveguide  for  a thickness  of  4500ft. 

The  measured  diffraction  pattern  for  this  Luneburg  lens  thick- 
ness profile,  using  a 2mm  Gaussian  input  beam  at  6328ft,  is  depicted  in 
Figure  3 using  a pair  of  prism  couplers,  an  F/3.5  Tropel  diffraction- 
limited  relay  lens  providing  a magnification  of  2.5,  and  a traveling 
slit  detector  of  lu  width.  The  diffraction-limited  beamwidth  in  the 
waveguide  should  be  2.2m  at  the  1/e  level;  whereas,  the  measured  width 
is  3.3u.  It  should  be  noted  that  the  relay  lens  diffraction  pattern  is 
convolved  with  the  Luneburg  lens  diffraction  pattern  limiting  the  ob- 
served result  to  3p.  The  high  level  first  side  lobe  is  due  to  spherical 
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aberration  introduced  by  the  output  prism  coupler.  The  large  shoulder 


lobes  are  in  part  due  to  scattering  effects,  particularly  at  the  edge 

of  the  input  coupler. 
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H«dius  of  the  Luneburg  Lens 


Figure  1.  Computed  F/4  Luneburg  Lens  Thickness  Profile  with  Mask  Edge 
Shape.  Insert:  Measured  Lens  Edge  Profile. 


Figure  2.  Computed  Ray  Traces,  Wavefront  and  Intensity  Diffraction 
Pattern  Using  Mask  Profile  of  Figure  1. 


Figure  3.  Measured  Diffraction  Pattern  of  F/4  Luneburg  Lens  Derived  From 
Mask  Edge  Shape  Illustrated  in  Figure  1. 
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EFFECT  OF  OPTICAL  WAVEGUIDE  SCATTERING  ON  THE  INTEGRATED 
OPTICAL  SPECTRUM  ANALYZER  DYNAMIC  RANGE* 

J.  T.  Boyd  and  D.  B.  Anderson 

ABSTRACT 

The  dynamic  range  limited  by  optical  scattering  of  the  Integrated 
optical  spectrum  analyzer  has  been  considered.  Calculations  of  dynamic 
range  as  a function  of  waveguide  attenuation  and  surface  roughness 
characteristics  are  presented.  Dynamic  range  Is  found  to  decrease  with 
Increasing  attenuation,  to  be  nearly  Independent  of  the  nature  of  the 
surface  roughness,  and  to  be  Independent  of  the  Fourier  transform  lens 
focal  length.  Adjacent  channel  rejection  Is  found  to  fall-off  slowly. 

The  effect  of  lens  scattering  Is  shown  only  to  be  Important  for  cases  In 
which  waveguide  scattering  Is  very  small  and  lens  scattering  Is  very  high. 

I.  INTRODUCTION 

The  Integrated  optical  approach  to  spectrum  analysis  of  electronic 
signals  in  the  hundreds  of  MHz  frequency  range  has  been  described  In 
previous  papers. x»2»3  In  these  previous  discussions  the  Important  perfor- 
mance parameter  of  signal  dynamic  range  has  been  discussed  only  qualita- 
tively. This  dynamic  range  Is  expected  to  be  limited  by  optical  scattering 
or  sensor  array  dynamic  range.  We  consider  herein  the  limitation  due  to 
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optical  scattering.  Although  some  progress  has  been  realized  recently 
in  reducing  scattering  in  optical  waveguides,4  scattering  is  still  ex- 
pected to  be  more  of  a limiting  factor  than  detector  dynamic  range. 

Numerical  results  presented  in  this  paper  support  this  conclusion. 

To  evaluate  the  effect  of  scattering  on  dynamic  range  in  a configuration 
such  as  that  of  the  integrated  optical  spectrum  analyzer,  a knowledge  of 
scattering  behavior  in  the  plane  of  the  waveguide  surface  due  to  an 
optical  wave  having  a finite  width  in  this  plane  is  required.  Theoretical 
analyses  of  scattering  in  optical  thin-film  waveguides  are  generally 
two-dimensional,  in  that  it  is  assumed  that  the  optical  wave  and  waveguide 
have  uniform  properties  in  the  direction  perpendicular  to  the  propagation 
direction  in  the  plane  of  waveguide  surface.5*6  The  emphasis  is  on 
obtaining  a solution  which  is  sufficiently  simple  so  as  to  be  readily 
amenable  to  an  engineering  calculation,  rather  than  on  performing  a three 
dimensional  analysis  of  scattering  in  optical  waveguides.  It  is  hoped 
that  the  numerical  results  presented  herein  will  stimulate  experimental 
investigation  of  scattering  in  optical  waveguides  in  the  plane  of  the 
surface  similar  to  previous  experimental  two-dimensional  investigations. 7 


II.  ANALYSIS 


A chip  configuration  of  the  integrated  optical  spectrum  analyzer 
is  shown  in  Fig.  1.  The  lens  regions  are  circular  in  view  of  the  fact 
that  the  generalized  Luneburg  and  geodesic  types  currently  appear  most 
promising.8*9  The  size  and  positions  of  these  regions  are  consistent  with 
desired  values  of  time-bandwidth  products  and  frequency  resolution.2*3 
The  compact  structure  in  Fig.  1 includes  a semiconductor  laser  coupled  to 
a silicon  optical  waveguide  structure.  Coherent  laser  light  propagating 
in  a thin-film  waveguide  is  expanded  and  collimated  by  a waveguide  lens 
system  prior  to  being  incident  on  a surface  elastic  wave  which  has  been 
excited  by  the  incoming  signal.  Diffraction  of  light  by  the  surface 
elastic  wave  occurs  in  such  a way  that  light  is  diffracted  at  different 
angles  corresponding  to  different  signal  frequency  components.  For  a 
range  of  signal  frequencies  and  wave  intensities  the  diffracted  light 
amplitude  is  linearly  proportional  to  the  signal  surface  elastic  wave 
amplitude.  A waveguide  lens  is  then  positioned  so  as  to  perform  a 
one-dimensional  Fourier  transform.  Signal  spectral  information  which  is 
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elastic  wave  region  then  becomes  represented  by  the  transverse  spatial 
variation  of  light  in  the  Fourier  transform  plane.  The  intensity  of  the 
light  spatial  distribution  in  the  transform  plane  corresponds  to  the  signal 
power  density  spectrum  and  can  then  be  converted  to  an  electrical  signal 
by  means  of  an  integrated  slab  waveguide  - CCD  imaging  array  such  as  that 
discussed  previously. 10 

Optical  scattering  affecting  spectrum  analyzer  performance  occurs 
in  the  waveguide  surrounding  the  lens  region,  in  the  lens  region,  and  in 
the  transition  between  these  two  regions.  However,  scattering  in  the 
waveguide  region  adjacent  to  the  detector  is  expected  to  dominate  other 
contributions  in  that  the  acceptance  angle  of  a detector  array  element 
for  scattered  radiation  diminishes  as  the  distance  between  the  array  element 
and  the  scattering  center  increases.  This  expectation  is  verified  by 
numerical  results  presented  in  the  next  section. 

Dynamic  range  as  limited  by  waveguide  scattering  is  defined  as  the 
ratio  of  the  power  reaching  a single  sensor  element  for  a single  frequency 
electronic  signal  input  to  the  power  scattered  to  the  adjacent  sensor 
element.  It  is  assumed  that  the  focal  length,  aperture,  and  detector 
array  element  spacing  are  chosen  so  that  the  width  of  the  focused  optical 
beam  corresponding  to  a single  frequency  input  equals  the  array  element 
size.  Dynamic  range  is  thus  defined  as  a measure  of  the  number  of  resolvable 
signal  levels  for  a narrowband  input.  For  wideband  signals  more  scattering 
occurs  and  thus  fewer  resolvable  levels  will  exist.  However,  dynamic  range 
for  spectrum  analysis  is  conventionally  defined  in  the  narrowband  sense  as 
above,  while  the  additional  effects  associated  with  wideband  signals  are 
classified  under  the  category  of  intersymbol  interference. 

We  thus  begin  by  considering  operation  of  the  integrated  optical 
spectrum  analyzer  of  Fig.  1 for  a single  electronic  frequency  input.  Light 
incident  on  the  surface  elastic  wave  region  is  thus  deflected  at  only  a 
single  angle  and  upon  passage  through  the  transform  lens  arrives  for  the 
most  part  at  a single  sensor  element.  Some  light  will  reach  adjacent 
sensor  elements  due  to  the  inherent  truncation  in  the  Fourier  transform 
associated  with  a finite  optical  aperture.  However,  we  will  assume  for 
the  present  that  for  a suitable  optical  transverse  beam  variation  along 
the  waveguide  surface  (truncated  Gaussian,  Hamming,  Hanning,  etc.),  side- 
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lobe  levels  will  be  minimal  and  thus  essentially  all  of  the  optical 
information  other  than  that  scattered  reaches  a single  sensor  element. 

In  an  actual  integrated  optical  spectrum  analyzer  it  would  be  desirable 
if  scattering  could  be  suppressed  to  such  an  extent  that  dynamic  range 
was  limited  by  the  sidelobe  levels  inherent  in  signal  truncation.  In 
propagating  from  the  Fourier  transform  lens  to  the  detector  array  the 
optical  beam  is  attenuated  due  to  interface  scattering,  bulk  scattering, 
and  bulk  absorption.  As  scattering  occurring  at  the  waveguide  surface 
is  dominant  in  a high  quality  optical  waveguide,  we  assume  that  all  wave- 
guide attenuation  is  due  to  this  contribution.  The  ratio  of  the  total 
scattered  power  to  the  signal  power  reaching  the  single  sensor  element  is 
approximately 

R1  = easL  -1  (1) 

where  a$  is  the  waveguide  attenuation  factor  (power)  due  to  surface  scattering 
and  L is  the  focal  length.  In  (1)  variations  in  path  length  have  been 
neglected,  but  (1)  is  expected  to  be  reasonable  for  small  scattering  and 
moderately  large  f-numbers. 

Of  the  fraction  R^  of  radiation  scattered,  portions  will  be  scattered 
out  of  the  waveguide,  into  higher-order  waveguide  modes,  and  into  the  same 
mode  propagating  in  a different  direction.  A high  quality  waveguide  surface 
is  assumed  so  that  back  scattering  will  be  small  in  comparison  to  forward 
scattering.  Radiative  scattering  will  not  be  collected  by  a sensor  array 
integrated  into  the  waveguide  structure.10  This  is  one  of  the  advantages 
of  such  an  integrated  structure.  We  assume  that  the  waveguide  will  support 
only  a single  transverse  mode  of  a given  polarization  so  that  scattering 
into  higher-order  modes  is  not  possible.  Thus,  only  a fraction  of  radiation 
scattered  into  the  same  waveguide  mode  propagating  in  other  directions  will 
reach  sensor  elements  besides  the  one  corresponding  to  the  signal  input.  To 
determine  dynamic  range,  this  fraction  must  be  determined  for  the  sensor 
element  adjacent  to  the  one  receiving  signal  power. 

As  a first  step  in  this  determination,  the  ratio  R2  of  power  scattered 
into  the  same  waveguide  mode  propagating  in  a different  direction  within 
the  waveguide  to  the  sum  of  this  power  plus  the  scattered  field  radiated 
from  the  guide  having  the  same  k component  along  the  waveguide  surface 
transverse  to  the  incoming  beam  can  be  determined  from  an  intuitive 
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generalization  of  the  coupling  coefficients  determined  by  Marcuse.5*6 
Marcuse's  treatment  of  waveguide  scattering  is  limited  to  two  dimensions, 
ignoring  the  transverse  dimension  y along  the  waveguide  surface  which  is 
of  present  interest.  Scattering  with  propagation  vector  components  in 
this  direction  require  surface  roughness  spectral  components  in  this 
direction  for  phase  matching.  Intuitively,  we  expect  that  the  power 
scattered  with  propagation  vector  components  ky  will  be  proportional  to 
the  ky  component  of  the  spectral  density  of  the  y variation  of  the  surface 
roughness  function.  For  a given  ky  this  factor  should  be  the  same  for 
both  radiation  modes  and  the  single  slab  mode  so  that  R2  is  independent 
of  this  spectral  component  and  determined  directly  from  Marcuse's6  equations 


(9.4.10)  and  (9.4.11)  to  be: 

4(icd/2  , tc/a~)  An 

R = 1 + j f 

ir(n2  + 3)  cos  (Bd-<t>.| ) 


An  + n 


In  (2)  d is  the  waveguide  thickness,  n2  is  the  waveguide  refractive  index, 

n is  the  effective  refractive  index  of  the  lowest  order  TE  mode,  k,  a.,  e 
9 ^ 

and  <j>.j  are  waveguide  mode  parameters  defined  by  the  equations:11 


taned 


e(a1  + a3) 


6 - a1a3 


32  ♦ k2  - k2n2 

-cu  + K2  = k2n2 
3 0 J 

2 . 2 _ .22 

"“I  + K " konl 


tan^  = a.j/6,  tan^^  = a3^6' 

In  (2)  An  is  defined  as 


(3-a) 


(3-b) 

(3-c) 

(3-d) 

(3-e.f) 


where  i is  the  surface  roughness  correlation  length  and  kQ  = 2ir/\  where  x 
is  the  wavelength.  In  (3)  n1  is  the  substrate  refractive  index  while  n3 
is  the  overlay  refractive  index  (assumed  = 1 in  calculations).  In  deter- 
mining (2)  from  Marcuse's  coupling  coefficients,6  we  employed  an  approxi- 
mation to  an  integral  which  sums  contributions  from  all  radiation  modes. 
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We  have  previously  described  details  of  this  approximation. 12  In  deriving 
(2)  we  have  followed  Marcuse  and  assumed  that  the  random  surface  roughness 
can  be  described  by  the  autocorrelation  function 

R(u)  = a2  exp(-  |u|/0  (5) 

2 

where  u measures  the  separation  of  two  points  along  the  surface  and  a is 
the  mean  square  surface  roughness.  The  surface  roughness  spectral  density 
(G(n)l ^ is  the  fourier  transform  of  the  autocorrelation  function,  given  as 

O O 2 . / c\ 


|G(n)| 


An 7 

n2+(An)' 


where  n is  the  surface  roughness  spatial  frequency  normalized  to  kQ. 

Thus  far  we  have  determined  an  expression  R1  for  the  fraction  of  light 
scattered  and  R ^ the  fraction  of  the  total  scattered  field  that  is  scattered 
into  the  same  waveguide  mode  propagating  in  a different  direction.  We  must 
now  determine  the  fraction  R^  of  the  field  scattered  into  the  waveguide  which 
reaches  the  sensor  element  adjacent  to  the  one  receiving  the  signal.  Defining 
0 as  the  scattering  angle  as  measured  from  the  z axis  as  shown  in  Fig.  2,  the 
surface  roughness  spatial  frequency  required  to  phase-match  scattering  at  an 
angle  0 is 

n = 2ng  sin  | (7) 

The  scattering  intensity  of  this  direction  is  proportional  to  the  surface 
roughness  spectral  component  S(0),  where 

S(0)  = G(n)  2 (8) 

keeping  in  mind  that  0 and  n are  related  by  (7).  The  scattered  field  arriving 
at  a given  sensor  element  is  composed  of  the  sum  of  radiation  from  various 
points  on  the  waveguide  surface.  For  each  such  point  the  acceptance  angle  of 
a given  sensor  element  is  different.  With  this  in  mind  the  ratio  of  field 
scattered  to  the  adjacent  sensor  element  to  the  total  field  scattered  in  the 
waveguide  is 

rL 

J J S(0)d0dz 
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In  (9)  [0  (z),  02(z)]  represents  the  acceptance  angle  of  the  adjacent  sensor 
for  scattered  radiation  originating  a distance  z from  the  focal  plane.  The 
y variation  of  this  angular  range  has  been  neglected  as  a first  approximation. 
This  approximation  is  expected  to  be  valid  except  for  small  f-number  lenses. 

At  first  it  would  seem  that  the  integral  in  the  numerator  of  (9)  would  need  to 
be  weighted  to  account  for  the  decreased  width  of  the  beam  approaching  the 
focal  region.  However,  the  field  intensity  increases  so  as  to  exactly  com- 
pensate for  this  effect. 

The  expression  R3  give  in  (9)  can  be  further  evaluated.  Utilizing 
the  geometry  in  Fig.  2,  0^(z)  and  02(z)  can  be  expressed  as 


©2(z ) = tan' 


rns-sQ/2 


<ns+s0/2’ 


(10-a) 


(10-b) 


where  s is  the  detector  array  period,  s is  the  elemental  detector  size, 

o 

and  n is  the  number  of  detectors  away  from  the  signal  detector.  To 
evaluate  dynamic  range  as  we  have  defined  it,  the  n-1  adjacent  detector 
would  be  considered.  Generally,  the  acceptance  angular  interval  [0^(z), 
02(z)]  is  quite  small  so  that  S(  ) can  reasonably  be  considered  constant 
over  this  range.  Thus  (9)  reduces  to 
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where 


0 = tan 
o 
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The  integrand  in  (11)  can  be  approximated  accurately  over  most  of  the  range 
of  z in  the  interval  [0,F]  by  assuming 


tan[02(z)  - o^z)]  « 02(z)  - o^z) 
while  over  the  range  where  (13-a)  is  significant 


tan  0 « 
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Incorporating  the  approximations  (13)  into  (11)  gives  the  result 

-L 


R3  = 


i 


s0  an 


2 2 
An  + 4n 


z2  + 


2 2 
n s 


r 2 2 2 
[z  + n s 


s2/4] 


dz 


(14) 

The  second  term  in  the  integrand  represents  the  effect  of  acceptance  angle 
variation  with  z while  the  first  includes  the  affect  of  the  proportionality 
of  scattering  power  to  the  surface  roughness  spectral  component  and  the 
variation  of  the  appropriate  spectral  component  with  z.  A plot  of  the 
integrand  in  (14)  is  shown  in  Fig.  3 for  the  nth  adjacent  sensor  element 
for  n-1,  2,  and  3.  The  magnitude  of  this  integrand  represents  the  relative 
scattering  power  received  by  the  detector  adjacent  to  the  one  receiving 
the  signal  as  a function  of  the  point  of  origin  of  the  scattered  radia- 
tion. Note  that  there  is  a position  from  which  the  scattering  to  adjacent 
photodetector  elements  originates  which  yields  a maximum  value.  This  is 
due  to  the  variation  in  acceptance  angle  and  by  the  increase  in  the  angle 
of  off-axis  scattering  required  as  the  displacement  between  detector  array 
and  scattering  center  decreases.  Note  that  peaks  corresponding  to  the 
adjacent  detector  elements  decrease  with  successive  elements  and  that  the 
maxima  are  progressively  broadened. 

In  view  of  the  definitions  of  R-j , R2>  and  R3,  the  ratio  of  optical 
power  received  by  the  first  adjacent  detector  next  to  the  signal  receiving 
one  to  the  received  signal  power  is  given  by  R-j R2R3  eacl1  ttiese 

ratios  given  by  (1),  (2),  and  (14).  Dynamic  range  DR  as  limited  by  wave- 
guide scattering  and  defined  earlier  is  thus  given  as 

DR  - -10  log10(R1R2R3).  (15) 

Dynamic  range  (limited  by  waveguide  scattering)  defined  in  (15)  has  been 
calculated  and  results  are  presented  in  the  next  section. 

III.  NUMERICAL  RESULTS 

The  expression  for  dynamic  range  (15)  determined  in  the  previous 
section  has  been  evaluated  as  a function  of  waveguide  attentuation  for 
typical  spectrum  analyzer  parameters  with  the  results  plotted  in  Fig.  4. 

Note  that,  as  expected,  higher  dynamic  ranges  are  associated  with  smaller 
waveguide  attenuation  and  thus  situations  in  which  there  is  less  waveguide 
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acattering.  Fig.  4 indicates  that  the  use  of  smaller  detectors  lalows  a 
larger  dynamic  range  due  to  the  decreased  acceptance  angle  for  scattered 
radiation.  However,  smaller  detectors  generally  have  a smaller  device 
dynamic  range  and  for  such  situations  this  may  be  the  limiting  factor  in 
overall  spectrum  analyzer  dynamic  range. 

For  the  parameters  in  Fig.  4 the  dynamic  range  calculation  has  been 
repeated  for  different  focal  length  lenses.  The  resulting  dynamic  range 
is  nearly  invarient  with  respect  to  focal  length.  For  longer  focal 
lengths  we  would  expect  R3  to  be  smaller  since  the  acceptance  angle  for 
scattering  of  the  adjacent  detector  is  much  smaller  over  most  of  the  range 
of  scattering  centers.  Calculations  agree  with  this  expectation  in  that 
Rg  is  seen  to  decrease  with  increasing  focal  length.  However,  this  decrease 
is  almost  exactly  offset  by  an  increase  in  R^  which  measures  the  ratio  of 
the  total  scattered  field  to  that  received  by  the  detector  element  corres- 
ponding to  the  signal.  The  net  result  is  that  the  product  ^^3  anc*  t*1us 
the  dynamic  range  remains  nearly  invarient  with  respect  to  Fourier  transform 
lens  focal  length. 

In  Fig.  4 dynamic  range  was  plotted  as  a function  of  waveguide 
attenuation  for  a fixed  surface  roughness  correlation  length.  Since  wave- 
guide attentuation  depends  on  both  correlation  length  and  mean  square  roughness, 
varying  the  attenuation  and  keeping  the  correlation  length  constant  implies 
that  the  mean  square  roughness  varies.  Since  these  three  parameters  are 
interdependent,  choosing  one  as  a variable  and  keeping  a second  constant 
implies  that  the  third  varies  with  the  first  in  order  to  maintain  the 
functional  relationship  between  the  three  parameters.  In  Fig.  5,  dynamic 
range  is  plotted  as  a function  of  the  reciprocal  of  the  normalized  surface 
roughness  correlation  length  An,  define  in  (4).  Both  the  case  in  which 
attenuation  is  held  constant  with  mean  square  roughness  varying  to  maintain 
the  functional  dependence  between  it,  attenuation,  and  An,  and  the  case 
in  which  mean  square  roughness  is  held  constant  with  attenuation  variable 
are  presented.  The  limit  of  very  small  correlation  lengths  or  large  values 
of  An  is  applicable  to  point  defects,  the  nearly  constant  dependence  of 
dynamic  range  on  An  for  constant  attentuation  coupled  with  curves  in  Fig.  4 
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clearly  indicates  that  dynamic  range  varies  with  waveguide  attenuation 
and  very  little  with  the  nature  of  the  surface  roughness  which  gives  rise 
to  this  attenuation. 

In  Fig.  6,  the  rejection  ratio  of  adjacent  channels  in  dB  is  plotted 
as  a function  of  the  nth  adjacent  channel.  The  curves  correspond  to  a common 
value  of  waveguide  attenuation  of  .ldB/cm,  implying  that  the  two  curves 
correspond  to  two  different  values  of  mean  square  roughness.  Rjection 
ratio  is  defined  as  the  ratio  of  the  power  reaching  the  intended  detector 
array  element  for  a single  frequency  electronic  signal  input  to  the  scattered 
power  reaching  the  nth  adjacent  detector  array  element.  Rejection  ratio  for 
n=l  is  thus  equal  to  the  negative  of  dynamic  range.  We  note  in  Fig.  6 that 
the  rejection  ratio  drops  faster  for  the  large  correlation  length  case 
(An  = .001).  Surfaces  characterized  by  localized  or  point  defects  thus  give 
rise  to  slower  decreased  in  rejection  ratio  as  compared  to  those  having 
smoother  variations.  Also  shown  in  Fig.  6 is  the  locus  of  the  sidelobe 
maxima  for  a truncated  Gaussian  beam  in  a diffraction  limited  system.  We 
note  that  the  rate  of  fall  off  for  the  long  correlation  length  case  is 
comparable  to  that  for  the  sidelobe  maxima,  whereas  for  the  small  correlation 
length  case  this  rate  is  much  more  gradual  than  that  for  the  sidelobe 
maxima. 

The  previous  calculations  include  only  scattering  contributions 
originating  from  the  waveguide  between  the  Fourier  transform  lens  and 
the  detector  array.  The  question  naturally  arises  as  to  the  effect  of 
scattering  originating  in  the  Fourier  transform  lens,  particularly  in 
view  of  the  fact  that  the  lens  scattering  rate  generally  exceeds  the 
waveguide  scattering  rate.  However,  the  lens  is  further  away  from  the 
detector  array  so  that  detector  element  acceptance  angles  for  lens  scattering 
centers  are  much  smaller  than  for  scattering  centers  in  the  intervening 
section  of  waveguide.  The  effect  of  dynamic  range  is  thus  not  as  much  as 
one  might  initially  expect.  We  have  used  the  same  formalism  previously 
described  to  include  the  effect  of  lens  scattering.  The  waveguide  lens 
region  has  been  modeled  as  a region  having  waveguide  properties  similar 
to  the  surrounding  waveguide  but  with  a higher  attenuation.  This  model 
should  accurately  characterize  the  effects  of  a geodesic  lens,  while,  for 
a Luneburg  lens,  the  model  over-estimates  the  effects  of  scattering  on 
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dynamic  range.  This  over-estimate  occurs  because  part  of  the  light 
scattered  within  the  waveguide  in  the  lens  region  is  scattered  into  higher 
order  modes  which  radiate  into  the  substrate  at  the  termination  of  the 
lens  region,  whereas  this  effect  is  not  included  in  the  calculations. 

Fig.  7 is  a plot  of  dynamic  range  as  a function  of  waveguide  attenuation. 
Curves  are  included  for  no  lens  scattering  (similar  to  previous  results), 
and  for  attenuations  of  1.0,  3.16  and  10.  dB/cm  due  to  lens  scattering. 

As  expected,  at  very  low  values  of  waveguide  scattering  lens  scattering 
dominates  dynamic  range.  Note,  however,  that  for  waveguide  attention  of 
.1  dB/cm,  only  the  10  dB/cm  lens  scattering  case  noticeably  degrades  dynamic 
range  and  only  then  by  3 dB. 

Computer  evaluation  of  in  (11)  and  the  approximation  of  R3  in 
(14)  were  carried  out  for  a wide  variety  of  parameters  w'th  very  accurate 
agreement  found  in  all  cases.  The  approximations  (13)  are  thus  valid  for 
practical  situations. 

IV.  SUMMARY 

An  analysis  of  the  effect  of  optical  waveguide  scattering  in  the 
plane  of  the  waveguide  surface  on  the  dynamic  range  of  the  integrated 
optical  spectrum  analyzer  has  been  performed.  Initially,  only  scattering 
in  the  waveguide  region  between  the  Fourier  transform  lens  and  the  focal 
plane  is  included.  This  is  expected  to  be  valid  since  acceptance  angles 
become  quite  small  for  scattering  centers  outside  this  region.  Numerical 
calculations  of  dynamic  range  as  a function  of  waveguide  attenuation  and 
surface  roughness  properties  were  presented.  The  dynamic  range  is  found 
to  strongly  depend  on  attenuation,  but  to  be  nearly  independent  of  the 
nature  of  the  surface  roughness.  Calculations  indicate  that  dynamic 
range  is  nearly  independent  of  the  focal  length  of  the  Fourier  transform 
lens.  Channel  rejection  in  the  nth  adjacent  channel  for  variable  n is 
determined  and  found  to  fall  off  quite  slowly  with  n.  The  effect  of  lens 
scattering  is  numerically  evaluated  and  shows  that  lens  scattering  is  only 
important  for  cases  in  which  waveguide  scattering  is  very  small  and  lens 
scattering  is  very  high.  This  calculation  confirms  the  expectations  that 
only  scattering  from  the  waveguide  region  between  the  Fourier  transform 
lens  and  the  focal  plane  is  of  major  significance. 
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Waveguide  Lenses 


Figure  2 


Detector  Element 


Photodetector  Array  Geometry  of  Spectrum  Analyzer  in  the 
Spatial  Frequency  Plane. 
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Surface  Roughness  g X = 9 pm  Detector  Spacing  * 32  pm 

Correlation  Length  = ft  pm  Detector  Size  - 25  pm 


Figure  3 Scattered  Power  in  Adjacent  Detectors  as  a Function  of  the 
Separation  of  the  Scattering  Source  from  the  Detector  Plane. 
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Figure  5 Dynamic  Range  as  a Function  of  An  = X/(2tt«.)  for  Constant 
Attenuation  and  Variable  Mean  Square  Surface  Roughness 
and  for  Constant  Mean  Square  Surface  Roughness  and 
Variable  Attenuation. 
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ABSTRACT 

A "first-cut"  approach  to  a wideband  real-time  RF  spectrum  analyzer 
is  described  using  an  integrated  optical  planar  waveguide  circuit  contain- 
ing Luneburg  expansion  and  transform  lenses,  a surface  elasto-optic  spa- 
tial modulator,  and  a waveguide-coupled  photodetector  array-all  formed 
on  a single  thermally-oxidized  silicon  substrate  with  a hybrid-coupled 
GaAs  laser. 

SUMMARY 

An  optical  approach  to  RF-spectral  analysis  requires  a coherent 
source,  appropriate  expansion  and  transform  lenses,  a spatial  modulator, 
and  photodetector  array  in  the  frequency  |j>lane.  Two  integrated  optical 
approaches  are  being  pursued  by  industry.  * They  are  based  upon  (1)  a 
hybrid  coupling  approach  for  the  GaAs  laser  and  the  detector  array  to  a 
diffused  waveguide  surface  on  a LiNbOg  substrate  containing  geodesic 
lenses  and  a surface  elastic  wave  Bragg-type  modulator  and  (2)  the  use 
of  a thermally-oxidized  silicon  substrate  employing  a multiplicity  of 

3 

overlay  films  for  the  planar  waveguide,  assorted  lenses  , and  a piezo- 
electric to  launch  surface  elastic  waves  and  with  the  oxide  removed  to 

4 

form  the  detector  array  in  the  substrate.  A "first-cut"  at  the  inte- 
gration of  all  these  components  on  a silicon  substrate  is  delineated. 


* This  research  was  supported  in  part  by  the  Air  Force  Avionics  Laboratory. 
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The  thermally-oxidized  silicon  substrate  provides  a smooth  low- 
loss  isolation  layer  for  a Corning  7059  waveguide  extending  throughout 
the  circuit  area  except  where  the  oxide  is  removed  for  optical  coupling 
into  the  detector  junctions,.  The  GaAs  laser  is  coupled  to  the  planar 

5 

waveguide  by  a butt-end  technique  using  a cleaved  edge. 

Thin  Ta205  overlay  films  are  employed  to  form  Luneburg  waveguide 
lenses  used  for  transverse  expansion  of  the  optical  beam  to  encompass 
the  spatial  modulator  time-aperture  and  to  refocus,  forming  the  Fourier 
transform.  These  Luneburg  lenses  are  circular-symmetric,  inhomogeneous 
and  are  void  of  all  aberrations  except  for  field  curvature.  They  are 
formed  by  RF-sputtering  through  a circular  mask  where  the  mask  edge 
shape  controls  the  overlay  film  radial  thickness  profile.  A synthesis 
technique  has  been  developed  to  calculate  the  mask  shape,  the  resulting 

7 

film  deposition  profile,  followed  by  ray  tracing  through  the  lens  , to 
yield  diffraction-limited  results  which  have  been  confirmed  experimentally 
for  F/4  lenses. 

8,9 

The  spatial  modulator  employs  a transverse  Bragg-type  surface 
elasto-optic  interaction  using  an  Nb^Og  overlay  waveguide.  Surface  waves 
are  launched  by  using  an  RF-sputtered  ZnO  piezoelectric  film  over  an 
interdigital  transducer  array  in  the  transducer  section  only  with  tapered 
edges  for  coupling  into  the  surface  elasto-optic  section.  The  highly 
oriented  ZnO  overlay  film  yields  an  electromechanical  coupling  coefficient 

p 10 

k = 0.015.  To  obtain  wide  bandwidth  (310-540  MHz),  both  beam  steering 

11,12 

and  angular  tilting  array  techniques  have  been  employed. 

13 

The  output  waveguide-coupled  photodetector  array  in  the  frequency 
plane  with  a linear  density  period  of  6 microns  provides  for  a frequency 
resolution  less  than  1 MHz.  An  "all  dry"  processing  approach  has  been 
employed  for  the  detector  using  UV  interferometric  fringes  to  generate  the 
detector  array,  waveguide  coupling,  and  metallic  contact  in  a single  litho- 
graphic step. 

The  use  of  a silicon  substrate  with  compatible  thin  films  and  pro- 
cessing procedures  familiar  to  the  semiconductor  industry  represents  a 
potential  significant  low-cost  approach  to  RF-spectral  analysis.  Waveguide 
losses  and  scattering  impose  an  ultimate  limit  upon  the  dynamic  range  and 
resolution.  Surface  elastic  absorption  across  the  spatial  modulator  optical 


aperature  imposes  a high  frequency  roll-off  which  affects  the  time-bandwidth- 

product  vfhich  may  be  partially  compensated  for  by  equalization. 
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APPENDIX  E 


INTEGRATED  OPTICAL  SPECTRUM  ANALYZER  SIGNAL  LEVEL  REQUIREMENTS 

By 

Professor  J.  T.  Boyd 
University  of  Cincinnati 

I.  INTRODUCTION 

o g 

Wideband  spectrum  analysis  of  electronic  signals  in  the  10  - 10  Hz 

frequency  range  utilizing  an  integrated  optical  approach  has  been  introduced 

1 2 3 

and  analyzed  in  previous  papers.  * ’ Recent  research  in  this  area  has  been 
directed  toward  utilization  of  an  Si  02/Si  substrate  for  fabrication  of 
optical  waveguide  lenses,  surface  elastic  wave  transducers  and  waveguide 
coupled  linear  detector  arrays  with  charge-coupled  devices  (CCDs).  This 
substrate  coupled  with  a semiconductor  laser  would  then  comprise  the  spectrum 
analyzer. 

The  use  of  Si  02/Si  offers  a number  of  advantages.  For  the  detector 
portion  of  the  device,  its  use  allows  fabrication  of  sophisticated  CCD  readout 
or  photodiode-amplifier  detector  arrays  which  provide  both  good  performance 
and  in  the  case  of  CCD  arrays  convenient  serial  signal  output.  Furthermore, 
such  detector  structures  can  be  integrated  into  the  optical  waveguide  structure. ^ 
This  integration  minimizes  scattering  which  occurs  in  waveguide-detector 
coupling^  ancj  thus  allows  superior  performance  as  compared  with  the  use  of 
an  external  detector  array. 

The  present  analysis  considers  signal  levels  for  which  a detector 
array  is  expected  to  operate.  An  array  of  photodiodes,  possibly  coupled 
to  individual  integrated  amplifiers,  provides  the  ultimate  in  spectrum 
analyzer  signal  output;  that  is,  instantaneous,  real-time  signal  output  from 
each  spectrum  analyzer  spectral  point.  We  thus  consider  the  signal  handling 
and  noise  properties  of  this  type  of  detector  in  Section  II  to  follow.  There 
are  a number  of  anticipated  applications  of  the  integrated  optical  spectrum 
analyzer,  however,  in  which  it  will  be  advantageous  to  sacrifice  some  of 
the  advantages  of  real-time  signal  output  for  each  spectral  channel  for  the 
convenience  of  spectral  signal  output  in  a serial,  time-multiplexed  format. 
Conversion  of  spectral  information  in  a parallel  format  represented  as  signal 
levels  along  an  array  of  photodiodes  to  a serial  format  is  simply  effected 
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by  the  incorporation  of  a CCD  into  the  detector  array.6  An  integrated  optical 
waveguide-CCD  image  array  suitable  for  the  above  purpose  has  previously  been 
demonstrated. ^ Detailed  analysis  of  the  signal  handling  capability  and  noise 
properties  of  this  device  are  considered  in  Section  III. 

The  knowledge  determined  in  Sections  II  and  III  concerning  signal  ranges 
over  which  the  detector  arrays  operate  allows  estimation  in  Section  IV  of 
the  incident  optical  power  required  to  fully  utilize  detector  dynamic  range. 
Estimation  of  detector  quantum  efficiency,  waveguide-detector  coupling 
efficiency,  scattering  losses,  deflection  efficiency,  and  laser  coupling 
efficiency  then  allow  calculation  of  the  laser  power  required  to  fully  utilize 
detector  dynamic  range. 

II.  INTEGRATED  OPTICAL  PHOTODIODE  DYNAMIC  RANGE 

To  determine  the  dynamic  range  of  a photodiode,  both  the  maximum  signal 
handling  capability  and  noise  level  of  the  photodiode  must  be  determined.  A 
photodiode  is  normally  operated  in  a series  circuit  with  a bias  voltage  and 
a load  resistor.  The  load  resistor  may  correspond  to  the  input  impedance  of 
an  amplifier.  The  optical  signal  is  measured  by  the  current  I flowing  in 
the  series  circuit,  given  as 


I = -I0  [exp(-eVr/kT)  - 1]  + IL 


(1) 


where  1^  is  the  current  flowing  due  to  the  presence  of  light  incident  on 
the  photodiode  I0  is  the  pn-junction  reverse  bias  saturation  current,  Vr  is 
the  reverse  voltage  appearing  across  the  pn-junction,  e is  electron  charge, 
k is  Boltzman's  constant,  and  T is  temperature.  The  voltage  Vr  is  related 
to  the  bias  voltage  V0  by 

Vr  = V0  - IRL  , (2) 


where  R|_  is  the  load  resistance.  The  current  Ij_  is  proportional  to  incident 
light  power  as: 


(3) 
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where  P is  the  incident  optical  power,  hv  is  the  photon  energy,  and  nq  is 
the  photodiode  quantum  efficiency. 7 Thus,  for  a wide  range  of  currents 
such  that  the  exponential  term  in  Eq.  (1)  is  negligible,  the  total  current 
is  also  proportional  to  light  power  through  Eq.  (1).  But  at  large  light 
levels,  the  current  flowing  increases  such  that  the  exponential  term  in 
Eq.  (1)  becomes  significant.  At  this  point  Eqs.  (i)  and  (2)  form  a tran- 
scendental equation  which  must  be  solved  for  I.  Solutions  indicate  that 
the  total  current  saturates  for  increasing  incident  light  power.  At  the  onset  of 
saturation  a departure  from  linearity  in  photodiode  response  thus  occurs. 

We  define  the  maximum  optically  excited  current  II  as  that  value  corres- 
ponding to  a 1%  deviation  in  response  linearity.  This  value  of  current 
is  a solution  of 


h max  = WrL  + (Vth/R|.)  ln  0°'2  *L  max/1©)  = VRL  . (4) 


where  Vth  = kT/e  is  the  thermal  voltage.  The  value  of  R|_  is  usually 
related  to  the  electrical  bandwidth  Af  of  the  photodiode  by 


Af 


_J 

2irR|_Cc| 


(5) 


where  Cj  is  the  junction  capacitance  of  the  photodiode.  Substituting  Eq.  (5) 
into  Eq.  (4)  gives 


max  = 27rV0Cc|Af  . (6) 

Thus,  the  maximum  photodiode  response  is  proportional  to  the  bias  voltage 
and  is  inversely  proportional  to  load  resistance  or,  equivalently,  proportional 
to  junction  capacitance  and  detector  bandwidth. 

Now  that  the  maximum  photodiode  photocurrent  has  been  determined, 
determination  of  dynamic  range  requires  calculation  of  photodiode  noise 
current.  Photodiodes  are  characterized  by  both  signal  shot  noise  and  thermal 
noise  associated  with  Rl-  Because  shot  noise  is  much  less  than  thermal  noise 
at  minimum  detectable  incident  optical  power  levels  and  because  it  is  signal 
dependent,  it  is  not  included  in  dynamic  range  calculations.  The  rms  current 
due  to  thermal  noise  (T^)*5  is 
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I 


(Ith),S  = = Af  ^8nCdkT  . (7) 

This  is  the  minimum  detectable  photocurrent  since  smaller  values  of  photo- 
current cannot  be  distinguished  from  noise. 

Dynamic  range  is  defined  in  the  usual  way  as  the  ratio  of  the  maximum 
signal  level  to  the  minimum  signal  level 


DR 


ratio 


*L  max 
min 


(8) 


where  Eqs.  (6)  and  (7)  were  used.  In  terms  of  dB  the  electrical  dynamic  range 
applicable  to  photodiode  current  is 


(9) 


while  the  optical  imaging  dynamic  range  is 


00) 


As  deflected  optical  field  is  linearly  proportional  to  acoustic  field  in  the 
range  of  acoustic  power  applicable  to  the  spectrum  analyzer,  the  optical 
dynamic  range  DR0  characterizes  the  spectrum  analysis  process.  Note  that 
photodiode  dynamic  range  is  independent  of  bandwidth.  We  consider  a numerical 
example  similar  to  our  integrated  optical  photodiode  array  previously 

O 

fabricated.  The  pn  photodiode  dimensions  can  be  expected  to  approach 
7 urn  x 50  urn  with  a substrate  impurity  concentration  of  1015  cm"^.  For 
V0  = 100  volts,  the  depletion  region  width  is  11.5  um,  giving  a junction 
capacitance  of  Cd  = 3.2  x 10'^  f.  For  V0  = 100  volts  spectrum  analyzer 
dynamic  range  is  DR0  = 50.4  dB.  This  value  is  comparable  to  that  achiev- 
able with  a CCD  image  array,  as  will  become  evident  from  the  calculations  in 
the  following  section. 
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III.  INTEGRATED  OPTICAL  CCD  IMAGE  ARRAY  DYNAMIC  RANGE 

The  CCD  image  array  integrated  into  an  Si02/Si  optical  waveguide  structure 
which  has  been  demonstrated^  is  shown  in  Fig.  1.  This  particular  device 
configuration  provides  for  charge  excitation  in  the  photodiode  regions  with 
the  gate  potentials  adjusted  so  that  photo-excited  charge  accumulates  under 
the  integration  gate  for  each  array  element.  The  presence  of  the  photo- 
diode provides  an  electrode-free  surface  for  waveguide-detector  coupling. 

With  the  presence  of  the  integration  gate,  this  device  configuration  is 
applicable  to  operation  of  the  integrated  optical  spectrum  analyzer  in 
either  a pulsed  or  an  integrating  mode.  Charge  integrates  in  the  potential 
well  formed  underneath  the  integration  gate  in  the  device  shown  in  Fig.  1.  At 
the  end  of  the  integration  period  T0,  the  transfer  gate  is  turned  on  sc  that 
the  charge  accumulated  in  each  array  element  is  transferred  simultaneously  into 
the  CCD  shift  register  portion  of  the  device.  This  parallel  transfer  of 
charge  occurs  in  a time  less  than  1/2  of  the  CCD  clock  period  so  that  the 
CCD  is  stationary  during  this  time.  Once  this  charge  has  transferred  into 
the  CCD,  the  transfer  gate  is  turned  off  and  the  charge  is  transferred  along 
the  CCD  to  the  end  where,  after  suitable  amplification,  a serial,  discrete- 
time signal  emerges  representing  the  spatial  light  distribution.  As  soon  as 
the  transfer  gate  is  turned  off,  charge  integration  again  begins  and  the 
cycle  is  repeated.  The  integration  time  is  thus  given  as 

To  = Nx  = N/fc  , (11) 

where  N is  the  number  of  CCD  stages  equal  to  the  number  of  array  elements, 
fc  is  the  CCD  clock  frequency  and  x is  the  CCD  clock  period  equal  to  the 
time  required  to  transfer  a charge  packet  from  one  stage  to  the  next  adjacent 
stage. 

We  shall  consider  the  following  three  cases  of  spectrum  analyzer 
operation: 

Case  I - This  is  a pulsed  mode  in  which  the  pulse  duration  is 
much  less  than  the  surface  elastic  wave  transit  time 
Te  across  the  optical  aperature.  Pulses  occur  once 
each  Te  seconds  occurring  when  the  signal  of  interest 
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is  positioned  in  the  optical  aperture.  For  a large 
number  of  spectral  points  (large  time-bandwidth  products), 
the  transit  time  Te  may  not  provide  ample  time  to 
completely  empty  a single  CCD  register.  For  such  situ- 
ations, we  shall  assume  a CCD  having  several  sections, 
each  with  separate  output  aligned  one  after  another 
along  the  Fourier  transform  plane. 

Case  II  - An  integrating  mode  in  which  the  integration  period  is 
on  the  order  of  the  transit  time  Te.  The  laser  operates 
CW  and  we  shall  use  an  integration  time  Te  in  our 
calculations.  As  in  the  previous  case,  the  integration 
time  Te  may  not  provide  ample  time  to  completely  empty  a 
CCD  register  corresponding  to  a large  number  of  spectral 
points.  For  such  situations  we  shall  assume  a CCD  having 
several  sections  and  several  separate  outputs. 

Case  III  - An  integrating  mode  in  which  the  integration  period  is  long 
compared  to  the  transit  time  Te.  The  laser  operates  CW, 
and  we  shall  assume  that  the  integration  time  is  determined 
as  the  time  required  to  shift  out  a single  section  CCD 
imaging  device. 

These  three  cases  are  applicable  to  a wide  variety  of  situations. 

In  what  follows,  calculations  will  be  presented  for  a Surface  Channel  CCD 
(SCCD)  similar  to  the  one  previously  fabricated^  as  opposed  to  a Bulk  Channel 
CCD  (BCCD).  Generally,  SCCD's  have  larger  signal  handling  capabilities  and 
larger  dynamic  ranges  than  BCCD's.  However,  BCCD's  are  characterized  by 
lower  noise  levels  and  a larger  maximum  operating  frequency.  The  calculations 
which  follow  could  be  altered  to  apply  for  a BCCD  by  replacing  fast  interface 
state  noise  applicable  to  SCCD's  by  a bulk  state  trapping  noise  and  replacing 
the  maximum  charge  handling  capability  of  a surface  channel  by  that  of  a bulk 
channel . 

To  calculate  the  CCD  image  array  dynamic  range,  the  maximum  signal 
handling  capability  and  noise  level  of  the  device  must  be  determined.  The 
maximum  signal  handling  capability  of  a charge  transfer  device  is  limited  by 
the  breakdown  of  the  Si O2  insulator  region.6  Gauss'  law  then  implies  that 
the  total  charge  e nmax  associated  with  the  maximum  electric-field  which 
the  Si02  region  can  withstand  is  given  as 
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e n 


max 


e e E A 
o ox  max 


(12) 


where  nmflX  is  the  maximum  number  of  charge  carriers,  e is  the  electron  charge, 

A is  the  area  of  the  electrode  in  question,  e0  is  the  free  space  permittivity 
(e0  = 8.85  x 10“14  farad/cm.),  and  eox  is  the  dielectric  constant  of  Si O2 
(eox  w 4.0).  Because  the  area  of  the  integration  gate  Aig  is  slightly  larger 
than  the  area  of  the  CCD  electrodes  Asg  for  the  devices  previously  fabricated,4 
the  maximum  charge  which  can  be  handled  by  the  overall  device  is  limited  to 
that  which  can  be  handled  by  the  CCD,  i.e.. 


Vox 


max 


E A 
max  sg 


(13) 


The  value  of  varies8^  over  a range  1 - 6 x 10^  v/cm. , depending  on  a 

ITiaX 

variety  of  processing  and  material  conditons.  For  the  present  calculations, 
we  shall  choose  the  value 


E«,  s 41  '°6  v/c,n-  • 


as  the  maximum  allowable  field.  Using  this  value  in  Eq.  (13)  gives 


(14) 


nmax^Asg  = 8,85  x 1q12  cm’2  * O5) 

This  value  is  more  conservative  than  that  given  by  Sequin  and  Tompsett.6 
The  electrodes  of  the  previously  fabricated  devices  have  a typical  size  of 
10  urn  x 125  pm,  giving  Asg  = 125  x 10"5  cm2.  Using  this  value  in  determining 
nmax  ^rom  a^ove  gives 

nmax  * * ’°8  ' <’6> 

The  value  in  Eq.  (15)  is  much  greater  than  values  of  the  surface  state  density 
Qss/e  for  the  Si 02/Si  interface,  justifying  the  omission  of  Qss  in  Eq.  (13). 

Now  that  the  maximum  signal  handling  capability  has  been  calculated, 
determination  of  dynamic  range  requires  calculation  of  the  noise  levels 
present  in  the  CCD  device  in  Fig.  1.  Noise  in  this  device  arises  from  a 
variety  of  sources  including  output  amplifier  noise,  noise  associated  with 
dark  current  fluctuations,  noise  associated  with  signal  input  circuitry. 


noise  associated  with  fluctuations  in  the  emission  and  trapping  of  carriers  by 
fast  interface  states,  and  noise  associated  with  fluctuations  in  bias  charge 
transfer  loss.  These  sources  of  noise  shall  be  considered  separately  in 
what  follows  with  the  net  effect  of  noise  eventually  determined.  As  the 
various  sources  of  noise  are  independent,  their  mean  square  fluctuations  in 
the  number  of  carriers  will  add  to  give  the  cumulative  effect.  In  addition 
to  the  above  sources  of  noise,  shot  noise  associated  with  the  input  optical 
signal  will  also  give  rise  to  signal -dependent  fluctuations  in  both  the 
imperfect  charge  excitation  and  transfer  processes.  However,  as  the  rms 
fluctuation  in  the  number  of  carriers  is  equal  to  the  square  root  of  the 
average  number  of  signal  carriers,  this  fluctuation  will  be  negligible  at 
minimum  detectable  signal  levels.  This  source  of  noise  will  thus  not  be 
included  in  the  dynamic  range  calculation. 

For  surface  channel  CCD's,  fast  interface  state  trapping  and  emission 
can  seriously  degrade  charge  transfer  efficiency.  This  is  particularly 
true  after  operation  for  a period  of  time  with  no  input  signal.  For  such 
situations  these  traps  empty  and  thus  when  signal  charge  does  appear  a 
significant  amount  of  it  is  trapped.  This  problem  can  be  overcome  by 
injecting  a continuous  small  bias  charge  referred  to  as  fat  zero  so  that 
trap  occupation  and  thus  free  charge  remain  constant  on  the  average.  Emission 
and  trapping  of  charge  continue  to  occur  but  at  balancing  rates.  Although 
this  technique  significantly  improves  transfer  efficiency,  fluctuations  in 
the  trapping  and  emission  of  carriers  is  still  present.  In  general,  these 
fluctuations  are  independent  for  each  storage  electrode  and  thus  mean  square 
fluctuations  add.  The  rms  number  of  noise  electrons  rTjs  due  to  fast  interface 
state  fluctuations  for  fat  zero  operation  has  been  shown  to  be^0’^ 

"is-  VkT  Nss  in  2 (2  NAJg  ♦ A,g)  , (17) 

where  k is  Boltzman's  constant,  T is  absolute  temperature,  Nss  is  the  fast 
interface  state  density  in  ev^cm-2,  and  A^g  is  the  area  of  the  integration 
gate.  Note  that  although  charge  from  all  but  one  of  the  parallel  channels 
travels  less  than  N stages,  fast  interface  state  noise  accumulates  for  all 
N stages  in  each  signal  charge  packet.  Although  values  of  Nss  as  low  as 
109ev'^cm"^  have  been  observed,  we  employ  a more  typical  value^  of  5xl09 
ev^cnf2  in  Eq.  (7)  and  find 
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08) 


n.$  = 47  /N  + 1.1 

where  the  area  of  the  integration  gate  is  taken  as  Ajg  = 2.87  x 10"^citi"^. 

The  noise  contribution  of  fast  interface  states  as  given  by  Eqs.  (17)  and 
(18)  is  an  upper  limit  as  this  noise  is  suppressed13  at  low  CCD  signal 
frequencies.  Low  CCD  signal  frequencies  correspond  to  analysis  of  broad  band 
input  signals.  Narrow  band  input  signals  after  undergoing  spectrum  analysis 
give  rise  to  CCD  signal  variations  approaching  the  Nyquist  limit  so  that  no 
suppression  of  the  above  noise  contributions  is  expected  for  such  input 
signals. 

Output  amplifier  noise  generally  makes  a significant  contribution  to  over- 
all device  noise.  Because  charge  packets  in  CCD's  correspond  to  low-level 
electrical  signals,  amplification  is  required.  An  MOS  reset  amplifier  inte- 
grated onto  the  same  chip  as  the  CCD  has  been  found  to  be  most  satisfactory 
for  this  amplification.  The  reset  process,  however,  introduces  thermal  noise 
generally  on  the  order  of  200  carriers. ^ White,  et.  al.,15  introduced  the 
correlated  double  sampling  technique  to  remove  this  noise,  however.  We  shall 
assume  that  correlated  double  sampling  is  employed  so  that  this  noise  is 
removed  from  the  CCD  of  present  interest.  The  remaining  primary  contribution 
to  output  amplifier  noise  is  thermal  noise,  as  1/f  noise  is  suppressed  at  CCD 
operating  frequencies  of  present  interest.  Careful  design  of  the  output 
amplification  circuitry  has  yielded  a device  having  a measured  rms  noise 
conti rubti on  of  27  carriers  compared  to  a theoretical  design  of  25.  We  shall 

choose  a more  conservative  value  for  the  rms  carrier  fluctuation  n„,  of 

oa 

"oa  * 100  ’ (,9> 


arising  from  the  output  amplifier. 


Noise  associated  with  dark  current  generation  in  CCD's  is  particularly 
significant  in  imaging  devices.  Fluctuations  in  dark  current  is  a shot  noise 
process  so  that  the  rms  number  of  noise  carriers  associated  with  these  fluctu- 
ations n^  is  equal  to  the  square  root  of  the  mean  number  of  carriers  generated; 
i.e. , 


nd 


1 e fc 


(20) 
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In  Eq.  (20)  J0  is  the  dark  current  level  which  will  be  taken  to  be  J0  = 6 x 10"^ 
A/cm2  corresponding  to  measurements,14  although  lower  values  have  been 
reported.12  Using  a clock  frequency  of  10  MHz,  Eq.  (20)  becomes 


nd  = .39  /N 


(21) 


As  in  the  case  of  fast  interface  state  noise,  Eq.  (20)  includes  the  fact  that 
dark  current  accumulates  throughout  the  entire  device  even  though,  with 
parallel  entry,  most  signal  charge  is  transferred  over  only  a portion  of  the 
device. 

Input  thermal  noise  is  significant  both  from  the  parallel  injection  of 
optically  excited  charge  and  from  serial  injection  of  a bias  charge  or  fat 
zero.  Again,  1/f  noise  is  not  significant  at  the  high  operating  frequencies 
of  present  interest.  In  general,  the  rms  number  of  carriers  n^  associated 
with  input  thermal  noise  is  given  as®»14 


in 


kT(a.  C.  + a-  C.  ) 
in  in  lg  ig 


(22) 


where  C-n  is  the  input  gate  capacitance,  C.  is  the  integration  gate  capacitance 
and  ain  and  a.j  are  parameters  of  the  order  of  unity  which  depend  on  the 
specific  input  conditions.  We  shall  assume  both  equal  one.  The  input 
thermal  noise  due  to  fat  zero  injection  can  be  reduced  by  reducing  the  width 
of  the  serial  input  gate.14  This  reduces  the  signal-handling  capability  of 
the  input  gate  but,  since  only  low-level  bias  charge  is  injected  serially, 
it  is  of  no  consequence.  For  a serial  input  gate  having  dimensions  12  urn  x 
10  pm.  the  capacitance  of  the  input  gate  and  the  integration  gate  are: 


Cin 


e e A. 

= .04  pf 

Lox 


(23a) 


EE  A . 

■ °vox-  ig  = 1.02  pf  , 

xox 


(23b) 


where  we  assumed  the  oxide  thickness  tQX  = .1  pm.  Using  these  values  of 
capacitance  in  Eq.  (22)  gives 

"in  * 4,4 

for  the  rms  number  of  input  noise  carriers. 
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(24) 


The  significant  sources  of  noise  have  now  been  quantified.  As  these 
sources  are  independent,  their  mean  squares  add.  Denoting  mean  square  values 
as  nx2,  where 


and  x is  any  of  the  previous  subscripts,  the  total  rms  number  of  noise  carriers 


- _ 2 2 . 2 . 

nt  ' >nis  + noa  + nd  + nin  ' 


Now,  using  Eqs.(18),  (19),  (21),  and  (24)  in  Eq.  (26)  gives 

1 

nt  = [ (414)2  + (.39^  )2  + (47/5+T  )2  + (100)2]  2 (27) 

The  total  noise  level  nt  can  now  be  used  to  calculate  dynamic  range. 

Dynamic  range  is  defined  as  before  as  the  ratio  of  the  maximum  signal 
level  to  the  noise  level;  i.e.. 


ratio 


= .9 


where  Eqs.  (16)  and  (27)  need  to  be  used  for  evaluation.  The  factor  .9  is 
introduced  to  account  for  fat  zero  utilization  of  10%  of  the  maximum  device 
signal  handling  capability.  In  terms  of  dB,  the  electrical  dynamic  range 
(CCD  current)  and  the  optical  imaging  dynamic  range  are 


DRe  = 20  log10  .9 


DRq  - 10  log1Q 


As  noted  earlier,  the  optical  dynamic  range  DRQ  is  the  one  characterizing 
operation  of  the  integrated  optical  spectrum  analyzer.  Table  I presents 
this  dynamic  range  information  for  each  of  the  three  cases  discussed  earlier 
as  a function  of  the  number  of  resolvable  spectral  points  in  the  signal 
Fourier  transform. 
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TABLE  I CCD  DYNAMIC  RANGE 
f = 10  MHz  10%  Fat  Zero 

t 


Case  I 
Pulsed  Mode 

Pulsewidth  «Transit  Time 


Case  II 
CW  Mode 

Integration  Time= 
Transit  Time 


Case  III 
CW  Mode 

Integration  Time> 
Transit  Time 


No.  of 
Spectral 
Points 


CCD 

Optical 

Dynamic 

Range 

(dB) 


No.  of 
Separate 
Outputs* 
Required 


CCD 

Optical 

Dynamic 

Range 

(dB) 


No.  of 
Separate 
Outputs* 
Required 


CCD 

Optical 

Dynamic 

Range 

-(dB?. 


No.  of 
Separate 
Outputs* 
Required 


32 

53 

1 

53 

1 

53 

1 

64 

53 

2 

53 

2 

52 

1 

128 

53 

4 

53 

4 

52 

1 

256 

53 

8 

53 

8 

51 

1 

512 

53 

16 

53 

16 

49 

1 

1024 

53 

32 

53 

32 

48 

1 

♦Fewer  outputs  would  be  required  for  those  cases  for  operation  at  higher 
clock  frequencies.  The  dynamic  range  varies  slowly  with  frequency  in  this 
frequency  range. 

IV.  INTEGRATED  OPTICAL  FOURIER  TRANSFORM  LASER  POWER  REQUIREMENTS 

The  determination  of  the  range  of  signal  levels  over  which  the 
detector  operates  allows  determination  of  the  laser  power  required  to 
fully  utilize  detector  dynamic  range.  If  laser  power  PQ  is  coupled  into 
the  waveguide  with  coupling  efficiency  nlc,  the  initial  power  in  the 
optical  waveguide  is 


Vo 


(31) 


This  optical  power  undergoes  attenuation  at  a rate  ag  as  it  propagates  in 
the  optical  waveguide,  attenuation  at  a rate  a1  as  it  propagates  in  lens 
regions,  reflection  losses  at  reflectors  (r  = product  of  reflectivities  at 
each  reflector),  and  transmission  losses  at  waveguide-lens  and  other  dis- 
continuities (t  = product  of  transmissivities  at  each  such  discontinuity). 
The  resulting  undeflected  optical  power  in  the  transform  plan  is  thus 


rt 


^lc'o 


(32) 


where  L1  and  Lg  are  the  mean  propagation  lengths  in  the  lens  and  waveguide 
regions,  respectively. 

In  the  operation  of  the  spectrum  analyzer,  a portion  of  this  optical 
power  is  deflected  to  the  detector  area.  We  define  as  the  fraction  of 
optical  power  deflected  with  nj  max  as  the  maximum  deflection  efficiency, 
corresponding  to  maximum  elastic  wave  power.  For  the  present  calculation,  we 
consider  single  frequency  surface  elastic  wave  excitation  having  maximum 
intensity,  so  that  all  deflected  optical  power  (except  side  lobe  contributions) 
is  incident  on  a single  detector  element.  The  maximum  deflected  optical  power 
Pine  corresponding  to  maximum  elastic  wave  power  incident  on  a single 
detector  element  is  thus 


inc 


nd  max  nsl 


-«iLi 


’agLg 


rt 


^o 


(33) 


where  nsi  is  the  fraction  of  optical  power  in  the  main  lobe  of  the  focused 
optical  beam.  The  maximum  elastic  wave  intensity  should  ideally  have  a 
value  which  yields  maximum  deflection  efficiency  while  maintaining  linearity 
between  the  deflected  optical  field  and  the  elastic  wave  field.  Criteria 
for  the  maintenance  of  linearity  are  somewhat  arbitrary.  For  the  present 
we  impose  the  criteria  that  the  deviation  from  linearity  should  not  exceed 
1%.  The  corresponding  deflection  efficiency  corresponding  to  the  maximum 
elastic  wave  power  should  thus  be  n<j  max  = -06. 

For  photodiode  detection  the  photodiode  current  corresponding  to  the 
incident  power  in  Eq.  (33)  is 


i.  = f!!a 

inc  hv 


inc 


(34) 


similar  to  Eq.  (3).  As  the  incident  deflected  power  Pinc  corresponds  to  the 
maximum  optical  power  which  would  ever  reach  a detector  element,  the  current 
line  should  equal  the  maximum  current  I|_max  from  Eq.  (6)  which  a photodiode 
can  pass  while  remaining  within  1%  of  photoresponse  linearity  in  order  to 
fully  utilize  detector  dynamic  range.  Thus,  imposing  the  condition 


*inc  = max 


(35) 


71 

\ 
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in  Eq.  (34)  and  using  Eq.  (33)  in  Eq.  (34),  we  can  solve  for  the  required 
laser  power  PQ  to  be 

a-iL-,  a L 

27rhvV  C jAf  e ' e 9 9 

p = (3, 

0 eVdmaxns1nlc  rt 

Equation  (36)  has  been  evaluated  for  the  following  parameters: 


X = .9  pm  (hv  = 2.21  x 10’19  j) 

a-j  = .46  cm"^  (2  dB/cm.  attenuation  in  lens  region) 

ctg  = .115  citH  (.5  dB/cm.  waveguide  attenuation) 

l_1  = 2.5  cm  (propagation  length  in  lens  region) 

Lg  = 10  cm  (propagation  length  in  waveguide  region) 

Oq  = .8  (waveguide  detector  quantum  efficiency  = .8  has  been  measured5) 

ndmax  = (maximum  deflection  efficiency  which  maintains  nonlinear 
deviation  < 1%) 

nsl  = *®5  (estimated  fraction  of  power  in  main  lobe) 

n]c  = -45  (laser-waveguide  coupling  efficiency  = .8  has  been  experimentally 
real ized16) 

r = .9  (net  fractional  reflector  losses  = 1 - r) 
t = .9  (net  fractional  transmission  losses  = 1 - t) 

V0  = 10  volts  (photodiode  bias) 

C,j  = 3.2  x ICT^f  (photodiode  junction  capacitance  corresponding  to 
photodiode  discussed  in  Section  II) 


The  only  parameter  in  Eq.  (36)  not  specified  above  is  the  photodiode  bandwidth 
Af.  For  real-time  output  the  spectrum  analyzer  can  only  resolve  signal 
changes  for  times  greater  than  the  inverse  of  its  resolution;  i.e.,  Te/TBW, 
where  TBW  is  the  time-bandwidth  product  of  the  spectrum  analyzer.  The 
photodiode  bandwidth  which  will  allow  the  detector  to  resolve  the  above 
changes  is 


Af  = 


(37) 


For  Te  = 3 psec  (corresponding  to  a 1 cm  optical  aperature)  and  TBW  variable, 
we  have  evaluated  Eq.  (36)  and  the  results  are  shown  in  Table  II. 
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TABLE  II.  LASER  POWER  REQUIREMENTS-PHOTODIODE  DETECTION 


TWB 

(Time-Bandwidth  Product) 
and  the  Number  of 
Spectral  Points 


Required  Laser 
Power  (mW) 


32 

64 

128 

256 


1.8 

3.6 

7.2 

14.4 


For  CCD  detection  the  number  of  carriers  excited  at  a detector  element 
corresponding  to  the  incident  power  in  Eq.  (33)  and  an  integration  time  of  Tx 
seconds  is 


(38) 


where  hv  is  the  incident  photon  energy  and  nq  is  the  detector  quantum  efficien- 
cy. The  integration  time  Tx  corresponds  to  the  laser  pulse  duration  for  Case  I 
operation,  to  the  transit  time  Te  for  Case  II  operation,  and  to  the  time  T0 
required  to  empty  a single  output  CCD  for  Case  III  operation.  The  incident 
deflected  power  Pjnc  corresponds  to  the  maximum  optical  power  which  woul-1  ever 
reach  a detector  element  so  the  number  of  excited  carriers  n should  be  equal 
to  the  maximum  number  which  the  detector  can  handle;  i.e.,  nmax  in  Eqs.  (13) 
and  (16)  in  order  to  fully  utilize  detector  dynamic  range.  Thus,  imposing 
the  condition 


n " nmax  » 


(39) 


in  Eq.  (38)  and  using  Eq.  (33) 
laser  power  P0  to  be 

a,Li  a L„ 

hvn  e 1 e 9 9 
p _ max 

° " VdmaxnslnlcrtTx 

We  have  evaluated  Eq.  (40)  for 
nmax  9iven  in  Eq.  (16)  and  N = 
is 


in  Eq.  (38),  we  can  solve  for  the  required 


(40) 

the  above  parameters  along  with  f c = 1 0 MHz, 

TBW  variable.  The  variable  integration  time  Tx 


215 


( .2  ysec  Case  I 


Tx  = < 3.2  ysec  Case  II  (41) 

( 10'1  Nusec  Case  III 

The  results  of  the  calculation  of  required  laser  power  using  Eq.  (40)  2**e 
tabulated  in  Table  III  for  various  values  of  N. 


TABLE  III.  LASER  POWER  REQUIREMENTS  - CCD  DETECTION 
P0  = Laser  Power  Required  (mW) 


N 

Number  of 
Spectral  Points 

Case  I 

Pulsed  Mode 

Pulse  Width=.2usec 

Case  II 

CW  Mode 

Integration  Time=Elastic 
Wave  Transit  Time 

Case  III 

CW  Mode 

Integration  Time= 
Time  Required 

To  Empty  CCD 

128 

73.6 

4.6 

1.15 

256 

73  5 

4.6 

.58 

512 

73.6 

4.6 

.29 

1024 

73.6 

4.6 

.14 
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Figure  1.  Integrated  Optical  Waveguide  and  CCD  Detector  Array 
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